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ABSTRACT 



Magnetic reconnection in plasmas, a process characterized by a change in connec- 
tivity of field lines that are broken and connected to other ones with difi^erent topology, 
owes its usefulness to its ability to unify a wide range of phenomena within a single 
universal principle. There are newly observed phenomena in the solar corona that 
cannot be reconciled with two-dimensional or steady-state standard models of mag- 
netic reconnection. Supra-arcade downfiows (SADs) and supra-arcade downfiowing 
loops (SADLs) descending from reconnection regions toward solar post-fiare arcades 
seem to be two different observational signatures of retracting, isolated reconnected 
fiux tubes with irreducible three-dimensional geometries. This dissertation describes 
work in refining and improving a novel model of patchy reconnection, where only a 
small bundle of field lines is reconnected across a current sheet (magnetic discontinu- 
ity) and forms a reconnected thin fiux tube. Traditional models have not been able 
to explain why some of the observed SADs appear to be hot and relatively devoid of 
plasma. The present work shows that plasma depletion naturally occurs in fiux tubes 
that are reconnected across nonuniform current sheets and slide trough regions of 
decreasing magnetic field magnitude. Moreover, through a detailed theoretical anal- 
ysis of generalized thin fiux tube equations, we show that the addition to the model 
of pressure-driven parallel dynamics, as well as temperature-dependent, anisotropic 
viscosity and thermal conductivity is essential for self-consistently producing gas- 
dynamic shocks inside reconnected tubes that heat and compress plasma to observed 
temperatures and densities. The shock thickness can be as long as the entire tube and 
heat can be conducted along tube's legs, possibly driving chromospheric evaporation. 
We developed a computer program that solves numerically the thin fiux tube equa- 
tions that govern the retraction of reconnected tubes. Simulations carried out with 
this program corroborate our theoretical predictions. A comparison of these simu- 
lations with fully three-dimensional magnetohydrodynamic simulations is presented 
to assess the validity of the thin fiux tube model. We also present an observational 
method based on total emission measure and mean temperature to determine where 
in the current sheet a tube was reconnected. 
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1. INTRODUCTION 

1.1. The Nearest Star 

The Sun sustains life on Earth. We are connected to it in an essential way. It 
powers our world, warms it to a temperature that allows liquid water to be, and 
drives our weather, evaporating water from the seas and producing winds. Rain and 
snow irrigate our lands that together with sunlight allow vegetation to grow. Plants 
generate oxygen that we breathe and feed the animals that we eat. Solar power, 
directly or indirectly, fuels machines, energizes the lights of our homes, and powers 
our cars. It is an amazing cycle. 

We move trough space at the right distance from the Sun for life, the way we 
know it, to thrive. Had the Earth been closer, we might have been a scorching Venus 
inferno; farther away, a frozen Mars. If there is life somewhere else, we do not know. 
We know there is life here, and other planets in our solar system do not seem to 
support it; we might be inclined to assume that habitable conditions are rare. But 
we only have one data point, and it is overly optimistic to think that we completely 
understand the circumstances that make life possible. Is planet Earth unique? I 
doubt it. There are so many planets in the universe, that the likelihood of finding 
Earth-like planets seems non trivial. What we definitely know is that the Sun plays a 
determining factor for life. Other stars may do the same in their planetary systems. 
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From the universe's viewpoint, our Sun is a fairly ordinary star. Not the brightest. 
Not the biggest. It is not especially variable or active star, it does not have a peculiar 
magnetic activity or age. Our solar system is located in a spiral arm of the Milky way 
galaxy, it is one more among the other stars and planets. Figure 1 shows an artist's 
concept of our galaxy and the location of our Sun. There are countless galaxies; 
we are only a speck in the vastness of the cosmos. In addition, the entire human 
existence is much less than a sigh in the universe's history. Setting this existential 
lightness aside, it is prudent to understand as much as we can from the star that is 
the foundation of our existence. We would like to know how it works, why it changes, 
and how these changes influence us here. 

The Sun is also the key to understand other stars, which are also huge incandes- 
cent balls of plasma. Many physical processes that occur elsewhere in the universe 
can be examined in detail on the Sun. We know its age, radius, mass, and lu- 
minosity (Priest, 1982). Through helioseismology we have also been able to learn 
detailed information about its interior (for a thorough review on helioseismology, see 
(Christensen-Dalsgaard, 2002)). We now know that we live in the extended solar 
atmosphere. This outer part of our nearest star is called solar corona and is perme- 
ated by magnetic fields generated inside the Sun that extend through the entire solar 
system. The Sun's visible sharp edge, the photosphere, is something of an illusion. 
It is the level beyond which the gas in the solar atmosphere is thin enough to be 
transparent. 



Figure 1. Artist's concept illustrating a view of the Milky Way. Our sun lies near a 
small, partial arm called the Orion Arm, or Orion Spur, located between the Sagit- 
tarius and Perseus arms. Image credit: NASA/JPL-Caltech. 

Our star blows itself away sending more than a million tons of electrons and 
protons (Feldman et al., 1998) into space every second (Withbroe & Noyes, 1977). 
This seemingly eternal and stormy solar wind, first predicted by Parker (1958), which 
streams past the Earth at averaged speeds of more than 400 km s"^ (Gosling et al., 
1976), fills the entire solar system. The solar wind was first observed in 1959 by the 
Lunik III and Venus I satellites (Priest, 1982), began to be studied in detail by the 
Mariner II probe in 1962, and has been studied more recently with Ulysses spacecraft 
(Aschwanden, 2005). Disturbances in the solar wind shake the Earth's magnetic field 
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and pump energy into the radiation belts (van Allen, 1981). These belts have harmful 
effects on satellites and humans working in space (Schrijver, C. J. & Siscoe, G. L., 
2010). 

The Sun's proximity allows us to study it with astonishing resolution. Fortunately, 
this not only makes it a physical laboratory where plasma theories can be tested, but 
also makes it a remarkable source of newly observed phenomena. New models are 
needed to explain observed features that cannot be explained with standard theories. 
This dissertation is an attempt to unify under the umbrella of a novel model some of 
such observations. 

1.2. The Biggest Explosions in the Solar System 

The Sun generates the biggest explosions in the solar system. Enormous eruptions 
of solar coronal material known as coronal mass ejections (CMEs) send clouds of hot, 
magnetized plasma out into planetary space (reviews on CMEs can be found, for 
example, in Forbes (2000) and Klimchuk (2001)). They apparently result from a 
rapid, large-scale restructuring of magnetic fields in the low solar corona. After a few 
solar radii, they have reached a constant velocity between 400-600 km s~^ (gradual 
CMEs, Sheeley et al. (1999)), and may exceed 750kms~^ (impulsive CMEs, Sheeley 
et al. (1999)). If directed toward Earth, CMEs typically reach us one to five days 
after their eruption. They interact with Earth's magnetosphere and may produce 
powerful geomagnetic storms, intense aurorae, and electric power blackouts. They 
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shake the Earth's magnetic field causing current surges in power fines tfiat destroy 
equipment and knock out power over large areas. Most of tfie energy released by 
tfie eruption goes to tfie expelled material, tfiis translates in an average of 10'^° to 
10^^ erg of kinetic energy (Vourlidas et al., 2000). Energetic protons accelerated by 
intense solar explosions can inflict serious damage on any astronaut caugfit in space 
witfiout adequate sfiielding (National Researcfi Council, 1996). Figure 2 sfiows a 
CME erupting from tfie Sun as observed by tfie LASCO C3 coronagrapfi aboard tfie 
Solar & Heliospfieric Observatory (SOHO) spacecraft (Domingo et al., 1995). 

Anotfier dramatic, magnetically energized type of solar explosion are solar flares. 
Tfiey are sudden outbursts tfiat rip tfirougfi tfie atmospfiere above sunspots, releas- 
ing unimaginable amounts of energy, fieating solar corona plasma to temperatures 
between 10-20 MK (Ascfiwanden, 1999), and sometimes exceeding 40 MK (Lin et al., 
1981; Holman et al., 2003; Caspi & Lin, 2010). In just 100 to 1000 seconds, tfie blast 
can release 10'^^ erg in a relatively compact region (Ascfiwanden, 2005; Crosby et al., 
1993). Figure 3 sfiows an image in 195 A of a flare observed by tfie Transition Region 
and Coronal Explorer (TRACE) (Handy et al., 1999). 

In solar flares, most of tfie magnetic energy is mainly transferred into accelerated 
particles tfiat subsequently emit intense X-ray, 7-ray, and radio radiation, but tfie 
exact acceleration mecfianisms befiind tfiis process remain poorly understood (Lin & 
Hudson, 1971, 1976; Emslie et al., 2004, 2005; Caspi, 2010). Tfiis radiation from solar 
flares may fieat up tfie Eartfi's upper atmospfiere, cfianging its electrical properties. 
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Figure 2. Coronal mass ejection, planets, and Pleiades. A Coronal Mass Ejection 
erupts from the Sun whose size and position are indicated by the white circle. This 
image was captured by LASCO C3 coronagraph aboard SOHO, on May 15, 2000. 
LASCO is able to take images of the solar corona by blocking the light coming directly 
from the Sun with an occulter disk, creating an artificial eclipse within the instrument 
itself. The shadow crossing from the lower left corner to the center of the image is 
the support for the occulter disk. The field of view encompass approximately 31 
diameters of the Sun. A rare alignment of four planets is also shown; the bright lines 
bisecting their disks is an instrumental effect due to the intense sunlight reflected 
from the planets. The Pleiades star cluster is 408 light-years away from the Sun. 
Image credit: Courtesy of SOHO/LASCO C3 consortium. SOHO is a project of 
international cooperation between ESA and NASA. 

disrupting radio navigation or telecommunications, and making it expand farther into 
space. Friction between the expanded atmosphere and satellites may shorten mission 
lifetimes. As we become more dependent upon satellites in space, electronics, and 
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telecommunications, we will increasingly feel the effects of this space weather. The 
need to understand and predict flares, as well as CMEs is clear. 




Figure 3. Flare observed by TRACE. On 21 April 2002, TRACE observed an XI. 5 
flare in Active Region 9906 at the solar limb (circular horizon). The image is rotated 
so that north is to the left. The observations, at a high cadence in 195 A (a mixture of 
Fe XII at 1.5 million degrees and Fe XXIV at some 10 million degrees) show localized 
dark blobs that move downward into the bright haze (one of them is shown with a 
white arrow). An arcade of cooling loops form, with dark, cooled material draining 
back towards the solar surface. Earth size is shown to scale. Acknowledgement: 
TRACE is a mission of the Stanford-Lockheed Institute for Space Research, and part 
of the NASA Small Explorer program. 
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1.3. Magnetic Reconnection 

The only plausible source of energy for the colossal explosions mentioned in the 
previous section is the strong magnetic field in the solar corona (Forbes, 2000). Even 
though the plasma in the lower solar corona is hot (approximately 1 MK), the mag- 
netic pressure there greatly exceeds the gas pressure. Solar flares occur in active 
regions associated to sunspots, where magnetic fields are clumped together in intense 
bundles. Flares and CMEs are synchronized with an 11-year cycle (Hoyt & Schatten, 
1998a,b) of solar magnetic activity. The number of active regions, with their ener- 
getic and magnetized loops, varies from a minimum to a maximum and back to a 
minimum with a period of approximately 11 years. The frequency of solar explosions 
are greatest at the maximum of the cycle and lowest at the minimum of the cycle 
(Norquist, 2011). 

The coronal plasma can generally be approximated as a perfectly-conducting fluid 
(Priest, 1982; Low, 2003) which thus obeys the principles of ideal magnetohydro dy- 
namics (MHD). This zero-resistivity approximation immediately leads to the"frozen- 
in" condition that guarantees that fluid particles and field lines move together (Priest, 
1982, 2000; Birn & Priest, 2007). Therefore, plasma is confined by the magnetic field 
and shaped mostly into coronal loops (like the ones shown in Figure 3). 

The magnetic energy necessary to power solar explosions is stored in the form 
of non-potential magnetic field energy (Priest & Raadu, 1975; Hu & Low, 1982; 
Aly & Amari, 1989; Birn & Priest, 2007). An equilibrium field with no currents 
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except for a current sheet (magnetic discontinuity) is the lowest energy state (for 
perfectly conducting evolution) subject to a photospheric flux constrain (Longcope, 
2001), therefore they are stable under ideal MHD evolution. The current sheet raises 
the field energy with respect to a potential field with the same photospheric flux. If 
the perfectly conducting assumption breaks, the field reorganizes itself toward lower 
energy states. This restructuring of the magnetic field, called magnetic reconnection, 
is characterized by a change in connectivity of field lines that are broken and con- 
nected to other ones with different topology. In this way, large amounts of stored 
magnetic energy can be released as thermal, kinetic, and radiative energy, as well as 
accelerated particles (Mandrini, 2010). 

Magnetic reconnection owes its appeal to its ability to unify a wide range of 
phenomena within a single universal principle (Forbes & Isenberg, 1991; Priest, 2000). 
It is not only present in many theories that undertake to explain the initiation and 
evolution of flares or CMEs (Lin & Forbes 2000; Forbes & Acton 1996; Forbes et al. 
2006; Schrijver 2009; see also Forbes 2000 for energy estimations in CMEs), but it is 
also involved in the interaction of the solar wind with Earth's magnetosphere, where 
substorms originate (Pudovkin & Semenov, 1985). (For a recent review on magnetic 
reconnection in space plasmas and laboratory, see Yamada et al. (2010).) 

Some kind of diffusion mechanism initiates the connection of field lines from one 
side of a current sheet to field lines on the other side. This transfer of flux from one 
domain to another decreases the energy of the field. The exact nature of this process 
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is not yet well understood, but it has recently become clear that it must produce an 
electric field restricted to a small portion of the current sheet (Biskamp & Schwarz, 
2001; Birn et al., 2001). 

A schematic fiare-like scenario is depicted in Figure 4, where an arcade similar 
to the one in Figure 3 is represented by gray arc-shaped loops. Reconnection occurs 
within the current sheet above the arcade of post-fiare loops, represented in the figure 
by a solid rectangle. In this figure, skewed magnetic fields are assumed on each side 
of the current sheet. Field lines on the side of the current sheet that is closer to the 
viewer (solid black lines) form a finite angle with field lines on the back side of the 
current sheet (gray solid lines with endpoints in the top and bottom walls of the box). 

Magnetic reconnection was first quantitatively modeled in two dimensions and 
steady-state configurations with uniform resistivity (Sweet, 1958; Parker, 1957), but 
the energy release in this model was too slow to explain the rapid development ob- 
served in fiares and CMEs (Parker, 1963). This problem was solved by assuming a 
two dimensional steady electric field localized to a small diffusion region with four 
slow-mode shocks (SMSs) attached to it. At these shocks, the magnetic field is de- 
fiected and partially annihilated, and the plasma is heated and accelerated (Petschek, 
1964). The compressible plasma case was studied by Soward & Priest (1982). 

Localized, non-steady models began with the work of Semenov et al. (1983), fol- 
lowed by Biernat et al. (1987), and Nitta et al. (2002) all of which were two dimen- 
sional. For these spatially localized (in only one dimension) and short-lived episodes. 
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Figure 4. Flare current sheet cartoon. The sohd rectangle represents the current sheet 
plane located above the post-flare arcade (gray arc-shaped loops at the bottom). The 
black lines depict some magnetic fleld lines on the side of the current sheet that is 
closer to the viewer, and the dark gray ones correspond to fleld lines on the back side 
of the current sheet. The lower plane of the box represents the solar surface. The 
small sphere in the current sheet shows a patchy reconnection region. Only a small 
bundle of fleld lines that intersect this region are going to reconnect. The dashed 
rectangle in the plane of the current sheet corresponds to the region shown in Figure 
9. 
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four slow-mode shocks (SMSs) extend from the reconnection site as well, but they 
close back together in a finite region that increases in size as time goes by (as opposed 
to the infinite length SMSs from steady state). These teardrop-shaped SMSs continue 
growing and retracting as they sweep up mass, even if reconnection ceases (Semenov 
et al., 1998). Therefore, this conversion of magnetic energy into kinetic and thermal 
energy continues and the dissipated energy can exceed the ohmic dissipation from the 
reconnection electric field, indicating that a small scale, short lived event can affect 
the global structure of magnetic fields. 

All of the models above considered reconnection between magnetic field lines 
assumed to be perfectly antiparallel. A departure from this special case, called recon- 
nection between skewed fields, introduces a magnetic field component in the ignor- 
able direction, called a guide field; such models are sometimes called 2. 5- dimensional 
(Figure 4 is an example of this kind of configuration). In this case, it is necessary to 
supplement the SMSs with rotational discontinuities (RDs; Petschek & Thorne, 1967; 
Soward, 1982; Heyn et al., 1988; Skender et al., 2003), as shown in Figure 5, from 
Longcope & Bradshaw (2010). In the special case of uniform anti-parallel fields, these 
two distinct shocks merge into a single SMS of the "switch-off type" (y-component 
of the magnetic field is annihilated), often associated with the Petschek model. The 
effect of the RDs is to produce bulk flow in the ignorable direction in order that the 
velocity change, magnetic fields, and shock normal may all lie in a single plane; this is 
a requirement of the SMS often called "co-planarity" (Colburn & Sonett, 1966). This 
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Figure 5. Petschek reconnection of skewed magnetic fields. Reconnection occurs 
within a portion of a large-scale current sheet (far left). Two perspectives of the 
reconnection region are shown: the end-on view (center, whose horizontal scale is 
exaggerated for clarity) and the face-on view (right). Flux transfer occurs along a 
line X (gray) along the ignorable direction (z). Shocks, RD', SS' (SMS), SS (SMS), 
and RD (thick dashed lines) originate there, dividing the reconnection region into 
sub-regions labeled 0', 1', 2, 1, and (inside boxes). The central region, 2 (shaded 
gray) contains plasma which has been heated and compressed by the SSs (SMSs). 
Reconnected field lines (thin solid lines) bend at each shock. A single flux tube is 
called out in gray in each panel. It passes though the shocks at points R' , S' , S, and 
R, left to right. Credit: Longcope & Bradshaw (2010). 



new flow component is confined to the region between the RD and SMS. Non-steady 
reconnection between skewed magnetic fields was analyzed by Semenov et al. (1992) 
and Biernat et al. (1998) for dayside magnetospheric events. 

Common to all the reconnection models is the presence of Alfvenic outflows orig- 
inating at the reconnection site. In steady state models, these form coherent jets. In 
transient reconnection models, any outflow is more accurately described as a distinct 
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retraction event rather than as a steady jet. It is the RDs rather than the SMSs 
which accelerate the reconnection flux to the Alfven speed in an outflow jet. 

While all of the foregoing reconnection models predict shocks, few investigations 
of them have explicitly included viscosity or thermal conductivity. This omission is 
noteworthy since strictly speaking the irreversible nature of any shock demands a non- 
ideal effect; in the coUisional plasma of the solar corona this is usually assumed to be 
viscosity. For strongly magnetized plasma, viscous diffusion and thermal conduction 
are highly anisotropic. They are directed along the magnetic field since perpendicular 
coefficients are many orders of magnitude smaller than their parallel counterparts 
(Braginskii, 1965). Such strong anisotropy poses challenges for models, especially for 
numerical solutions. In addition, transport coefficients in a fully ionized plasma have 
a strong dependence on temperature (~ T^) making transport far more significant 
at high post-shock temperatures than it is assumed to be in general (Spitzer, 1962). 
It has therefore been recognized as a potentially significant factor in flares (Cargill 
et al., 1995; Forbes & Malherbe, 1986). In particular, the heat generated at SMSs is 
predicted to be conducted along reconnected field lines ahead of other reconnection 
effects, possibly driving chromospheric evaporation (Forbes et al., 1989). 

Nevertheless, reconnection models generally neglect viscosity altogether to say 
nothing of accounting for its strong temperature dependence and directional anisotropy. 
Craig et al. (2005) and Litvinenko (2005) studied the effects of temperature-independent 
viscosity in various reconnection disturbances. Only recently, Craig (2008) and Craig 
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& Litvinenko (2009) included anisotropic and temperature-dependent viscous dissi- 
pation in a three dimensional steady-state model at X-points. Yokoyama & Shibata 
(1997, 2001) included anisotropic thermal conduction in two dimensional simulations 
of reconnection, and Seaton & Forbes (2009) presented a detailed theoretical analysis 
of these simulations. In most of these two-dimensional, anti-parallel models, the heat 
conduction front extends in front of SMSs, as expected. 

Direct observation of the actual reconnection process has been elusive. Its inter- 
nal structure remains unknown, even for well-diagnosed laboratory experiments (Katz 
et al., 2010). In order to explain the short times involved in the mentioned astro- 
physical explosions, fast reconnection (Petschek, 1964) ought to occur in a minuscule 
diffusion region on a CS where non-ideal field-line transport coefficients are locally 
enhanced. This requires measurements at prohibitively high-resolution. Several fast 
reconnection models include localized magnetic resistivity enhancement on the CS 
(Ugai & Tsuda, 1977; Scholer & Roth, 1987; Erkaev et al., 2000; Ma & Bhattachar- 
jee, 2001; Biskamp & Schwarz, 2001). 

Only in the last 15 years or so has technology started to close the gap between the- 
ory and observations. Nonetheless, only indirect measurements of solar corona recon- 
nection have been possible. For example, the cusp-shaped, hot and dense flare loops 
(Tsuneta et al., 1992; Forbes & Acton, 1996) consistent with the classical CSHKP 
model (Carmichael 1964; Sturrock 1968; Hirayama 1974; Kopp & Pneuman 1976; 
see Figures 1 and 2 in Forbes & Acton 1996 for a modern version including thermal 
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conduction) are considered observational evidence of reconnection. These arcades, 
slowly growing and whose footpoints separate as their heights rise, are interpreted as 
closed magnetic loops formed after two-dimensional, steady-state reconnection. They 
pile up on top of each other to form the arcades. The earlier loops cool down as new 
loops lie on top of them, and an apparent footpoint motion is observed. The loops 
are hot, as expected from the conversion of magnetic energy into thermal energy from 
reconnection. 

There are several other indirect observations of reconnection, like reconnection in- 
flows (Yokoyama et al., 2001; Narukage & Shibata, 2006; Lin et al., 2005) with speeds 
between a few km s"^ to approximately 106 km s"^, as well as Alfvenic (460 — 3500 
km s"^) reconnection outflows (Lin et al., 2005; Wang et al., 2007) and reconnection 
in/out pairs (Sheeley & Wang, 2007; Savage et al., 2010). Loop shrinkage also seems 
to indicate reconnection. Forbes & Acton (1996), using images of flare loops taken 
with SXT (Tsuneta et al., 1991) aboard Yohkoh (Ogawara et al., 1991), showed re- 
connected cusped loops shrinking to form post-flare arcades. This process takes hours 
and is different from the fast passage (1-2 minutes) of field lines through the re- 
connection outflows (jets). Hard X-ray (HXR) sources near the apex of soft X-ray 
(SXR) loops indicate some high-energy process, possibly electron acceleration in that 
region (Masuda et al., 1994; Krucker & Lin, 2008; Tomczak, 2009). The acceleration 
may occur at the reconnection site or might be the result of the interaction of the 
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downward reconnection outflow jet with the top of the SXR arcade, producing su- 
perhot plasmas (Shibata et al., 1995). Plasmoid (magnetic island) formation (Forbes 
& Priest, 1983), as a consequence of reconnection, has been widely observed (e.g., 
Shibata et al. (1995); Nishizuka et al. (2010)). For a comparative study of plasmoid 
in reconnecting CSs (solar and terrestrial contexts), see Lin et al. (2008). 

1.4. Patchy Reconnection 

Even though most of the observations mentioned in Section 1.3 can be under- 
stood by standard, steady-state, and two-dimensional reconnection theories, there 
are newly observed phenomena that cannot be reconciled with these representa- 
tions. Supra-arcade downflows (SADs) and supra-arcade downflowing loops (SADLs) 
(Savage et al., 2010; Savage & McKenzie, 2011) seem to be two different observa- 
tional signatures of retracting, isolated reconnected flux tubes with irreducibly three- 
dimensional geometries. 

SADs, dark voids descending from reconnection regions, were first reported by 
(McKenzie & Hudson, 1999) and since then have been observed by several instruments 
(McKenzie, 2000; McKenzie & Hudson, 2001; Innes et al., 2003b,a; Asai et al., 2004; 
McKenzie & Savage, 2009; Savage et al, 2010; Savage & McKenzie, 2011). They 
are directly related to the path of bright loops that appear lower down in the corona 
(Sheeley et al., 2004). These downward- moving features are believed to represent flux 
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tubes, created by patchy reconnection, retracting downward to later form the post- 
eruption arcade (McKenzie & Savage, 2009). They descend through a supra-arcade 
fan (see Svestka et al. 1998 for a description of one of such fans) toward an underlying 
flare arcade and decelerate close to the arcade apex (Sheeley & Wang, 2002; Sheeley 
et al., 2004). The white arrow in Figure 3 points to one of these dark voids. 

SADs are spatially localized as dark pockets, followed by a dark wake or lane. 
These downflows are observed as high as 40 — 60 Mm above the top of soft X-ray ar- 
cades and interpreted as magnetic reconnection localized outflows (McKenzie, 2000). 
Therefore, SADs appear to be the result of intermittent or patchy reconnection high 
in the corona. Three-dimensional and time-dependent models are needed to explain 
them. Figure 6 shows examples of SADs, as observed by TRACE (Handy et al., 1999; 
Schrijver et al., 1999; Golub et al., 1999) in 195 A passband. The image to the left 
of the figure shows an arcade for the 2002 April 21 X-class flare (Wang et al., 2002; 
Gallagher et al., 2002; Innes et al., 2003a), where two SADs are enclosed by a white 
dashed rectangle. Usually, arcades like this one can be described as two dimensional 
or 2. 5- dimensional, suggesting the presence of a wide spread CS above them (as shown 
in Figure 4). The image to the right is a time sequence for the location delimited by 
the dashed rectangle in the left panel. The white arrow on the top right snapshot 
indicates a dark pocket (SAD's head) that moves toward the arcade. This pocket is 
followed by a dark wake. The white arrow at time 1:44:53 indicates a second dark 
pocket (SAD) moving in the same direction. 
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Figure 6. 2002 April 21 flare. Left panel: rotated TRACE 195 A image, taken at a 
time when dark voids can be seen descending toward the flare arcade. The dashed 
rectangle encloses a handful of them. Right panel: time sequence for the window 
delineated by the white dashed rectangle in the left panel. The white arrow in the 
1:43:33 snapshot shows the position of one of the dark pockets. The other arrow 
points to a second localized void. 

A three-dimensional model of patchy reconnection (Klimchuk, 1997), accompa- 
nied by magnetohydrodynamic (MHD) simulations, has been presented by Linton 
& Longcope (2006). They studied the relaxation of solar-like low plasma-/? (ratio 
between thermal and magnetic pressure) reconnected flux tubes. Reconnection is as- 
sumed to occur across a uniform CS with skewed magnetic fields. Even though the 
initial configuration is 2.5-dimensional, reconnection is initiated within a small patch. 
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producing two distinct bundles of reconnected field lines, flux tubes; the ensuing 
reconnection dynamics is irreducibly three-dimensional. 

In that study, resistivity is briefly enhanced by an unspecified mechanism within 
a localized and short-lived region (an illustrative reconnection patch is depicted in 
Figure 4 as a small sphere). This creates a pair of bent, thin reconnected flux tubes 
that retract in opposite directions. They found that the perpendicular dynamics of a 
reconnected tube was well described by equations for thin flux tubes (TFT; Spruit, 
1981). The equations differed from those used traditionally in convection zone models, 
by their application to the low-/3 corona. RDs propagate along the legs of the tubes. 
Only the perpendicular dynamics of the tubes was analyzed by Linton & Longcope 
(2006), who also argued that the descending coronal voids are reconnected flux tubes 
descending from a flare site high in the corona. 

Longcope et al. (2009) extended that study to include the tube's parallel dynam- 
ics. The initial sharp angle at the reconnection site creates bends that move along 
the legs of the tubes as RDs. These bends rotate the magnetic field and generate 
parallel supersonic inflows toward the center of the tube. There, the collision of these 
inflows launches outwardly two gas-dynamic shocks (GDSs) that heat and compress 
the plasma without changing the magnitude of the magnetic field. The GDSs are dis- 
connected from the resistive region and are features of the ideal relaxation following 
reconnection. The rate of reconnection is not relevant in this case since the diffusion 
region is short-lived and does not participate in the post-reconnection evolution. 



21 

The GDSs are parallel shocks and close relatives to the nearly parallel SMS in 2.5- 
dimensional steady models. In fact, the post-GDS temperatures in the TFT model are 
almost identical to the post-SMS temperatures in steady-state magnetohydro dynamic 
MHD models. In both kinds of models, plasma heating occurs through a two-stage 
process notably different from commonly used switch-off models, applicable to purely 
anti-parallel reconnection. First, the RD converts magnetic energy to bulk kinetic 
energy in the form of Alfven-speed flow partly parallel to the field line. The parallel 
component of this flow, present in both, is then thermalized in a second shock of large 
hydrodynamic Mach number (see Figure 1 in Longcope & Bradshaw 2010 for a single 
flux tube as it passes through these shocks). 

Patchy reconnection seems to agree with some observations of reconnection. Long- 
cope et al. (2005), through EUV observations and modeling, analyzed an emerging 
active region in the vicinity of an existing one. They studied non-fiaring reconnection 
along a separator overlying the volume between them. They characterized 43 loops 
observed by TRACE that resulted from reconnection between the active regions, and 
assumed a correspondence between them and modeled flux tubes. The average diam- 
eter of the loops is 3.7 Mm with very small variation, suggesting a sporadic or patchy 
reconnection. The estimated flux of each one of them is 4 x 10^^ Mx for an assumed 
magnetic field of 37 G. Figure 7 reproduces Figure 19 of that paper and illustrates the 
reconnection model applied to the active regions. The three-dimensional CS is shown, 
as well as a pair of reconnected flux tubes created by patches of reconnection. The 
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Figure 7. Illustration in perspective of a model of reconnection. The model assumes 
reconnection between solar active regions 9570 and 9574. The vertical separator 
current sheet is shown as a skeleton, and a pair of reconnected flux tubes, 4 x 10^* Mx 
each, created by patches of reconnection on the separator current sheet, lie underneath 
it. Credit: Longcope et al. (2005). 



time history of the loops (observed at 171 A) with very similar diameters suggests a 
highly intermittent reconnection process. 

Although what triggers reconnection is still an open question, several promising 
mechanisms have been proposed. For example, the MHD CS tearing mode instability 
(Furth et al., 1963; Sturrock, 1966; Forbes & Priest, 1983; Priest, 1985; Biskamp, 
1986) may initiate bursty reconnection. Shimizu et al. (2009) showed numerically 
that a two-dimensional fast magnetic reconnection in a simple one-dimensional CS 
can be unstable to three-dimensional resistive perturbations. Therefore, it is possible 
to randomly eject three-dimensional magnetic loops along a CS atop a post-flare 
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arcade, resembling SADs. Catastrophe models explain how magnetic configurations 
are stable for a long time and suddenly release their energy (Cassak et al., 2005). 
The Hall reconnection model, extensively studied in the Geospace Environmental 
Modeling (GEM) magnetic reconnection challenge (Birn et al. 2001, and companion 
papers), yields energy release rates consistent with the ones observed in the solar 
corona (Drake et al., 2006). 

Even though three-dimensional outflows can be generated from purely two-dimensional 
conflgurations, the assumption of perfectly anti-parallel pre-reconnection magnetic 
flelds seems too restrictive for the complex solar corona, where fleld lines are likely 
to form flnite angles across CSs. Furthermore, the addition of a guide fleld is directly 
related to the temporal evolution of reconnection. Drake et al. (2006) showed that 
ambient guide magnetic fleld controls whether magnetic reconnection, once it is estab- 
lished, remains steady (anti-parallel reconnection) or becomes bursty (reconnection 
with a guide fleld) . Through particle simulations in a system consisting of two Harris 
CSs (Harris, 1962), they show that magnetic reconnection in a large-scale system with 
a guide fleld will break up into many islands, as a result of tearing instability. Guide 
flelds with a shear angle of 127° are sufficient to produce secondary islands that grow 
to finite amplitude, producing well-developed fiux tubes. 

Besides the assumption of being retracting reconnected fiux tubes, there are sev- 
eral unanswered questions about SADs. To the best of our knowledge, there has not 
been a satisfactory explanation of why they seem to be hot and relatively devoid of 
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plasma. McKenzie & Hudson (1999), using an image-ratio method between the SXT 
A112 and AlMg filters (Tsuneta et al., 1991), reported temperature ratios between 
the voids and the surrounding rays of the fan of 9.1 : 7.9. The wakes seem to have 
a slightly lower ratio: 8.6 : 7.9. The density ratios between the voids and the rays 
ranged between 1.0 : 4.0 to 1.3 : 2.3, depending on the assumptions about the line- 
of-sight depth. The darkness of the voids in X-ray and extreme-ultraviolet (EUV) 
images and spectra, together with the lack of absorption signatures in EUV (Innes 
et al., 2003b) is best explained as pockets of very low plasma density. Innes et al. 
(2003b) ruled out the possibility that the dark tracks are caused by coronal rain (cold 
plasma falling gravitationally). 

There is another piece of the puzzle that has not been completely understood: 
SADs' reported speeds (30 - 500 km s"^) are smaller than the assumed Alfven speed in 
the corona (McKenzie & Hudson, 1999; McKenzie, 2000; Asai et al, 2004; McKenzie 
& Savage, 2009; Savage et al., 2010). This contradicts any standard reconnection 
theory. While few if any of these events have measured speeds matching the local 
Alfven speed, they are generally interpreted as evidence of reconnection. Either these 
outflows/downfiows were Alfvenic at some earlier time or all reconnection theories 
must be revised. A possible explanation for reconnected tube's deceleration is given 
by Linton & Longcope (2006). They showed that a drag effect from the plasma around 
the tube can reduce its speed considerably. The external fluid deforms and accelerates 
as the tube moves through it; this requires energy and therefore slows the tube. The 
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ratio between the speed observed in their simulations and the perpendicular Alfven 
speed (expected speed from the TFT model) increases with the reconnection angle 
and ranges approximately between 0.1 and 0.5 for the low plasma-/? (ratio between 
thermal pressure and magnetic pressure) case. This "added mass" effect accounts for 
most of the flux tube drag at high plasma-/?, but is not enough, however, to explain 
the speed reduction at low plasma-/? (the expected speeds taking into account added 
mass effects are double the one observed in the simulations). They conclude that 
there is an additional source of drag in this case, but they do not speculate about its 
origin. 

1.5. Dissertation Layout 

This dissertation describes work done in refining the patchy reconnection model 
presented by Linton & Longcope (2006) and Longcope et al. (2009), in order to answer 
some of the questions stated above. 

Chapter 2 presents generalized TFT equations for the dynamics of reconnected 
flux tubes that include pressure-driven parallel dynamics as well as temperature- 
dependent, anisotropic viscosity and thermal conductivity. These additions are es- 
sential for self-consistently producing the GDSs predicted analytically by Longcope 
et al. (2009) using conservation laws. 

A model for plasma heating produced by time-dependent, spatially localized re- 
connection within a flare current sheet separating skewed uniform magnetic fields 
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(same reconnection scenario studied by Longcope et al. (2009)) is presented. The 
reconnection patch creates flux tubes of new connectivity which subsequently retract 
at Alfvenic speeds from the reconnection site. Heating occurs in gas-dynamic shocks 
which develop inside these tubes. The evolution of the reconnected tubes is studied 
through numerical solutions of the TFT equations, for realistic solar coronal condi- 
tions. 

The strong gas-dynamic shocks are generated by compressing plasma inside recon- 
nected flux tubes. They develop large velocity and temperature gradients along the 
tubes, rendering the diffusive processes dominant. These diffusive processes determine 
the thickness and internal structure of the shocks that evolve up to a steady-state 
value, although this condition may not be reached in the short times involved in a 
flare. This internal structure is analized using methods similar to those of Thomas 
(1944), Grad (1952), Gilbarg & Paolucci (1953), and Kennel (1988). For strong shocks 
at low Prandtl numbers (ratio of viscosity to thermal conductivity) , typical of the so- 
lar corona, the GDSs consist of a long thermal front where the temperature increases 
and most of the heating occurs, followed by a narrow isothermal sub-shock where 
plasma is compressed. The length of each of these sub-regions and the speed of their 
propagation is theoretically estimated. The thickness of the shocks is comparable to 
the entire length of the tubes, far greater than the term "internal structure" might 
suggest. The analysis is novel in its focus on solar flare conditions, including very low 
Prandtl numbers and transport coefficients with strong temperature-dependence. 
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Fowles (1975) showed that for strong shocks, it is impossible to assume inviscid 
behavior, no matter how large thermal conduction is. Coroniti (1970) studied dissipa- 
tion discontinuities in shock waves for temperature-independent transport coefficients. 
The seemingly harmless assumption of temperature-independent coefficients leads to 
shock thicknesses on the order of the particle's mean free path, raising concerns about 
the adequacy of a fluid treatment: when the thickness of the shock is of the order of 
particle mean free path, the fluid model breaks down (Campbell, 1984). However, if 
transport coefficients are given an accurate temperature-dependence, we have found 
that the thickness of the shocks exceeds the particle mean free paths by at least an 
order of magnitude. Our simulations also show that the ratio between the heat flux 
and the nominal free-streaming heat flux is smaller than 0.004. These two results 
allow us to treat our plasma collisional fluid. 

The work described in Chapter 2 has been published in the Astrophysical Journal 
(Guidoni & Longcope, 2010). 

In Chapter 3, through a simple but realistic fluid model of patchy reconnec- 
tion, it is shown that plasma depletion naturally occurs in flux tubes that are recon- 
nected across Syrovatskii-type (non uniform) CSs (Green, 1965; Syrovatskii, 1971) 
with skewed magnetic flelds (similar to the one shown on flgure 4). In this kind 
of CS, the background magnetic pressure has its maximum at the center of the CS 
plane and decreases toward its edges. The reconnection patch creates two V-shaped 
reconnected tubes that shorten as they retract in opposite directions, due to magnetic 
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tension. One of them moves upward toward the top edge of the CS, and the other 
one moves downward toward the top of the underlying arcade. 

Rotational discontinuities (RDs) propagate along the legs of the tubes and gener- 
ate parallel supersonic flows that collide at the center of the tube. There, gas-dynamic 
shocks that compress and heat the plasma are launched outwardly. A long lasting, 
dense central hot region created by the GDSs moves outward from the reconnection 
region. 

While the tube descends toward the top of the arcade, plasma density is decreased 
by tube's lateral expansions in response to the decreasing external background pres- 
sure. This effect may decrease plasma density by 30% - 50% of background levels. 
Therefore, retracting reconnected flux tubes may present elongated regions relatively 
devoid of plasma, resembling SADs. 

A reconnected tube will arrive at the top of the arcade that will slow it to a stop. 
Here, the perpendicular dynamics is halted, but the parallel dynamics continues along 
its legs; the RDs are shut down, and the gas is rarefied to even lower densities. After 
the tubes lie on top of the arcade, rarefaction waves continue decreasing plasma 
density. The hot post-shock regions continue evolving, determining a long lasting 
hot region on top of the arcade. We provide an observational method based on total 
emission measure and mean temperature, that indicates where in the CS the tube 
has been reconnected. 
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Current sheets themselves are notoriously difficult to observe directly, yet they 
are a necessary component of all magnetic reconnection models. Nevertheless, there 
are observations providing evidence (Schettino et al., 2010; Savage et al., 2010). Sy- 
rovatskii-type CSs may be ubiquitous in the solar corona. Edmondson et al. (2010) 
show that even a simple multipole topology would form very stable CSs bounded by 
two Y-type nulls, where an X-type null deforms under stresses due to motions in the 
photosphere. They did not include a guide field, but suggested that if added, the 
expected three-dimensional islands formed after reconnection may resemble longer 
fiux-tube-like magnetic islands. 

The work described in Chapter 3 has been published in the Astrophysical Journal 
(Guidoni & Longcope, 2011). 

In Chapter 4, we test the validity of the TFT assumptions by comparing results of 
the model with those from fully three-dimensional solutions of the general MHD equa- 
tions. In that chapter, we present three-dimensional MHD simulations carried out in 
a volume filled with plasma including a Green-Syrovatskii (Green, 1965; Syrovatskii, 
1971) CS. The reconnection episode is triggered by imposing a short-lived and local- 
ized magnetic resistivity patch on the top half of the current sheet. The results of 
these simulations are compared to the corresponding thin flux tube simulation. 

The work described in Chapter 4 will be submitted for publication in the near 
future. 
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As a final remark, we should mention that the correct answer to the disparity be- 
tween observations and models regarding the speed of localized reconnection outflows 
is still not known, but most likely the slowing down is due to the interaction between 
tubes and their surroundings. In the next chapters, for simplicity, this interaction 
is not included; the reconnected tubes are assumed to be completely isolated from 
their background. Therefore, we restrict ourselves to standard reconnection scenarios 
where the outflows are Alfvenic. 
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2. SHOCKS AND THERMAL CONDUCTION FRONTS IN RETRACTING 



The layout of the present chapter is as follows. The next section is a brief review 
of the general transport effects in magnetohydrodynamics. The following section 
describes the thin flux tube assumptions that lead to the final thin flux tube equations 
presented in the third section. Section 2.3 presents simulations of the retraction of a 
thin flux tube reconnected in a localized region inside a current sheet with uniform, 
skewed magnetic fields for a case with low plasma-/?, and parameters relevant to the 
solar corona. In the fifth section, an analytical analysis of the inner structure of the 
shocks is presented. Finally, in the last section, the main results of the paper and 
possible observational implications of the inner structure of the shocks are discussed. 



When resistivity, radiation and gravity are neglected, the magneto-hydrodynamic 
(MHD) equations for a charge-neutral plasma (we will assume sufficient coUisionality 
to justify their use) are (Shu, 1991, chap. 22) 



RECONNECTED FLUX TUBES 



2.1. MHD Equations 
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Here, p, P, B, v, T, s, ei, U, and q are the density, plasma pressure, magnetic field, 
velocity, temperature, specific entropy, specific internal energy, viscous stress tensor, 
viscous heat, and thermal conduction flux of the fluid, respectively. The advective 
derivative is defined as D/Dt = d/dt + v ■ V. In Equation (2.2), dyadic notation 
(Goldstein, 1959, chap. 5) is used for the second right-hand term. 

When the plasma is strongly magnetized by a field taken in the z direction, the 
stress tensor takes a simple gyrotropic form (Braginskii, 1965) 

n = -r^W^..(-^-^+i^), (2.5) 

where the zz component of the rate-of-stress tensor is 

dv 2 

W.. = 2^--V-v. (2.6) 

oz 3 

Here, r] is the largest of the dynamic viscosity coefficients, determined essentially by 
the ions (the electron-to-ion viscosity ratio is rje/rji ~ 0.01). The strongly magnetized 
limit of the conductive heat flux is 

q = -t,V\\{kBT), (2.7) 

where Vy = zz ■ V is the derivative parallel to the magnetic field, and k is the parallel 
thermal conduction coefficient, determined mostly by the electrons (the electron-to- 
ion heat conductivity ratio, Ke/ztj ~ 24). 
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We close the system using the ideal gas law 
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where 7 is the adiabatic gas constant, m is the average particle mass (for a fully 
ionized Hydrogen plasma, it is approximately equal to half the mass of the proton), 
n is the electron number density, cy = kB/{'y — l)rn is the specific heat at constant 
volume, and sq is a constant. 



2.2.1. The Thin Flux Tube Assumptions 

We consider an untwisted, thin tube of magnetic flux 0, like the one depicted in 
Figure 8, embedded in a strongly magnetized external fluid. The external magnetic 
field, thermal pressure, and density will be denoted Be, Pe, Pe, respectively. These 
quantities may depend on position, but not time, since we will assume the background 
configuration to be in equilibrium. The tube is considered so thin that the background 
quantities remain unperturbed as the tube moves through the external plasma. 

We will assume that the tube can be represented by an internal field line, R(/), the 
axis (see Figure 8). The field is untwisted in the sense that all field lines are parallel 
to the axis. The time it takes for a fast magneto-sonic wave to establish pressure 



2.2. Thin Flux Tube Equations 
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Figure 8. Frenet coordinate system. The axis of a thin magnetic flux tube is repre- 
sented by a central magnetic field line (dashed line) parameterized by vector R(/). 
The unit vector b is parallel to the local magnetic field, unit vector p points in the 
direction of the curvature vector k, and unit vector n completes the system. The two 
ellipses are the cups of a small tube segment of mass Sm. 

balance with its surroundings scales as Tjm = a(/)/ ^/c'^{l) + vl{l), where Va is the 
Alfven speed, and Cs is the sound speed. Since the fast-magnetosonic time is negligible 
with respect to the Alfven time over a characteristic scale Lch, ta = Lch/va{l), we 
assume instantaneous, local pressure balance 



R2 r2 
P+— = Pe + ^, 



(2.11) 
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where P and B are pressure and magnetic field inside the tube, respectively. 

It is convenient to define a Frenet coordinate about the axis R, parameterized by 
its arc-length /, as shown in Figure 8. The Frenet unit vectors defining the coordinate 
system are 

p(/) = k(/K(/) = ^r,(/), (2.13) 
n(/) = p(/) X b(/). (2.14) 

The unit vector b is parallel to the axis, p is the normal vector and n is the bi-normal 
completing the system. Local derivative operators associated with the unit vectors 
are 

d\\ = h-V = du 
V_L = -b X (b X V). 

We assume that the fiux tube is sufficiently isolated from its surrounding that no 
viscous stress crosses its surface. This is justified by the small cross- field transport 
coefficients. The stress tensor is therefore 



n = -^(b-5)bb. (2.15) 



dl ^ 

For parameters typical at the high corona, T = 1 MK, n = 10^ cm"^, i? = 10 G, 
L = 100 Mm, the Reynold's number is -R„ ~ 3 x 10'^, indicating a relatively small 
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contribution from viscosity. For such high Reynolds numbers, the viscous stress tensor 
may become significant only at steep velocity gradients, like shocks. Shocks in thin 
fiux tubes were predicted by (Longcope et al., 2009). 

The small Prandtl number of an ionized plasma, P,- = rj/nm ~ 0.01, indicates 
that the effect of thermal conduction is larger than viscosity. Temperature gradients 
will then be even smaller at shocks than velocity gradients. 

2.2.2. Thin Flux Tube Equations 

The thin flux tube equations are derived by applying the MHD equations to a 
small piece of the tube with mass (5m, like the one shown in Figure 8. The position of 
every point in the tube can be expressed as r' = R+r. Vector r belongs to the circular 
area perpendicular to the tube axis, and can be expressed as r = r sin 6' n+r cos 6' p = 
rr, where r and 9 are the polar coordinates associated to the n-p plane. 

The local radius of the tube, a(/), is small compared with scales of variation along 
the tube, permitting Taylor-expansion of any state variable / as follows 

fc=0 fc=0 i=0 ^ ^ ^ 

The second equality was obtained expanding the binomial term in the local polar co- 
ordinates. The shortness of the tube element allows us to assume that the piece's axis 
lays on a plane and its radius of curvature is constant. For this case, the differentials 
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of area and volume for the piece can be expressed as 
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An example of the expansion (2.16) is the mass of the tube piece. Computing the 
volume integral using expression (2.19) for the volume differential, and expanding the 
density around the center of the tube, gives 
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where (9(a (/')) represents an error of fourth order in the local radius. If the arc-length 
interval, 5/, is small then the above equation can be expressed as 



5m{l) ^ p{l)Tia^{l)5l. 



(2.20) 
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Expanding the dot product B ■ b, gives a similar expression for the magnetic flux 



B ■ rfAii 

I 

= B{l)Tia^{l) + 0{a\l))-B{l)na^{l). (2.21) 

For a given tube piece, the ratio between the mass and the magnetic flux, = 
(5m/</) = {p/B)5l, is constant. This deflnition, along with Equations (2.20) and (2.21), 
allow us to deflne the integrated mass per flux, /i = / {p/B)dl, to parameterize the 
tube instead of /. The differential with respect to fj, is 

± = .2 22) 

In order to obtain the thin flux tube mass equation we integrate over the piece's 
volume every term of the corresponding MHD equation, and keep only the flrst non- 
vanishing terms. To the flrst order, the velocity at any point in the tube is the sum of 
the velocity of the center of the tube plus a radial term Vr- This radial velocity at the 
surface of the tube is equal to the change in local radius {D / Dt)a{l) . Finally, applying 
the fundamental theorem of calculus to recover integrals along the arc-length, and 
expression (2.21), the flnal result is 

D fB\ ^ dv , , 

In the integration of the momentum Equation (2.2), the surface integrals for the 
perpendicular area vanish since B ■ dA^ = and 11 ■ dA^ = (in this last equality 
we have used deflnition (2.15)). For the flrst term of the right-hand side, constraint 
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(2.11) is applied and expanded around the center of the tube yielding 



Dt 




(2.24) 



where every term is evaluated at the center of the tube. In the above equation, the 
viscous momentum term is a perfect derivative in the integrated mass. There is no 
viscous momentum interchange between the tube and its surroundings. 

To determine the form of the viscous heating, we compute the rate of change in 
kinetic energy as the dot product between the momentum Equation (2.24) and the 
velocity. Balancing this loss with a heating rate gives 



With the above definition, and similar treatment than for the mass and momentum 
equation, the thin flux tube energy equation becomes 



where we see that the thermal conduction term is also a perfect derivative in the 
integrated mass per flux (there is no thermal conduction across the tube). 

For our simulations, it is more convenient to re-write the energy equation by 
defining a new variable Pc in the following way 




(2.25) 




(2.26) 




(2.27) 
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In an adiabatic case, Pc is constant in time. The change in plasma entropy with 
respect to the initial state (assumed to be in equilibrium) is directly related to Pc, 
where As = cy In [Pc(/W, t)/Pe(At)]- With definition (2.27), the energy equation is 
transformed to an equation for P^ 



(2.28) 



Dt ' ' ' \ p 

When density is known, the above equation allows us to update pressure, and use it 
in the momentum equation. 

Heating and cooling occur where there is a change in entropy. From Equations 
(2.3), (2.26), and (2.28), the volumetric heating rate can be written as 



Dt 7-1 VPe/ Dt 
dw\ ^ „ d f hi d 



In this equation, the viscosity term is always positive, and therefore always heats 
the plasma. On the other hand, the thermal conduction term can have either sign, 
depending on the temperature's arc- length second derivative. 

The transport coefficients in an ionized plasma depend on temperature as (Spitzer, 
1962) 

,5 



VcT-^, (2.30) 
KcTI, (2.31) 



where rjc and are constants. 
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For thin flux tubes, there is no need to solve the induction equation, since the 
magnitude of the magnetic field inside the tube is determined by constraint (2.11), 
and its direction is given by the position of the tube piece governed by the momentum 
equation (resistivity has been neglected, and therefore the magnetic field follows the 
fluid). 

It is possible to describe the TFT equations in dimensionless variables, by defining 
a characteristic magnetic field magnitude, Bq, density, po, and length, Lq. For ex- 
ample, these characteristic variables may represent the magnetic field magnitude and 
density at the reconnection point, x/j, and the length of a typical reconnected tube, 
respectively. Then, Bq = B^i^fj) and po = pei'x.R). With these values, we can con- 
struct an associated Alfven speed vao = Bq/ y/Anpo, and an associated temperature 

= v^Qin/kB, to non-dimensionalized velocities and temperatures. 

The new dimensionless variables are 

,,* , , Bp * p J.* 4- ^ AO 1* I ID* p 47r p* p 47r ti* T 

V* = ^, B* = B#, X* = k* = kLo, h* = b, 

VAO ' -Bo ' Lo^ ' 



_ r]* {T* f/^ , with r/* = 7]J^ , 

PoLqVao IC \ ) 1 IC I'- poLoVAO 

n*\5/2 



_ f.* (T^fr^ ^ith ^e = K, ""^fy^\ (2.32) 

poLoVAO C \ / ^ C L poLoVAO ^ ' 



With these new variables, the dimensionless thin flux tube equations have exactly the 
same form as Equations (2.23), (2.24), (2.26), and (2.28), except that the Boltzmann 
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constant is no longer present. From now on, we will refer to the thin flux tube 
equations as the dimensionless ones. 

2.3. Solutions of Reconnection Dynamics 

2.3.1. Low Plasma- /3 Uniform Background Case 

The results of the last section are quite general. In this section, we will simplify 
the equations to apply them to the low-/3 solar corona. For simplicity, we will also 
assume uniform background magnetic field. 

We will consider skewed, uniform external magnetic field forming an angle 2C, on 
opposite sides of a static, infinite current sheet, as shown in Figure 9. The current 
sheet is assumed to be infinitely thin, represented by the a; — |/-plane. The initial pres- 
sure is also assumed to be uniform everywhere. A spatially and temporally localized 
reconnection event has occurred (by an unspecified physical mechanism), connecting 
some field lines from one side of the current sheet to some on the other side. The 
amount of flux is determined by the size of the reconnection region. This bundle of 
field lines forms two V-shaped flux tubes (as shown in panel (b) of Figure 9) and 
magnetic tension at their cusps causes them to retract. These tubes will remain 
as coherent entities, preserving their magnetic flux </), while they move. No further 
reconnection is assumed. 

For this uniform case, the characteristic magnetic field magnitude and density at 
the reconnection point are the same as everywhere else. Then, Bq = B^, po = Pe, 
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Figure 9. Current sheet schematic, (a) The two tubes represent the front (sohd 
contour) and back (dashed contour) pre-reconnection tubes. Thick arrows indicate 
the direction of the magnetic field. The sphere in the middle represents the small 
reconnection region. The tubes are initially part of the uniform background equilib- 
rium, (b) Reconnected V-shaped tubes a short time after reconnection. The angle 
between the initial tubes is 2C,. Due to magnetic tension they will retract in opposite 
directions. The end-points of the flux tubes are located near the edges of the current 
sheet. 

and Lq is some characteristic length. The perpendicular gradient of the external total 
pressure is zero. 

In the solar corona, the ratio between thermal pressure and magnetic pressure, 
P = SttP/B"^, is small. For a plasma with magnetic field Bg = 10 G, number density 
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He = 10^ cm"^, and temperature Tg = 1 MK, /3 is approximately equal to 0.007. 
Consequently, we will neglect the components of the thermal pressure gradient which 
are perpendicular to the tube's axis with respect to the magnetic pressure gradient in 
that direction. The initial small internal pressure might increase due to compressions 
along the tube, but we will consider it small compared to the internal magnetic 
pressure at all times. In the parallel direction, there is no force related to the magnetic 
field, therefore the parallel derivative of the thermal pressure cannot be neglected. 

With the above considerations, and making the temperature dependence of the 
transport coefficients explicit, we can re- write the constraints and thin flux tube 
equations in the following way 
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where 



St — Si + Sk + Smi 



(2.37) 
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with specific internal, kinetic, and magnetic energies 

The specific magnetic energy of the tube doubles because there is work done by the 
background magnetic field expanding into volume vacated by motion of the tube. We 
have assumed that the external field returns to its original value after the tube has 
passed. In more realistic scenarios this value may be slightly different. 

Conservation of energy improves stability of any numerical solution. Toward this 
end, we will neglect the small last term in Equation (2.36) (proportional to /3), that 
is not part of the perfect derivative in the integrated mass. With this simplification. 
Equation (2.28) remains unchanged, and the change of the tube's total energy per 
fiux is equal to 



D (Et 



, (2.39) 

ML 



47r Be Be \ dl J 

where /xq and pi represent the end-points of the tube, respectively. If the-end points 
are fixed, and there is no temperature gradient there, the total energy of the tube 
is conserved. The neglected term came from the change in kinetic energy, calculated 
from Equation (2.35), and to be consistent, the correspondent perpendicular term 
will be neglected in the momentum equation, as well. Then, the change in total 
momentum per fiux of the tube is equal to 
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With the neglect of the proportional to /3 small term, and the fact that Pe is uniform, 
the ideal part of Equation (2.35) reduces to the thin flux tube equation presented by 
Longcope et al. (2009). Therefore, gas-dynamic shocks and rotational discontinuities 
are expected along the reconnected flux tubes. 

The volumetric heating rate when the magnetic field magnitude is uniform be- 
comes 



In the above equation, viscosity always contributes to heat the plasma. Depending on 
the curvature of temperature to the 7/2 power, and the value of the Prandtl number, 
the thermal conduction term may contribute to plasma cooling or heating. 

2.3.2. Simulations 

We developed a computer program, called Dynamical Evolution of Flux Tubes 
(DEFT), which solves the thin flux tube equations for the retraction of the two 
reconnected thin flux tubes described in the last section. The program takes into 
account the effect of the transport coefficients, including their strong dependence 
on temperature. It also can perform simulations using real coronal dimensionless 
parameters, as well as implement a strong anisotropy in the transport coefficients 
(heat conduction and momentum transport are allowed only in the direction parallel 
to the magnetic field). This is far more difficult to achieve in fully three dimensional 
MHD programs. The DEFT program can be implemented for different background 




(2.41) 
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conditions as well, but we will restrict ourselves in this paper to the low-/3 uniform 
background case from the last section. In this case, the two tubes moving in opposite 
directions are symmetric. In cases with non-uniform background it might not be so. 
We use the conservative form of the thin flux tube equations introduced in the last 
section. 

The DEFT program uses a staggered mesh where each tube piece is represented 
by grid points at its ends and by its mass per flux. The code implements a Lagrangian 
approach where each tube piece is followed, guaranteeing mass conservation. Due to 
mirror symmetry, only half of the tube is simulated, the other half is a mirror image 
in the x-direction. 

The initial state corresponds to two opposite V-shaped flux tubes, connected at 
the reconnection region, as shown in Figure 9. The initial angle between the magnetic 
field lines at each side of the current sheet is equal to 2C. Each tube is divided in N 
segments separated by iV -|- 1 tube-point positions. The initial tubes are clearly out of 
equilibrium due to the sharp angle between the field lines. To avoid introducing length 
scales at the limit of resolution, the central parts of the initial tubes are smoothed. 
The end-points of the tube are assumed to be fixed, and to have no temperature 
gradient with the external boundaries (this last condition ensures no heat transfer 
from the end-points). 

Since the reconnected tubes originated as background plasma they are initialized 
with uniform density, p = p^, and thermal pressure P = P^- The mass per flux of each 
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piece is equal to 6m/^ = pe6le/Be, where 6le is the segment's initial length. The mag- 
netic field magnitude, B^, is constant everywhere throughout the entire simulation. 
The initial velocities are zero (equilibrium state), even for the central bend (we are 
not interested in the details of the initial transient evolution). Viscosity and thermal 
conductivity are calculated using Equations (2.30) and (2.31), where T = P/p (unit- 
less quantities). The unit parallel vectors are calculated at segment's centers from 
definition (2.12), and curvature vectors are calculated at grid points from definition 
(2.13). Arc-length derivatives are computed using centered differences. 

At every time, the + 1 positions and velocities, as well as the A^ pressure 
pre-factors Pc, are advanced using the current state variables in Equations (2.35) 
and (2.28), following a simple Eulerian approach, for example v(t -|- 6t) = v(t) + 
{D/Dt)\(t)St. The time interval is chosen to satisfy the Courant-Friedrichs-Lewy 
condition (Courant et al., 1967). 

All the terms in Equation (2.35) are evaluated at grid points, while the terms in 
Equation (2.28) are evaluated within segments. The state variables and the unit tan- 
gent vector b are defined inside the segments. They are averaged between neighbors 
when needed at grid points. A similar approach is used when vectors defined at tube 
points are needed at segment centers. From the advanced positions, velocities, and 
pressure pre-factors Pc, it is possible to define all the advanced state variables. The 
density of each segment at a given time is calculated as p = {6m/ ^) {Be/ 61), where 
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61 is the updated length of the segment. Pressure follows from Equations (2.27) and 
(2.28), and the process is repeated until the end of the simulation. 

As the two tubes retract they move in opposite directions, but due to the sym- 
metry of the background configuration their dynamics are identical. Therefore, we 
present results only for the downward moving one. The simulated equations are 
dimensionless, but for concreteness the results are re-dimensionalized using charac- 
teristic magnetic field magnitude, density, and length. For this illustrative purpose, 
we choose Bq = 10 G, density po = nofn with electron number density uq = 10* 
cm"^, and m equal to half the mass of the proton. The initial angle between the 
pre-reconnection field lines is 2( = 90 degrees. The associated Alfven speed and tem- 
perature are vao — 2.2 Mm s"^, and Tff^ ^ 290 MK. We will consider a plasma at Tg 
= 1 MK, and the plasma-/? in this case is equal to (3 = 2Te/T^ ~ 0.007. For a typical 
length of Lq = 100 Mm, the Reynolds number is i?^ ~ 2982. An Alfven time (time to 
move one characteristic length at the Alfven speed) would be approximately 46 s. As 
mentioned before, the Prandtl number of the solar corona is typically equal to 0.01, 
and we will use this value in our simulations. The total number of segments in the 
tube is 4000. 

The simulation is carried out for approximately 3200 s, to allow for the full de- 
velopment of shocks (discussed in detail in the next section). The initial length of 
the tube was chosen to be unrealistically long in order to achieve the steady state for 
the shocks. As we will see in the next section, the thickness of the shocks, when fully 



50 

developed, is of the order of the length of an entire coronal loop. The fluid inside the 
tube evolves toward the steady state, which it probably would not attain within a 
realistic flare time. 

The left panel of Figure 10 shows the downward moving tube at different times. 
The tube is seen to be composed of three main segments, one horizontal segment 
moving downward at the Alfven speed related to the y-component of the magnetic 
field at the reconnection point {vAey = VAeSin{(), with VAe = vao), and two legs from 
the unperturbed parts of the initial tube. The two corners where these three pieces 
come together move at the Alfven speed along the legs of the initial tube and will 
be called bends (triangle shapes in the figure). In the right panel, the ?/-position of 
the right bend is plotted as a function of time for the same times as in the left panel. 
The solid line correspond to a line with a slope equal to —VAey 

The magnitude of the magnetic field remains constant along the tube due to 
constraint (2.33), therefore the area of the tube also remains constant (Equation [ 
2.21]). At the bends, the initially stationary plasma is deflected along the bisector 
of the angle between the two adjacent straight segments (average direction of the 
curvature force), as shown in Figure 11. In this figure, a time approximately equal 
to 2550 s was chosen to show the shape and internal structure of the tube long after 
reconnection. 
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Figure 10. Left panel: simulation results for the temporal evolution of a down- 
ward moving reconnected tube. Each horizontal segment corresponds to a different 
time. Vertical ticks account for the thickness of the two outward moving gas-dynamic 
shocks. The width of the shock evolves up to a steady value. Triangles indicate bend 
positions. Right panel: right bend y-position as function of time (same times chosen 
for the left panel). The solid line corresponds to a line with a slope equal to —VAey 
Both panels share the same y-axis scale. 

The velocity of the deflected plasma (region 2 in Figure 11) is (see Longcope et al., 
2009) 

V2 = -2vAe sin^ ^0 X - VAe sin(C)y- (2.42) 

This velocity is depicted in the top panel of Figure 11 (the left side is a mirror image). 

The parallel component (x direction) of this velocity is proportional to the Alfven 
speed, and is therefore generally supersonic [vAe = Cse\/'2/ll3). The tube gets shorter 
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Figure 11. Top panel: dashed line represents initial tube. Solid line represents the 
tube at time equal to 2550 s (vertical scale is compressed). Horizontal arrows delimit 
regions 1, 2 and 3 (post-shock plasma, inflow region, and unperturbed leg of the tube, 
respectively). Arrows show velocity directions in the tube (not to scale). Plasma is 
deflected in the bisector direction at the bends. Large and small vertical ticks indicate 
the beginning and end of the width of the shock, respectively. Bottom panel: the left 
side depicts temperature profile in half of the tube, and the right side, the plasma- 
P profile for the same time as in the top panel. Dashed lines show the theoretical 
steady-state Rankine-Hugoniot post-shock values. Crosses are the same as vertical 
ticks in the top panel. In both panels triangles indicate the position of the bends. 
There is no change in the state variables at the bends (Alfven waves). The tube has 
mirror symmetry in the x-direction. Both panels share the same x-axis scale. 



as it moves, creating parallel supersonic flows at the bends. These supersonic flows are 
directed toward the center of the tube where they collide. The small, but non-zero 
thermal pressure there prevents plasma from piling up and stops it in the parallel 



direction. The two colliding flows generate two gas-dynamic shocks (the magnetic 
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field remains constant across these shocks) that move outward from the center of the 
tube. In the top panel of Figure 11, velocity discontinuities in the parallel direction 
between region 1 (post-shock region) and region 2 (inflow region), are evident. The 
vertical ticks indicate the beginning and end positions of these shocks. In the bottom 
panel of the figure, the left side depicts the temperature profile along the tube, and 
the right side the plasma-/? profile, both for t = 2550 s. We see that temperature 
rises almost one order of magnitude from the pre-shock value, and the plasma-/? more 
than 30 times its initial value, and neither changes at the bends since these corners 
are Alfven waves. 

From Equation (2.39) and the boundary conditions at end-points of the tube, we 
see that the total energy of the tube is conserved (as shown in Figure 12). Even 
though, temperature changes significantly at the shock, only a small fraction (less 
than 10 %) of the available magnetic energy is converted to thermal energy (Figure 
12). The magnetic energy of the tube decreases with time as the tube shortens (the 
magnitude of the field remains constant in the tube, but its length decreases), and is 
mostly converted to kinetic energy at the bends. 

2.4. Gas-Dynamic Shock Internal Structure 

Transport coefficients (viscosity and thermal conduction) and their temperature 
dependence determine the inner structure of the gas-dynamic shocks. As the flux 
tubes retract, they follow a time- dependent evolution until they reach a steady-state 
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Figure 12. Energies' temporal evolution for the simulation presented in Section 2.3.2. 
Each line represents the temporal evolution of the corresponding specific energy inte- 
grated over the entire tube. The percentage of magnetic energy converted to thermal 
is small (less than 10 %). The rest is converted to kinetic energy. 



characterized by a transition between upstream and downstream states dictated by 
the Rankine-Hugoniot conditions (Rankine, 1870; Hugoniot, 1887). The diffusive 
terms smooth out the transitions but do not alter the upstream or downstream values 
provided these regions are asymptotically uniform {d'vjdl — )■ and dT/dl — )■ 0). 
These asymptotic values are governed by conservation laws observed by the diffusion. 

Gas-dynamic shocks in reconnected thin flux tubes were introduced by Longcope 
et al. (2009), in the case of ideal thin flux tube equations. It was shown there that 
only two kind of discontinuities are supported by reconnected thin tubes embedded in 
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a low plasma-/? uniform background: one that corresponds to a bend moving at the 
Alfven speed, where temperature, density and pressure are continuous, and a second 
one that corresponds to a gas-dynamic shock within a straight section of the tube. 
Both of these discontinuities were present in our simulation. In this section, we will 
focus on the later one, where heating and cooling occur (the heating at the bends is 
negligible) . 

To determine the effect of diffusive terms in the inner structure of the shocks, 
we follow an approach similar to the one used by several authors (Thomas, 1944; 
Grad, 1952; Gilbarg & Paolucci, 1953; Kennel, 1988). We consider a straight section 
{dh/dl = 0) and further assume the presence of a transition (shock) moving to the 
right at constant speed Vs- The inflow region will be denoted by subscript "2", and 
the post-shock region by subscript "1". In our reconnection scenario, the inflow region 
corresponds to the plasma moving supersonically after being accelerated by the bend 
(region 2 in Figure 11), and the post-shock region corresponds to the center of the 
tube (region 1 in the same figure). 

2.4.1. Steady-state Solution 

The Mach number (the ratio between the fluid velocity and the local sound speed) 
in the inflow region is determined by the initial plasma-/? and the initial angle between 
the field lines (Longcope et al., 2009), namely 




(2.43) 
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and the Mach number in the reference frame of the shock is 



M2,s = M2 + — = J^, (2.44) 

Cs,2 V ^-^2 

where 0^,2 is the sound speed in region 2, and U2 and T2 are the pre-shock velocity in 
the shock reference frame and temperature, respectively. For our simulation in the 
last section, the Mach numbers were M2 — 3.85 and M2^s — 5.3. 

In the reference frame of the shock, the inflow velocity corresponds to plasma 
moving at U2 = ^3.^2 ~ ^'s (< 0), and the post-shock flow is Ui = —Vg (plasma is 
stopped in the parallel direction). The ratio between the shock speed and the inflow 
speed magnitude, |fx,2|, is given by 



(2.45) 



\v^^2\ V Mi 16 4 
If constant states are assumed, d/dt = for all the quantities calculated in the moving 
frame of this discontinuity. In this reference frame, the steady-state, conservative, 
unitless TFT equations for uniform background case and low beta can be written as 

piUi = P2U2 = pu, (2.46) 
Piul + Pi = p2ul + P2 = pu^ + P-v-gj^, (2.47) 

Plul + UiPi^^ = P2UI + U2P2^^ = (2.48) 

7 — 1 7 — 1 

3 ^27 ^ du ^ dT 

= pu + uP - 2n?7— - 2k—, 

7 — 1 ol ol 
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The second equalities in each of the above equations refer to the state variables 
inside the shock. The usual Rankine-Hugoniot conditions can be found from the first 
equalities in each of the above equations. 

It is convenient for our analysis of the internal structure of the shock, to de- 
fine a characteristic velocity u* = ^ {u2 + T2/U2) = (W2/2) (l + 1/7M|^). With this 
definition, the Rankine-Hugoniot conditions can be written as 

P2 ^ Ml ^ (47-»77 +1)-U2 

Pi U2 U2 ' 

P . p,50) 

I2 -t2 

The plasma-/? ratio can be calculated as Pi/ P2 = {pi/ P2)(Ti/T2). In Figures 11, 15, 
and 14, horizontal dashed lines show the theoretical steady-state jump values for our 
simulation. Even though the post-shock plasma-/? is more than 30 times the initial 
value, it is still below unity. 

The second equalities in Equations (2.47) and (2.48) are differential equations 
for the shock's internal state variables. Their boundary conditions are given by the 
steady-state Rankine-Hugoniot conditions. Combining these two equations we obtain 

2(7 - 1) dT 2T - n2(7 - 1) + 4(7 - l)u*u + ^^(7 + 1) - 4jU2U* 



Pr u du — 2uu* + T 

h{T, u) 



(2.51) 



where we have defined functions h(T, u) and g{T^ u) as the numerator and denomina- 
tor of the above differential equation. indicates the Prandtl number. This equation 
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is valid for any temperature dependence of the transport coefficients, as long as both 
are the same; and can be numerically integrated for different Prandtl numbers and 
shock strengths. Numerical solutions of this equation in conjunction with Equations 
(2.47) or (2.48) provide a method to determine the thickness of the shock. 

In Figure 13, numerical solutions to the above equation are plotted for the same 
Mach number as in our simulation from last section (temperature and velocity are 
normalized to the pre-shock values). The pre- and post-shock values are fixed by 
the strength of the shock, independently of the transport coefficients, therefore we 
will call this kind of diagrams fixed point diagrams (they correspond to a velocity 
contrast-temperature plane described by Kennel (1988), or Navier-Stokes direction 
fields by Grad (1952)). Each solid line corresponds to a numerical solution to the 
above equation for a given Prandtl number (increasing thickness of the line represents 
increasing value in the Prandtl number). In this figure, the mixed dotted and dashed 
line shows the curve h{T, u) = that corresponds to the limit of infinite Prandtl 
number P^- As the Prandtl number increases from low values, the solutions (solid 
lines) approach this limit case. The dashed line represents the g{T,u) = curve, 
and asterisks indicate grid points (inside the shock region) from our simulation (P^ = 
0.01) for the same time chosen for Figure 11. Simulations coincide with the steady- 
state theoretical predictions. There are several grid points in each sub-region of the 
shock, which shows that the DEFT program can resolve the internal structure of the 
shock. 
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Figure 13. Fixed Point Diagram for super-critical Mach number equal to 5.3 (the same 
one than for simulation in Section 2.3.2). Temperature (ordinate) and velocity (ab- 
scissa) are scaled to the pre-shock values. For this case, > Top and bottom 
horizontal dotted lines represent post- and pre-shock temperatures, respectively. The 
leftmost vertical dotted line corresponds to the post-shock speed, and the rightmost 
dotted line the pre-shock speed. The dashed line represents the curve g{T, u) = 0, 
and the mixed dotted and dashed line depicts the h{T, u) = curve (the latter is the 
limit solution for an infinite Prandtl number). Each solid line represents a numerical 
solution for Equation (2.51) for different Prandtl numbers (increasing thickness of 
the line corresponds to increasing value of the Prandtl number). The shown Prandtl 
numbers are 0.01, 0.03, 0.1, 0.3, 1.0, and 3.0. The second vertical dotted line from the 
left shows the characteristic speed u* (maximum of the curve g(T,u) = 0), and the 
second dotted line from the right shows velocity Uq (separation between sub-shock 
and heat front for Pr = 0). Asterisks corresponds to grid points inside the shock for 
the simulation presented in Section 2.3.2 (P^ = 0.01) for time = 2550 s (the same 
time chosen for Figure 11). The internal structure of the shock is well resolved by 
the DEFT program (there are several grid points in each shock section). 



Velocity u* corresponds to the maximum value of the curve g(T,u) = in the 
T-u space, as shown in Figure 13. When the magnitude of this speed is smaller than 
the magnitude of the post-shock speed ui, the temperature increases monotonically 
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from T2 to Ti, without presenting an isothermal region, for all Prandtl numbers 
(Landau & Lifshitz, 1959, §Section 88). This occurs in weak shocks, called sub- 
critical, which have Mach numbers M2,s smaller than a critical Mach number M2,cr = 

v/(37-l)/(3-7)7- 

On the other hand, if the shock is strong, the fixed point diagram looks like the 
one depicted in Figure 13. In this case, if the Prandtl number is small (as in the 
solar corona and our simulation), the solution is composed of two regions (Landau 
& Lifshitz, 1959, §Section 88). In the first region, called the heat front or thermal 
conduction front, temperature is increased approximately along the g{T, u) = curve, 
from T2 up to approximately the steady-state value Ti. The second region, called 
isothermal sub-shock or isothermal discontinuity (or jump), consists of an almost 
isothermal {dT/du = 0) change of velocity from uq = Ti/ui to the steady-state value 
Ml. This last region corresponds to an isothermal compression of the plasma; most of 
the compression across the shock occurs in this region (see top panel of Figure 14). 
The heat front increases the temperature of the plasma without changing the density 
significantly. The compression ratio between the initial density p2 and the density 
at the transition between the sub-shock and the heat front po is bounded between 
the values (7 + l)/2 (limit for a infinitely strong shock) and (37 — l)/(7 + 1) (at the 
critical Mach number), for a Prandtl number equal to zero. For a gas constant equal 
to 7 = |, this corresponds to a density ratio between 1.33 and 1.5. 
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Figure 14. Top panel: electron number density evolution in the shock frame for times 
t = 63 s, 159 s, 319 s, 638 s, 1913 s, and 3156 s. Middle panel: temperature evolution 
in the shock frame for the same times as in the top panel. Bottom panel: volumetric 
heating rate in the shock frame for the same times as in the top panel. In these three 
panels, plus signs indicate the position of the heat front, crosses indicate the position 
of the isothermal sub-shock, and dark dots indicate the position of the transition 
between the isothermal sub-shock and the thermal front. The dotted lines indicate 
the earliest time. Horizontal dashed lines are the post-shock Rankine-Hugoniot cor- 
respondent values. 



2.4.2. Physical Size 

The thickness of the two internal shock regions, sub-shock and heat front, is 
determined by Equations ( 2.47) and ( 2.48), respectively. For the case of strong 
shocks at low Prandtl number, the thickness of the isothermal sub-shock is determined 
mostly by viscosity, and the heat front length by the thermal conduction, as shown 
below. To calculate these thicknesses analytically, we approximate our solution by 
the case where Pr is zero. In Figures 10 and 11, large and small vertical ticks indicate 
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the beginning and end of the thickness of the shock that correspond to the width of 
the sub-shock in addition to the width of the thermal front. 

In the bottom panel of Figure 15, the evolution of the plasma-/? of the tube is 
shown (in the reference frame of the shock). The thermal front extends to the right, 
and the isothermal sub-shock extends to the left. The top panel shows the evolution 
of the length of each region. 

The length (unitless) of the isothermal sub-shock, A/ss(i-!.o)! is calculated replac- 
ing T = Ti in Equation (2.47), then 



where integration is calculated from velocity ui to velocity uq = Ti/ui = 2u* — ui (this 
latter speed corresponds to the velocity at the transition between the sub-shock and 
the heat front as shown in Figure 13). The plasma-/? at this point can be calculated 
as /3o = P2U1U2/T2. In the bottom panel of Figure 15, the horizontal dotted line 
represents /3o- This value is used to approximate the position of the transition between 
the sub-shock and the heat front (zero horizontal coordinate in this figure). 

The integrand in Equation (2.52) diverges at both limits, therefore we will slightly 
modify them. We will define a smaller centered interval corresponding to a 95% of the 
original integration interval, and the integration will be performed with new limits. 
Defining a small interval Mint = {uq/ui — 1) /40, the analytical solution for the length 




(2.52) 
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of the sub-shock is prescribed as 



AZss 



Li°"(T2) log(u - 1) 




(2.53) 




UO 



where L'°°(T2) = r](T2) / p2Vth(T2)Lo is the unitless mean free path of the ions at 
temperature T2 (all the lengths are non-dimensionalized to characteristic length Lq); 
and Vth(T2) is the thermal speed of the ions at T2. For the simulation presented in 
Section 2.3.2, ^'""(Ts) ~ 0.9 x IQ-^. 

Expression (2.53) presents values for two shown in the last part of the 

equation. The first one corresponds to temperature independent rj, and the second 
one corresponds to 77 ~ T^. The length of the sub-shock is considerably different 
depending on which model is used for viscosity. For the simulation we presented in 
the last section (temperature dependence of viscosity given by Equation (2.30)), the 
ratio between post- and pre-shock temperatures was almost one order of magnitude, 
therefore the sub-shock length would be approximately 10^/^ times the one predicted 
by a constant viscosity model. 

For the parameters of our simulation, the bottom part of Equation (2.53) gives a 
length of the sub-shock approximately equal to 22 Mm, which is in good agreement 
with our simulation (see top panel of Figure 15). In this figure, the length of the 
sub-shock (negative positions) increases up to a steady-state value well approximated 
by the analytical solution for Prandtl number equal to zero. 
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Figure 15. Internal shock structure in the reference frame of the shock, for the 
simulation presented in Section 2.3.2. The origin is the transition between the sub- 
shock and the heat front, indicated with a solid circle. Both panels share the same 
X-axis. Bottom panel: plasma-/? as function of position in the tube for different times 
(64 s, 159 s, 319 s, 478 s, 956 s, 1913 s, 2550 s, 2869 s, and 3188 s). The dotted line 
represents the earliest time. Horizontal dotted line represents the value of /3o that 
corresponds to temperature Ti, velocity Uq, and density po- Horizontal mixed dotted 
and dashed line represent the theoretical post-shock values. Plus signs indicate the 
position of the heat front (they correspond to small vertical tick marks in top panel of 
Figure 11). Crosses indicate the position of the isothermal sub-shock (they correspond 
to large tick marks top panel of Figure 11). Top panel: temporal evolution of the 
heat front (solid line) and sub-shock (dashed line) lengths. Plus signs and crosses 
indicate the same positions as in the bottom panel. Leftmost vertical dotted line is 
the theoretical value for the length of the sub-shock calculated from Equation (2.53), 
and rightmost vertical dotted line is the theoretical value for the length of the heat 
front. 

When a constant viscosity model is used, the length of the sub-shock is smaller 
than the mean free path of the ions for temperature ratios larger than approximately 



65 

7 (see Figure 16). This is not true for the case of more reahstic temperature depen- 
dencies hke the one used in our simulation (rj ~ Ts). For this case, the length of 
the sub-shock is always more than one order of magnitude larger than that of the ion 
mean free path, supporting the use of fluid equations. The length of the sub-shock 
diverges when Uq = Ui (critical CcLSG^ , clS shown in Figure 16 (mixed dotted and dashed 
lines), and there is no sub-shock for Mach numbers below the critical one. 

To determine the length of the heat front, we will consider the portion of g{T^ u) = 
between velocities uq and U2, and replace the corresponding expression for the ve- 
locity u = u* (^1 + .^1 — j^^^ in Equation (2.48). If we define a unitless variable r = 



a/1 - T/{u*Y, and integral limits r2(M2,,) = yl - and ri(M2,,) = y 1 - 
the integral equation for the (unitless) thickness of the heat front becomes 



'7-l\ n K{T)TdT 

P2U2 V7 + 1/ in (r ± ri)(r - T2) 
where the sign in the denominator corresponds to super-critical cases, and the 
"— " sign corresponds to the sub-critical cases. In the sub-critical cases the above 
integral diverges at both integral limits, and for super-critical cases, only at T2. The 
critical case corresponds to ri = 0, and the integral is continuous at that point. To 
define a unique length for all the cases, we will integrate expression (2.54) in a centered 
interval corresponding to a 95% of the original integration interval by defining a small 
interval Tint = (t2 — ri)/40, analogous to the one use for the isothermal sub-shock. 
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Figure 16. Different lengths scaled with the pre-shock electron mean free path 
j^eiec^rp^^^ as function of the ratio of the post-shock to the pre-shock temperature. 
The second axis at the bottom shows the corresponding Mach number M2^s- Ver- 
tical dotted lines represent cases M2^s = 1 (no shock), M2,s = M2^cr (critical case), 
and M2,s = 5.3 (shock strength for the simulation presented in Section 2.3.2). The 
lengths in this figure correspond to the same Reynolds and Prandtl number used for 
the simulation presented in Section 2.3.2. Bottom horizontal dotted line represents 
the ratio of the ion to electron mean free path at T2, and top horizontal dotted line 
is the unity (electron mean free path at T2). Solid lines depict the length of the heat 
front A/hf (calculated from Equations (2.55) and (2.54)), for the cases labeled to the 
right. Mixed dotted and dashed lines represent the length of the isothermal sub-shock 
A/ss (calculated from Equation (2.53)), for the same cases, and dashed lines show 
the mean free path of electrons and ions at post-temperature Ti. 



The case of temperature independent thermal conductivity can be easily inte- 
grated since k{t) is a constant and can be taken outside the integral. The analytical 
solution for this case is 



A/ 



HF 



7-l^l /^L^i^'=(T2 
7 + 1 



27M2,. 



rilog(r ±ri) ± r2 log(r2 

n ± T2 



T2--rint 



n+Tint 



(2.55) 
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Here, L°''^'^(T2) = n{T2) / p2Vth{T2) Lq is the unit-less mean free path of the electrons at 
pre-shock temperature T2, and rrip and me represent the proton and electron masses, 
respectively. For the simulation presented in Section 2.3.2, L'^^^'^(T2) ^ 0.0021. 

The temperature-dependent thermal conductivity case has a slightly more compli- 
cated integral where k{t) = Kc\u*\^{1 — r^)^. For infinitely strong shocks {M2^s = co), 
expression (2.54) has a simpler form since ri(oo) = and r2(oo) = 1 (the diver- 
gence at T2 is removed). In this case the integral converges to a finite number, and 
we see that the length of the heat front is 00 since the factor \u*\^/p2U2 (that can 
be taken outside the integral) increases with the Mach number of the shock. The 
opposite occurs for the constant thermal conduction case. 

For finite shocks. Equation (2.54) can be integrated analytically or numerically. 
For the realistic coronal parameters used in our simulation, the length of the heat 
front is approximately equal to 93 Mm, in good agreement with our simulation (see 
top panel of Figure 15). This length is of the order of the length of a typical coronal 
loop, and one order of magnitude more than the length of the mean free path of the 
electrons at the post-shock temperature. Therefore, for a real coronal situation, the 
heat front will probably never reach its steady-state value since it will arrive at the 
end-points before completing its development. At one Alfven time (46 s), the length 
of the heat front is already approximately 28.0 Mm long (the evolution of the heat 
front is shown in the top panel of Figure 15). 
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2.4.3. Justifying a Fluid Treatment 

In Figure 16, the lengths of the isothermal sub-shock and heat front for temper- 
ature independent and temperature dependent models are plotted as function of the 
post-shock to pre-shock temperature ratio. All the lengths are scaled to the electron 
mean free path at temperature T2. Solid lines depict the length of the heat front, 
mixed dotted and dashed lines show the length of the isothermal sub-shock, and 
dashed lines display the particle's mean free path at Ti. 

Despite its complicated analytical form, the length of the heat front for the tem- 
perature dependent case follows a simple curve that for relatively strong shocks has a 
shape similar to the electron mean free path at Ti (dashed line). This gives some in- 
sight into the physics involved. The electron mean free path increases with post-shock 
temperature, as (Coulomb potential), while the thermal conduction also increases 
with post-shock temperature as Tt. The interplay between particle's mean free path 
and transport coefficient's temperature dependence determines the thickness of the 
shock. Whenever the diffusive coefficient's growth with temperature is faster than 
the particle's mean free path, the thickness of the shock is going to be larger than the 
particle's mean free path. Heat is transported larger distances ahead of the shock as 
the temperature ratio increases, maintaining a heat front length always larger than 
the mean free path. A similar analysis can be done for the sub-shock and the ion 
mean free path. When temperature independent transport coefficients are used, the 
results are several orders of magnitude smaller than the corresponding mean free 
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paths, which in the past created some concern about the use of fluid equations for 
shock studies. 

The vahdity of assuming fluid equations, as opposed to kinetic theory, can also 
be a concern if the ratio between the diffusion heat flux and the value for streaming 
particles is larger than a certain value (Campbell, 1984). We tested the coUisionality 
of the plasma in our simulation, and obtained ratios smaller than 0.004. This, and 
the above paragraph conclusion, support the use of continuum fluid equations. 

For the short times involved in a flare, the shocks may not achieve their steady 
state. The evolution of the isothermal sub-shock and the thermal front toward this 
state can be approximated by an exponential (see top panel of Figure 15). With 
this approximation, the characteristic time for the development of the sub-shock is 
approximately equal to 532 s, and its initial (maximum) speed is less than 10% that 
of the ion thermal speed. The characteristic time for the heat front is of the order 
of 188 s (four times the Alfven time), and its initial development speed is 3% of the 
electron thermal speed. We have tested the speed of the thermal front for a range 
of reconnection angles between ( = 20° to ( = 45°. For all these cases, this speed is 
between 3% and 5% that of the electron thermal speed. 

The electron number density and temperature evolution toward the steady state 
of our simulation are shown in Figure 14. The electron number density approaches its 
steady state from above the Rankine-Hugoniot value (top panel), and temperature 
approaches from below (middle panel). A short time after reconnection, the two 
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colliding fluids accelerated at the bends get overcompressed at the center of the tube 
(generating densities well above the expected value for inflnitely strong shocks), and 
they relax toward the equilibrium value. From the top panel of the figure, it is evident 
that most of the compression occurs at the sub-shock. The middle panel shows that 
most of the increase in temperature occurs at the thermal front. From the bottom 
panel of Figure 15, we see that pressure achieves its steady state in a characteristic 
time much smaller than the ones for density and temperature. 

The bottom panel of Figure 14 shows the volumetric heating rate. Most of the 
cooling occurs at the sub-shock where the plasma is isothermally compressed. Some 
cooling and heating occur at the thermal front. In this region, the plasma increases its 
temperature without considerably changing its density. In the sub-shock region, the 
velocity changes are large (see Figure 13), and it would be expected to have a strong 
positive contribution to the heating, however the cooling and heating are dominated 
by thermal conduction since the Prandtl number is small. The cooling comes from 
the negative curvature of the temperature to the power of 7/2 (not shown in the 
figure) as described in Equation (2.41). Even for a small temperature gradient in the 
sub-shock, the heating is negative. Without viscosity, however, it is not possible to 
develop the needed sub-shock to connect the pre-shock state to the post-shock state. 
On the other hand, it is possible to only have viscosity and no thermal conduction 
(limit of Prandtl number equal to infinity). 
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2.5. Conclusions 

We have presented new general fluid equations that describe the dynamics of thin 
flux tubes. These tubes are assumed to be isolated from their surroundings, except 
by total pressure balance with their exterior (their small local radius permits almost 
instantaneous equilibrium with the surroundings via fast-magnetosonic waves). This 
assumption neglects any drag force due to the external fluid. The ideal part of these 
equations is applicable to a wide range of conditions for external magnetic fields and 
pressure. We also proposed realistic field aligned diffusive terms for viscosity and 
thermal conduction, valid for strongly magnetized plasmas, that include their strong 
dependence on temperature (~ T^). 

We solve these equations numerically for the post-reconnection dynamics of a thin 
flux tube embedded in an uniform background with small plasma- /3 (/3 = 0.007). This 
tube was assumed to have been reconnected at an impulsive, small patchy region in a 
current sheet with skewed magnetic fields. Due to magnetic tension, the tube retracts 
at Alfvenic speeds reducing its length. Realistic coronal parameters (large Reynolds 
number and small Prandtl number) were used in the simulations. 

As the tube evolves, the initially stationary plasma is deflected at the tube bends 
(rotational discontinuities) and accelerated to super-sonic speeds toward the center 
of the tube where they collide. The small, but non-zero thermal pressure at the 
center of the tube stops the inflows, and the two colliding flows generate gas-dynamic 
shocks that move outward from the center of the tube. This shows the importance of 
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including the parallel thermal pressure gradient in the thin flux tube equations since 
there is no magnetic force in that direction. The magnitude of the magnetic field 
remains constant across the shocks (these shocks differ from the Petschek case since 
they are not switch-off shocks). The presence of diffusive terms in the thin flux tube 
equations provides a mechanism for the shock development. 

In the simulation, temperature rises almost one order of magnitude from the 
pre-shock value (1 MK) to the post-shock value (~ lOMK) across the gas-dynamics 
shocks. Nevertheless, only a small fraction (less than 10 %) of the available magnetic 
energy is converted to thermal energy. The magnetic energy of the tube decreases 
with time as the tube shortens, and is mostly converted to kinetic energy at the 
bends. This suggest that the rate of the reconnection mechanism has nothing to do 
with the final temperature achieved. This final value only depends on the strength of 
the shock determined by the initial angle between the field lines. 

The inner structure of the shock is related to the transport coefficients and 
strength of the shock (Mach number). For solar coronal conditions (strong shocks 
at low Prandtl numbers), the inner structure of the shocks present two sub-regions, 
a heat front where most of the temperature increase occurs (its width determined 
by thermal conduction), and an isothermal sub-shock where the fluid is compressed 
(its width depends on the viscosity). For these coronal conditions, the existence of 
a shock transition with thermal conduction alone is not thermodynamically permis- 
sible. We estimate the thickness of these regions assuming zero Prandtl number and 
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solving the governing differential equations analytically. The thicknesses of the sub- 
shock and heat front are many times the mean free path of the ions and electrons, 
respectively. This result, and our calculations regarding the free-streaming heat flux, 
lead us to conclude that the plasma is sufficiently collisional to use fluid equations. 

The higher density post-shock region should be readily observable. Post-shock 
density is almost four times greater than the initial background density, leading to 
sixteen times more emission measure (emission measure scales as the square of the 
density). The size of this region increases at a rate equal to the speed of the shock 
(Equation (2.45)), and at one Alfven time would be of the order of 10 Mm. There 
are several types of observations showing bright features of this scale (Masuda et al., 
1994; Tsuneta et al., 1997; Warren et al., 1999; Jiang et al., 2006). The radius of 
our tube is assumed small compared with this length, suggesting that any observable 
features will consist of multiple reconnection events in a small region, and several thin 
flux tubes being retracted outward from that region. 

The jump condition for thermal pressure is achieved quite rapidly. On the other 
hand, density and temperature achieve their steady-state values on characteristic 
times longer than that for thermal pressure. The latter approaches from lower values 
than the corresponding Rankine-Hugoniot post-shock temperature, and the former 
from higher values. For large Mach numbers (large initial reconnection angles), this 
temporary density overshoot may be larger than the superior limit for the Rankine- 
Hugoniot post-shock density ratio (— = ^^). The value of this ratio may exceed 
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one order of magnitude which would manifest as even stronger X-ray emissions near 
the reconnection site, that would decrease in intensity as the tube evolves. Most of 
the cooling occurs in the isothermal sub-shock region, and most of the heating in the 
thermal front. 

For real coronal situations, the steady-state thickness of the shock might never be 
achieved since the shock will arrive at the end-points before completing its develop- 
ment (the steady-state thickness value is of the order of the length of the loop). For a 
wide range of initial pre-reconnection angles (we tested angles ranging from C = 20° 
to C = 45°) the speed at which the heat front develops is between 3% and 5% that of 
the post-shock thermal electron speed. The heat front barely compresses the plasma, 
since most of the compression occurs in the isothermal region. The speed at which 
these fronts move along the legs of the tube depends on the reconnection angle. For a 
magnetic field of i?o = 10 G, and electron number density uq = 10^ cm"^, the range in 
speed for the thermal fronts is between 450 km s~^ and 650 km s~^ for angles between 
( = 20 degrees and C = 45 degrees. For larger angles, it is possible for the heat front 
to overtake the bends. 

Heat fronts are possibly relevant for chromospheric evaporation in flares. This 
is not directly treatable in our model since it is confined to the current sheet whose 
edge lays high above the chromosphere (see Figure 4). We expect the heat fronts 
found in our model to continue along the field lines toward their footpoints in the 
chromosphere. The temperature at the heat front increases approximately up to the 
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post-shock Rankine-Hugoniot value, which could be one order of magnitude larger 
than the initial temperature of the loop. At these high temperatures chromospheric 
material can be evaporated, and X-ray emissions from the foot points might be ob- 
served. The upward flow of evaporated chromospheric material would eventually 
encounter the sub-shock that is descending along the legs of the tube. What would 
happen after is a subject for further investigation. 

We neglected radiation in our model. The radiative cooling time at the pre- 
reconnection initial conditions is several orders of magnitude larger than the Alfven 
time. The density increase by as much as a factor of four would by itself decrease the 
radiative cooling time. Temperature, however, also increases one order of magnitude 
which dramatically increases the cooling time since the radiative loss function de- 
creases with increasing temperature at temperatures higher than 1 MK, (Cook et al., 
1989; Martens et al., 2000; Rosner et al., 1978). Therefore, radiation is a negligible 
effect in the post-reconnection process. 

Our model assumed a very simple initial geometry: uniform background condi- 
tions with pre-reconnection field lines that form an angle at both sides of the recon- 
nection. In future work, we will investigate a more realistic model with, for example, 
a Green-Syrovatskii (Green, 1965; Syrovatskii, 1971) current sheet with a guide field. 
This model presents new phenomena as the magnetic field magnitude changes along 
the tube. 
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3. DENSITY ENHANCEMENTS AND VOIDS FOLLOWING PATCHY 

RECONNECTION 

The layout of this chapter is as follows. The first section describes the background 
configuration where reconnection is assumed to happen. Three different locations 
of reconnection across the CS are described. The next section presents the TFT 
equations that govern the evolution of the reconnected flux tubes, as well as analytical 
equilibrium solutions to the perpendicular part of the equations. Section 3.3 describes 
simulations of the evolution of reconnected tubes for each reconnection location. The 
parallel dynamics of the tube is described in detail before the tube arrives at the 
top of the underlying arcade. In Section 3.3.2, the arrival at the arcade apex is 
simulated by an overdamped spring force exerted by the compressed arcade that 
slows down and stops the tube. After this arrival, the parallel dynamics changes and 
strong rarefactions develop. In Section 3.3.3, the temporal evolution of the tube's 
total emission measure and mean temperature are shown, and a method is described 
to determine where in the CS a tube was reconnected. The last section presents a 
discussion of our results and their possible observational consequences. 

3.1. Flare Current Sheet in the Solar Corona 

For the location of the reconnection episode in the solar corona, we assume a 
realistic magnetic field configuration like the one shown in Figure 17. This initial 
equilibrium background presents a CS in the x — y-plane, located high in the corona. 
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and skewed magnetic fields on each of its sides. Tlie matliematical description of tliis 
configuration is usually called Green-Syrovatskii (Green, 1965; Syrovatskii, 1971) CS, 
or "double Y-type" due to the shape of its separatrices when projecting field lines in 
the y — 2;-plane (see panel (b) of Figure 17 in the present chapter, and Figure 3(b) in 
Syrovatskii's paper). We included a uniform guide field in the ignorable x direction to 
provide an arbitrary angle between the pre-reconnection field lines. We also assume 
uniform initial values of density pe, and pressure Pg, on each side of the CS. The CS 
extends infinitely in the x-direction, but is finite in the ?/-direction. The bottom edge 
of the CS marks the top of the flare arcade (shown in Figure 17 as arc-shaped loops) 
with footpoints in the solar surface (bottom panel of the box). We will denote the 
half length of the CS in the ^/-direction as Lg. 

The entire region is assumed to be free of electrical resistivity, with exception of a 
short-lived small patch, somewhere in the CS where reconnection occurs (small sphere 
in Figure 17). The horizontal and vertical positions of the reconnection site will be 
denoted by xr and y^, respectively. Any horizontal position of the reconnection site is 
equivalent due to the initial 2. 5- dimensional symmetry of the magnetic configuration 
described above. For simplicity, we will assume it to be at the origin of the horizontal 
axis, then xr = 0. We choose the y-axis origin to be located at the center of the 
current sheet. 

We will distinguish three reconnection cases, based on their location on the CS. If 
reconnection happens in the top half of the CS (0 < < Lg), it will be called TOP 
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Figure 17. Flare current sheet geometry, (a) Field line configuration for a Green- 
Syrovatskii current sheet with a guide field in the positive x-direction. The transpar- 
ent gray rectangle represents the current sheet in the x — ?/-plane, located above the 
post-fiare arcade (gray arc-shape loops at the bottom). The current sheet is finite 
in the ^/-direction, but extends infinitely in the x-direction (2. 5- dimensional symme- 
try). The bottom plane of the box represents the solar surface. Black lines depict 
some magnetic field lines on the side of the current sheet that is closer to the viewer 
{z > 0), and the gray ones whose endpoints touch the top side of the box correspond 
to field lines on the back side of the current sheet [z < 0). The small sphere in the 
current sheet shows a generic patchy reconnection region. Field lines that intersect 
this region form a small bundle that reconnects. The dashed rectangle in the plane 
of the current sheet corresponds to the region shown in Figure 18. (b) Different view 
of the same field lines in panel (a), showing the double Y-type configuration. 

reconnection, if it happens in the middle of the CS {i/r = 0), it will be denoted as 
CENTER reconnection, and finally if reconnection occurs somewhere in the bottom 
half of the CS, it will be referred to as BOTTOM reconnection. Figure 18 represents 
the region of the CS enclosed by a dashed rectangle in Figure 17 and shows examples 
of the three different reconnection positions. For the sake of concreteness, we have 
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assumed that the half length of the CS is = 100 Mm, a typical length in the solar 
corona. 

The magnetic field near the CS can be described as the limit as 2 —> of the 
skewed Green-Syrovatskii magnetic field, as follows 

Be = SlBex±yBey]j^~- 

= x5e,. ± for z ^ 0. (3.1) 

Here, Bex and Bey are constants. Bex is the guide field, and Bey is the ?/-component 
of the usual Green-Syrovatskii field. The positive (negative) sign corresponds to 
the front (back) side of the CS. For simplicity, we defined unitless prime- variables, 
y' = y/Le and x' = x/Le. The lines y' = ±1 represent the edges of the CS. The 
magnitude of the external field is maximum at the center of the CS and decreases 
toward the edg GS, clS does the local Alfven speed. 

Only field lines in the vicinity of the plane of the CS that intersect the resis- 
tive patch will be reconnected, forming two opposite thin flux tubes. In Figure 18, 
spherical resistive patches are represented by their circular projections in the plane 
of the CS. Reconnected flux tubes have a finite thickness, as also does the CS. The 
tube's cross-sectional area at the reconnection region. A, is determined by the radius 
of the reconnection patch. This area and the magnitude of the magnetic field at the 
reconnection site determine the flux of the tube, = BA. 
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The actual size of the reconnection patches in the solar corona has not been 
determined, but if SADs represent the cross-sectional area of reconnected flux tubes, 
this value can be estimated to be only a few megameters in radius, as shown in Figure 
6. McKenzie & Savage (2009) reported dark voids sizes of approximately 10^ km^, 
with significant range for variation. The areas of the dark pockets can only give an 
estimation because they may change as tubes move, due to pressure balance with the 
background. Typical flare coronal loops extend for a hundred megameters, therefore 
it is reasonable to assume that the reconnected tubes are thin compared with their 
length. 

Field lines in front of the CS that intersect a given reconnection region will connect 
at this site with field lines on the other side of the CS. Figure 18 shows a sample 
of three pre-reconnection field lines on each side of the CS that intersect a TOP 
reconnection site located at = 3/4. We will assume that the pre-reconnection 
flux tubes are parameterized by one central field line (their axes, thicker lines in 
the graph), that crosses the center of the reconnection region. In this figure, only 
pre-reconnection representative field lines are drawn for the CENTER case. 

The half angle Ci? between the pre-reconnection representative field lines at a 
given reconnection location is a parameter in our model and can be calculated as 

tan (a) = I^V^l^ = tan (Co)^!^, (3.2) 



where Co represents the half angle that field lines make at the center of the CS. Both 
angles are shown in Figure 18. 

Any pre-reconnection field line near the CS can be described as 

y'{x') = ± sin [tan {Co)x' + c] , z ^0, (3.3) 

with c being a constant that can be determined by any point that intersects the curve. 
For field lines that cross the reconnection point (0,y^), c = ± arcsin(?/^), where the 
top sign corresponds to field lines in front of the CS {z > 0, thick solid lines in Figure 
18), and the bottom sign corresponds to field lines at the back of the CS {z < 0, thick 
dotted lines in Figure 18). 

The x-position at which a tube (its representative field line) is the closest to an 
edge of the CS can be found replacing y' = ±1 in Equation (3.3) (the positive sign 
corresponds to the top edge of the CS, and the negative sign corresponds to bottom 
edge of the CS), to obtain 

"I — arcsin (y'j^) 



/ 

•^bottom ~r 



tan (Co) 
I + arcsin (y^) 



z^O, 

z^O. (3.4) 



tan (Co) 

At these locations, the tubes are tangent to the edge of the CS since the only remaining 
component of the magnetic field is the guide field (Equation (3.1)). The tubes have 
legs that continue downward to the solar surface and these sections of the tubes do 
not generally lie in the same plane as the CS (see Figure 17). 
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After reconnection, the reconnected flux tubes are symmetric in the x-direction, 
but the initial 2.5-dimensional symmetry of the system is lost. The localization of 
the reconnection episode makes the problem three dimensional. Fortunately, the 
reconnected tubes are mostly two dimensional as they lie in the plane of the CS, 
and their retraction is along this plane (magnetic tension force is parallel to the 
plane of the CS). In Figure 18, representative lines for reconnected flux tubes in 
a BOTTOM reconnection site are shown. The dashed line represents the initial 
configuration for a reconnected flux tube whose left side was part of a tube on the 
back side of the CS and whose right side was part of a tube in the front side of the 
CS. This new tube is sharply bent near the reconnection site and will move upward 
as a transient feature that slides along the CS, between the flux layers. We call these 
tubes UPWARD moving tubes. The mixed dotted and dashed line in the figure is the 
initial configuration of the DOWNWARD moving reconnected tube that will retract 
in the negative ^/-direction. The directions of motion for each tube are indicated in 
the figure by arrows pointing outward from the BOTTOM reconnection site. The 
UPWARD and DOWNWARD moving tube definitions are valid for all reconnection 
cases. Any reconnected UPWARD moving tube is a mirror image of a DOWNWARD 
moving tube if both tubes have the same absolute value y'^ reconnection location, but 
with opposite sign. They will move in opposite directions, but their dynamics and 
shapes are equal. One difference between them is that the downward-moving ones are 
connected to a large mass reservoir (the sun's surface), whereas the upward-moving 
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ones are not. Since the tubes are assumed to be initially in equilibrium with their 
surroundings, this is only relevant if the tube footpoints are somehow perturbed and 
then interact with their boundaries. Our present analysis and simulations stop before 
this occurs. Therefore, in this paper we will only focus on DOWNWARD moving 
tubes. 

The shape of a reconnected tube changes with time mainly as a result of magnetic 
tension. This motion can be characterized by its representative field line arc length 
parameterization, R(/,t). At each point of the tube, the local tangent unit vector 
and curvature vector can be calculated from the parameterization as b = dH/dl and 
k = dh/dl, respectively. Figure 18 shows an example of the parameterization vector 
for an initial UPWARD moving tube, at a generic arc length /. Also, at this point of 
the curve, b and k are shown. Pieces of the tube like this one, that are far enough 
from the reconnection site, are initially in equilibrium since they were part of the 
original background configuration. There, the force due to the curvature of the tube 
is balanced by the force due to the external gradient of the magnetic field. 
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Figure 18. Excerpt of the plane of the current sheet, with three possible reconnection 
sites. This graph corresponds to the dashed box in Figure 17. Pre- reconnection field 
lines in front of the current sheet [z > 0) are represented by solid lines, and pre- 
reconnection field lines on the back side of the current sheet {z < 0) are shown as 
dotted lines. The edges of the current sheet are located at y = ±Le = ±100 Mm. 
The top circle shows a TOP reconnection site at i/r = |Le. The set of field lines that 
intersect this region (pre- reconnection tube) is represented by a central line (thicker). 
At this reconnection site, representative field lines on each side of the current sheet 
are skewed by 90°. The half reconnection angle {(r = 45°) is shown to the side of 
the TOP reconnection circle. The center circle represents a CENTER reconnection 
site with C/j ^ 57°. For this site, only the pre-reconnection representative field lines 
are drawn. The circle in the bottom half of the CS plane, located at yji = — |lve, 
represents a BOTTOM reconnection site with (ji = 45°. Here, the initial shape of the 
already reconnected tubes is drawn, and arrows indicate their direction of motion. 
The dotted and dashed line shows a reconnected tube that will move downward 
(DOWNWARD moving tube) due to magnetic tension at its center, and the dashed 
reconnected tube will move upward (UPWARD moving tube). For this last tube, a 
parameterization vector R(/, 0) is shown for a given arc-length /. Also at this arc 
length, the unit vector b is shown, as well as the curvature vector k. For all three 
reconnection cases, ~ 57°, therefore field lines on each side of the current sheet are 
parallel to each other for all reconnection positions. 
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3.2. Low Plasma-/3 Thin Flux Tube Equations 

Reconnection is assumed to happen in a small and short-lived region, therefore 
the reconnected tubes are thin and their evolution can be described by the TFT 
equations from the previous chapter (Guidoni & Longcope, 2010). These equations 
assume total pressure balance between the tubes and their surrounding plasma, and 
the tubes are assumed to be isolated from the background (no viscous momentum and 
heat exchange). The plasma in these tubes satisfies the frozen-in field condition since 
the resistivity that may have originated the reconnection episode was only non-zero 
at the reconnection region. 

In the solar corona, the plasma is strongly magnetized, ergo its plasma-/3 is small. 
For this reason, thermal pressure can be neglected with respect to magnetic pressure. 
This approximation and the assumption of total pressure balance constrain the mag- 
nitude of the magnetic field inside the tube; at all times, it is equal to the magnitude 
of the external magnetic field Be (Guidoni & Longcope, 2010) and can be calculated 
from Equation (3.1). 

The cross-sectional area of each element of the tube changes as the element moves 
between the flux layers of the CS, satisfying A= The density of each tube element 
of mass 5m varies according to p = 6mB(./4>6l, where 61 is the length of the element. 

The magnetic field direction is completely described by the positions of the tube 
elements since the frozen-in condition guarantees that fluid particles and field lines 
move together. These positions can be obtained integrating the velocity satisfying 
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the low plasma-/? TFT momentum equation (Guidoni & Longcope, 2010) 



^ (Bl\ , (Bl\ ^ d 



b?7 d\ 



(3.5) 



Dt dl / V4vry dl 

Here, p, v, P, and 77 are the mass density, velocity, plasma pressure, and viscosity of 
the fluid inside the tube, respectively. D / Dt is the advective derivative in the tube el- 
ement's reference frame and V_l represents the gradient in the direction perpendicular 
to the flux tube. 

The inviscid part of the above equation has a parallel-to-the-magnetic-field term 
with a derivative of the thermal pressure that cannot be neglected with respect to any 
magnetic pressure term since magnetic forces act only in the perpendicular direction 
(as shown in the second and third terms of the right-hand side). The viscosity term 
has components in both directions but is usually small since viscous Reynolds numbers 
in the corona are very large. The exception to this rule are shocks, where velocity 
gradients become significant. For strongly magnetized plasmas, the viscosity depends 
on temperature as r/ = r/cTi with rfc being a constant. 

The contribution of thermal conductivity to the change in entropy of the plasma 
is approximately two orders of magnitude larger than the viscosity. The change in 
entropy of the plasma (assumed to be an ideal gas) of any given tube element can be 
expressed as As = In (Pc/Pe), where Pc is related to the tube element's pressure and 
density as 

p{t) = p,(t)f^y. (3.6) 



Here, 7 is the adiabatic gas constant. 

The TFT entropy equation can be transformed into an equation for P^., as follows 
(Guidoni & Longcope, 2010) 



9v\^ „ d / K d 



(3.7) 



Dt " ' \ p 

The temperature dependence of the thermal conductivity is the same as for the vis- 
cosity, K = K.cT2, with Kc being a constant. The ideal gas constitutive relation is 
P = pksT /m = 2nkBT, with ks being the Boltzmann constant, fn the average 
particle mass, and n the electron number density. 

The TFT mass conservation has the following expression (Guidoni & Longcope, 
2010) 

— {^] = ^b-— . 3.8 

Dt\p ) p dl ^ ' 

One of our goals is to explain how reconnected flux tubes may become regions of 
density depletion, like in the dark voids mentioned in the introduction. To that end, 
the above equation can be re-arranged to show that the density of each tube element 
may change due to three different effects, 

^ln(rt = -^,n(^)-^ + v,.k, (3.9) 

where v\\ is the parallel component of the velocity and is the perpendicular one. 

The first term in the right-hand side corresponds to a change in area of the tube as 
the tube element moves to regions of different confining magnetic field magnitude. In 
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addition, if neighboring elements are slowing down ahead of the direction of motion, 
a pile-up occurs and the density of the element increases (its length gets reduced, and 
its area is fixed by the external field). The opposite effect occurs if the neighboring 
elements are slowing down upstream. The second term in the right-hand side reflects 
this effect, which is similar to what happens in traffic jams. The third term in the 
right-hand side is related to the expansion or contraction of a curved tube in its 
perpendicular direction. One example of this effect is a circular tube with fixed area 
expanding or contracting radially. In this case, the curvature vector always points 
toward the center of the tube, and the perpendicular velocity is parallel (contraction) 
or antiparallel (expansion) to this vector. Each tube element expands or contracts 
resulting in a change in density. We will show the relevance of some of these terms 
in Section 3.3.1. 

3.2.1. Joined Equilibrium Solution 

A simpler magnetic background configuration than the one assumed in the present 
chapter was studied by Guidoni & Longcope (2010). There, the background field is 
skewed and uniform, and the shape of the reconnected TFTs, as they evolve, can be 
described by a joined equilibrium solution (JES). We will call this uniform background 
model GuidonilO. In this case, the reconnected tubes, at any instant, are composed 
by three straight segments joined at two corners (the bends) that move at the Alfven 
speed along the initial tube. There, the magnetic field is rotated without changing its 
magnitude (RD). Each straight segment of the tube is an equilibrium solution of the 
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inviscid perpendicular part of the momentum TFT equation, but the joined solution 
is not in equilibrium due to the sharp angle at the bends. The two lateral segments 
are at rest and the central segment moves at the projection of the Alfven speed in the 
?/-direction. Since the magnetic field magnitude and the initial density are uniform, 
the Alfven speed is the same everywhere. 

The actual parallel dynamics in GuidonilO does not satisfy an equilibrium solution 
(constant pressure and temperature) since strong GDSs develop inside the tubes. 
Nevertheless, outside the shocks, the density and temperature are approximately 
constant. To completely satisfy an equilibrium solution of the general TFT equations, 
the change in velocity should satisfy b ■ ^ = 0; a condition that is not satisfied at 
the shocks. 

For our current non-uniform magnetic configuration (Equation (3.1)), it is also 
possible to find an analytical solution for the inviscid perpendicular part of the mo- 
mentum TFT Equation (3.5), where each mass element is in equilibrium (-^ = 0). 
In this case, the analytical solution differs quantitatively with the actual tube dy- 
namics, as we will see in the next section, but the general aspects of the solution are 
maintained. 

It is straightforward to see that a tube satisfying the following equation constitutes 
such an equilibrium solution 

k = |^ = ^V,ln(E,^). (3.10) 
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A possible general equilibrium solution of the TFT equations is given by a tube 
having uniform density and temperature that satisfies the above equation, as well as 
b ■ ^ = 0. In this section, we will only focus on the shape of the tube given by 
Equation (3.10), as the parallel dynamics is not expected to be in equilibrium. 

Since the reconnected tubes lie in the plane of the CS, the derivative with respect 



to the arc length of the tube can be expressed as d/dl = (l/Vl + x''^)d/dy', where 
x' = dx'/dy'. With this expression, after some algebra. Equation (3.10) becomes 

a/^^^^ = x'ii + x'^) = W V^J- ^ ^ 

The initial background configuration satisfies the general MHD equation's equilib- 
rium. Its field lines (Equation (3.3)) also satisfy the above equation. Therefore, 
pre-reconnection tubes are in equilibrium with respect to the general TFT equations. 
Immediately after reconnection, every tube element is also in equilibrium, except the 
central element that is sharply bent at the reconnection site. 

The reconnected flux tubes are symmetric in the x-direction and smooth at their 
center. A solution to Equation (3.11) that satisfies these conditions can be expressed 



as 



y{x) = y.cos 



tan (Co)^;' 



(3.12) 



yi + (l_^/2)tan2 (Co) 
This solution is a one-parameter family of curves whose parameter is the ^/-position 
of the center of the tube, y' . 
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It is possible to construct a JES where the shape of the reconnected tubes is 
composed of three different equihbrium regions: two unperturbed equihbrium side 
sections given by Equation (3.3) and a central curve described by the above equation. 
The left panel of Figure 19 shows such solutions for a CENTER reconnection case 
with half reconnection angle (r = Co = 45°and DOWNWARD moving tube. Each 
dotted line corresponds to a different parameter (position of the center of the tube, 
represented in the figure by triangles). The unperturbed side sections are part of the 
initial reconnected tube (solid lines in the figure). For example, a complete JES for 
the parameter y'^ labeled "C" in the figure corresponds to two side sections of the 
initial tube (sections "AB" and "DE") plus a central symmetric curve ("BCD"). For 
a given center position, the intersection between the dotted curve and the solid line 
(the bends, shown as circles in the figure) can be obtained setting Equations (3.12) 
and (3.3) equal to each other. 

The whole set of these JESs describes a reconnected tube that descends from 
the reconnection region. At each time, the tube is shorter than the original tube. 
Even though each section of the joined solution is an equilibrium solution, the joined 
solution is not. The bends are clearly out of equilibrium, therefore a complete curve 
like "ABCDE" does not satisfy equilibrium Equation (3.11). Left panels of Figures 20 
and 21 show similar joined equilibrium solutions for BOTTOM and TOP reconnection 
cases, respectively. 
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Figure 19. DOWNWARD moving tube for the CENTER reconnection case 
and reconnection angle Cfl = Co = 45° (note that y-axis and x-axis are 
not drawn to scale). Both panels share the same y-axis. Left panel: solid 
line represents the shape of the reconnected tube at t = 0. Each dot- 
ted line corresponds to a different JES with center position parameter y'^ = 
0.0, -0.08, -0.15, -0.23, -0.30, -0.38, -0.47, -0.55, -0.64, -0.79, and -0.98 (shown 
as triangles), respectively. Circles indicate the bends for the right side {z < 0) of the 
tube. Letters "ABCDE" join a complete JES curve for y'^ = —0.79. The horizontal 
dashed lines show extensions of the tube that simulate the legs that connect it to 
its footpoints at the photosphere. Right panel: DEFT simulation of the same initial 
reconnected tube as in left panel. Each central curved line represents the shape of 
the tube at a different time. The time interval between each curve is approximately 
5 s. Circles represent the theoretical positions of the bends (assumed to move at the 
local Alfven speed). Time in seconds is shown for some selected tube configurations 
to guide the eye. Large plus signs at the center of the tube indicate the beginning 
of the gas-dynamic shocks (heat front) and the smaller plus signs indicate the end 
of the shock (sub-shock). The solid black squares indicate the positions of a tube 
element as it moves along the tube. The arrow points to the location where the tube 
becomes tangent to the edge of the current sheet (^bottom) for the front side of the 
tube [z > 0). To the left of this point, the tube is extended as a straight line to 
simulate the leg that connects it to the photosphere (it corresponds to the horizontal 
dashed line to the left in the left panel). A similar straight line is added to the other 
end of the tube. 
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Figure 20. DOWNWARD moving tube for the BOTTOM reconnection 
case and reconnection angle = 45° and (o ~ 57°, with the same 

format as Figure 19. Left panel: the center position parameters are 
yj = -0.75,-0.78,-0.81,-0.84,-0.87,-0.90,-0.93,-0.97, and -1.01, respec- 
tively. Right panel: the time interval between each curve is approximately 2 s. 



For comparison purposes, all cases analyzed in the present section have a half 
reconnection angle = 45° (note that this is different for the CENTER reconnection 
case shown in Figure 18 that has ^ 57°; this choice visually simplified the graph 
making pre-reconnection field lines parallel to each other on each side of the CS, for all 
the reconnection cases). From now on, CENTER reconnection case implies = 45°. 



94 



Joined Equilibrium Solution Sinnulotion 
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Figure 21. DOWNWARD moving tube for the TOP reconnection case 
and reconnection angle = 45° and (o — 57°, with the same for- 

mat as Figure 19. Left panel: the center position parameters are y'^ = 
0.75,0.67,0.59,0.50,0.40,0.29,0.17,0.05,-0.08,-0.27,-0.43,, and -0.67, respec- 
tively. Right panel: the time interval between each curve is approximately 5 s. 

In the CENTER and BOTTOM reconnection cases, the concavity for the JESs 
is upward for all the central regions. On the other hand, for the TOP reconnection 
case, the concavity of the central regions changes from downward to upward when the 
center position changes from positive to negative (at y'^ = 0, the equilibrium central 
region corresponds to a straight horizontal line). 
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Equilibrium solutions (3.12) are valid as long as the position of the center of the 
tube is above the edge of the CS {y'^ > —1, or yc > 100 Mm). This corresponds to 
the arrival of the center of the tube to the top of the arcade. Since our analysis was 
restricted to field lines near the plane of the CS, equilibrium solutions beyond this 
point have no real significance. 

3.3. Simulations 

The low plasma-/? TFT Equations (3.5), (3.7), and (3.8) are non-linear, hence 
a time-dependent analytical solution of these equations is hardly a possibility. The 
JES, although far from being a general solution, presents the correct general aspects 
of the shape of the tube, as we will see in this section. The parallel dynamics are a 
different matter. The strong compressions of the plasma due to the shortening of the 
tube generate GDSs along the tubes that present discontinuities in the state variables 
(Guidoni & Longcope, 2010). 

To study the full dynamical solution of the TFT equations, we developed a com- 
puter program called Dynamical Evolution of Flux Tubes (DEFT) which solves the 
TFT Equations (3.5) and (3.7) in dimensionless form for the two reconnected tubes 
(DOWNWARD and UPWARD moving). Mass conservation is automatically satis- 
fied since the DEFT program implements a Lagrangian approach where each mass 
element of the tube is followed. The program uses a one- dimensional staggered mesh 
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where each tube piece is represented by grid points at its ends. A rather detailed 
description of this computer program was presented in the previous chapter. 

In this chapter, we will present three different simulations that correspond to the 
three cases described in Section 3.2.1. For all of them, the half reconnection angle is 
Cij = 45°. Results depend quantitatively on this angle, but we chose this intermediate 
value to show the general aspects of the solutions that are common for a wide range 
of angles. For these chosen cases, the half angle that field lines make at the center of 
the CS is Co ^ 57° for the BOTTOM and TOP cases, and Co ^ 45° for the CENTER 
case. 

The central initial part of the tubes are clearly out of equilibrium due to the 
sharp angle at the reconnection site. To avoid introducing length scales at the limit 
of resolution, the central part of the initial tube is smoothed. The end points of the 
tube are assumed to be fixed, and to have no temperature gradient (this last condition 
ensures no heat transfer from the end points). 

For concreteness, we will present the results from all simulations in units that 
assume a magnitude of the background magnetic field at the reconnection site of 10 
G, an initial uniform electron number density of rig = 10® cm"^, and initial uniform 
temperature = 1 MK. The Alfven speed VaeR at the reconnection site is then 
approximately 2200 km s"^. As in the above section, the half length of the CS is = 
100 Mm. We assume m equal to half the mass of the proton. Each tube has 1100 
points. 
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The unitless numbers for this simulation correspond to plasma-/? ^ 0.007, Prandtl 
number (ratio between the viscosity and thermal conductivity) Pr = 0.01, and viscous 
Reynolds number i?^ = LeVaeRPe/v — 2982. These values are typical values for the 
high solar corona. 

The simulated evolution of a DOWNWARD moving tube reconnected at the 
CENTER site is shown in the right panel of Figure 19. This simulation has similarities 
with the corresponding JES (left panel). For each given time, the tube consists 
of three sections: two unperturbed ones on its sides and a central one that moves 
downward. Times shown were chosen to coincide with the same bend positions as in 
the corresponding JES (circles are located in the same positions). The concavity of 
the central portion has the same sign as the JES, although the JESs are more curved. 

The flattening of the center part of the tube is due to the perpendicular compo- 
nent of the viscous force density in Equation (3.5). In the flgure, the small crosses 
near the center of the tube correspond to the location of sub-shocks, where parallel 
velocities have large gradients and viscosity becomes relevant. At those locations, 
the perpendicular viscous force acts in the direction opposite to the local curvature 
vector; therefore, the total force in the curvature direction is reduced. This increases 
the curvature of the tube, resulting in a change of curvature sign at the center of the 
tube. A thorough analysis of the shocks is presented in the next section. 

In the simulation, the center of the tube moves with approximately constant 
speed equal to the Alfven speed related to the y-component of the magnetic fleld at 
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the reconnection site, which is of the order of 1500 km s"^. This can be understood, 
at least for earher times, by comparison with the GuidonilO model where this is also 
the case. The initial tube can be approximated as the initial conditions for skewed 
uniform magnetic fields (same as in GuidonilO). The magnitude of the local Alfven 
speed decreases toward the bottom edge of the CS, therefore there are sections of 
the tube, like the center one, whose perpendicular speed is super-Alfvenic. The ratio 
between some tube elements' perpendicular speed and the magnitude of the local 
Alfven speed can reach values above 1.1, although this value is relative since the 
direction of the background magnetic field is different than the one of the tube. 

In over a minute, the central part of the tube arrives at the lower edge of the CS 
which marks the beginning of the top of the arcade, located a.t y = —100 Mm. The 
evolution of the tube after its arrival at the top of the arcade is discussed in detail in 
Section 3.3.2. 

In the GuidonilO model, the bends are RDs that move at the Alfven speed. We 
expect similar behavior in our present case. The positions of the circles in the right 
panel of Figure 19 were theoretically calculated as moving at the local Alfven speed 
along the tube. They coincide with the actual bend positions in the simulation. They 
continuously slow down in the y-direction as they approach the edge of the CS. Their 
speed in the x-direction, determined by the guide field, is constant. Being faster than 
the bends in the ^/-direction, the tube's center is the first to arrive at the top of the 
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arcade. This difference in speed is what determines the general upward concavity of 
the moving part the tube. 

At the bends, the initially stationary plasma is deflected along the bisector of 
the angle between the two adjacent tube portions (average direction of the curvature 
force). The resulting speed is Alfvenic and has a parallel component toward the 
center of the tube. These bends, as in GuidonilO, are RDs where only the magnetic 
field direction and the plasma velocity are changed. Notably, the angle at the bends 
remains approximately constant during the entire evolution of the tube, and equal to 
180 — (r degrees (same as in the uniform case). 

The BOTTOM reconnection case simulation (shown in left panel of Figure 20) 
presents the same characteristics as the CENTER one, but the time it takes for the 
tube to arrive at the top of the arcade is much shorter (~ 16 s), as the reconnection 
site is closer to the apex of the arcade. 

On the other hand, the TOP reconnection case (right panel. Figure 21) presents 
some differences. The general concavity sign of the tube changes with time, and this 
change occurs much later in the evolution of the tube compared to the corresponding 
JES (left panel). Initially, the bends accelerate in the vertical direction (the local 
Alfven speed increases toward the center of the CS) and move faster than the center 
of the tube, leading to a downward concavity (same as in JES). After the bends 
reach the center of the CS, they start decelerating as they move toward the top of 
the arcade. The tube's center accelerates continuously, catching up with the bends 
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after approximately 70 s, where the tube is almost flat. Meanwhile, the perpendicular 
viscosity force density increases the curvature near the locations of the shocks, curving 
the center part in the opposite direction with respect to the general curvature of the 
tube. For this TOP case, the center of the tube is not the first part to arrive at the 
top of the arcade. The ratio between some tube elements' perpendicular speed and 
the magnitude of the local Alfven speed reaches values above 1.7. 

3.3.1. Parallel Dynamics 

The parallel dynamics along the reconnected tubes presents complex features such 
as shocks and coupled nonlinear waves. There is no simple equilibrium counterpart 
to compare it with, as there was for the perpendicular motion. Nevertheless, its 
complexity is what makes the parallel dynamics very rich in behavior, and this section 
is devoted to its description. 

The bends accelerate plasma supersonically toward the center of the tube where 
they collide and are stopped abruptly. As a result, two strong GDSs are launched 
outward. Guidoni & Longcope (2010) describe in detail their inner structure for the 
uniform and skewed magnetic field configuration. It consists of one thermal front, 
where most of the heating occurs, and an isothermal sub-shock, where most of the 
compression occurs. The thickness of the GDSs increases with time and can be of 
order of tens of megameters, which is comparable to the size of the entire reconnected 
tube and larger than the particle's mean free path. The large and small plus signs 



101 

in the right panels of Figures 19-21 represent the beginning (heat front) and the end 
(sub-shock) of the GDSs, respectively, for each simulation time. 

For the three reconnection cases presented in this chapter, the GDSs are ahead of 
the bends at early times after reconnection. Initially, the shocks develop their inner 
structure quite rapidly until they achieve a near-equilibrium thickness. Meanwhile, 
the central part of the tube is increasing in length due to the difference in vertical 
speed between the center and the bends. Therefore, the central portion eventually 
becomes larger than the distance covered by the shocks, and the bends move ahead. 
For the TOP reconnection case, this situation is reversed after the bends move to the 
lower half of the CS, where they slow down and the shocks overrun the bends again. 

The electron density has a sharp increase at the sub-shocks where plasma is 
extremely compressed. Here, the second term of the right-hand side of Equation 
(3.9) is the most relevant to determine the electron density. The top panel of Figure 
22 shows density profiles for the CENTER reconnection case where each solid line 
corresponds to a different time in the simulation. At early times, the density at 
the center of the tube is over six times the background density. Even though the 
theoretical density jump limit for steady-state GDSs is 4, we achieve higher values 
in our simulations because there is an initial transient overshoot of density at the 
center of the tube. The plasma requires a finite time to adjust to the correct steady- 
state shock strength. Density decreases as times goes by, although never reaching the 
equilibrium value. In the figure, we include an inset to the left that reproduces the 
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t = 60 s electron density profile, where each point of the graph corresponds to a grid 
point of the tube. This shows that shocks are well resolved by the DEFT program. 

Density enhancements by shocks are a common feature in any standard recon- 
nection model; this compression occurs at slow-mode shocks or switch-off shocks, 
which are close relatives of the gas-dynamic shocks of our paper (Longcope & Brad- 
shaw, 2010). Therefore, dense outflows are usually expected from reconnection. In 
our patchy reconnection case, central long-lasting, high-density post-shock regions, 
or plugs, expand and move outwardly from the reconnection site. Their increase in 
size can be seen in the right panels of Figures 19-21 as tubes descend. The section of 
the tubes limited by the small crosses corresponds to the plug. This feature should 
be observable by an X-ray instrument of sensibility suitable to the emission measure. 
In fact, bright X-ray features apparently shrinking have been observed (McKenzie, 
2000; McKenzie & Hudson, 2001). Sheeley & Wang (2002) also describe some coro- 
nal inflows as barely visible density enhancements that move sunward, with dark tails 
forming behind them (sinking columns). 

Between the bends and the sub-shock, the electron density decreases before in- 
creasing again at the sub-shocks. This decrease in density is caused by the lateral 
expansion of the tube as it moves into regions of lower magnetic field (first term in 
the right-hand side of Equation (3.9)), as well as by the expansion of the curved tube 
in its perpendicular direction, although the latter effect is smaller than the former. In 
the figure, the positive x-positions of the bends are shown as circles for each chosen 
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Figure 22. DOWNWARD moving tube for the CENTER reconnection case. All 
panels share the same horizontal axis. Top panel: electron density along the tube 
for different times in the simulation. The x-positions of the bends (only for the right 
part of the tube) are shown as circles for each chosen time. Lower densities in the 
central part of the tube correspond to later times in the simulation. The inset figure 
repeats the electron density profile for the time t = 60 s, where each point of the 
curve corresponds to a grid point of the tube, to show that the sub-shocks are well 
resolved by the DEFT computer program. Middle panel: temperature along the tube 
for the same times as in the top panel. Large and small crosses indicate the beginning 
(heat front) and end (sub-shock) of the gas-dynamic shocks, respectively. The inset 
figure repeats the temperature profile for the time t = 60 s, where each point of the 
curve corresponds to a grid point of the tube, to show that thermal fronts are well 
resolved by the DEFT computer program. The arrow points to the location where 
the back side of the tube {z < 0) becomes tangent to the edge of the current sheet 
(a^bottom)- Bottom panel: parallel (along the tube) speed profile for the same times as 
in the top panel. Circles are also the same as in the top panel. Times in seconds are 
shown for each bend position. 
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time. After 60 s, the plasma density in these regions is less than 70% that of the 
background density. If observed, this tube would present dark and elongated thin 
regions (low emission measure) that get darker, descend toward the arcade, and grow 
in length as the bends move along the legs. This resembles the sinking columns or 
the SADs mentioned in the introduction. By the time the tube arrives at the top of 
the arcade, it would be mostly dark, and approximately half of the dark region would 
be cold, and the other half would have a range of higher temperatures as the heat 
fronts move outwardly. 

The GDSs are strong and the temperature jump across them may exceed an order 
of magnitude. The CENTER reconnection case's temperature profiles are shown in 
the middle panel of Figure 22 for the same times as in the top panel. The temperature 
across the shocks initially increases with time, achieving values higher than 7 MK, 
and later decreases. The large and small cross signs are the positions of the beginning 
and end of the thickness of the shocks, respectively. The length of the thermal fronts 
increases with time as they move outwardly from the center of the tube. The inset 
to the left reproduces the t = 60 s temperature profile, where each point of the graph 
corresponds to a simulation point in the tube. 

The thermal pressure gradient in Equation (3.5) is crucial for the generation 
of the shocks. If it were erroneously neglected in this equation, the plasma would 
continuously pile up at the center of the tube without limit. However, this term is 
not the only factor in the determination of parallel dynamics. The parallel, inviscid 
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part of Equation (3.5) can be re-written in the following way 

Dvn Dh IdP , , 

-Dt = "^'Dt-paT' ^^-''^ 

showing two possible sources for the change in tube element's parallel velocity (except 
perhaps inside the sub-shocks where the viscosity contribution may also be large). The 
first term on the right-hand side is a fictitious force caused by the change in direction 
of the unit parallel vector as tube elements move. The solid dark squares in Figures 
19-21 show the position of a generic tube element as it moves along the tube. 

The parallel velocity profile for the CENTER reconnection case is shown in the 
bottom panel of Figure 22, for the same times as in the other two panels. The positive 
x-positions of the bends are also shown as circles for each chosen time. There, the 
plasma at rest is abruptly accelerated to supersonic speeds, therefore the slope of 
the parallel velocity curve is almost vertical at these points. The same occurs at 
the bends on the other side of the tube where the inflows are directed to the right 
(positive speed). Since the angle at the bends remains mostly the same as time goes 
by, the value of the speed achieved at the bends is also approximately constant. The 
last time shown in this figure is 60 s, before any part of the tube has arrived to the 
top of the arcade. 

The parallel velocity after the bends slightly increases due to a small predomi- 
nance of the first term in Equation (3.13). This term also nonlinear ly couples the 
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parallel dynamics (sound waves) to the perpendicular one (Alfven waves), which in- 
troduces dispersion. The oscillations in parallel velocity that follow the bends are 
not an artifact of the simulation; they are real dispersive waves due to this coupling. 
This effect was not present in GuidonilO since the tube segments were straight at 
all times. These waves are smoothed out by diffusive processes, mostly related to 
thermal conduction. 

The following slow parallel deceleration coincides with the thermal front positions, 
where the second term of Equation (3.13) becomes relevant. There, the inflows are 
decelerated when tube elements encounter the heat front, until they are suddenly 
stopped at the sub-shocks (sharp decrease in parallel velocity near the center of the 
tube). In the post-shock region, the flows are reversed in direction with respect to 
the post-bend velocity. 

The BOTTOM reconnection case presents the same general features as the CEN- 
TER reconnection case for the three quantities plotted in Figure 22. These two cases 
are similar because their dynamics corresponds purely to the lower part of the CS 
where the background magnetic field decreases monotonously toward the edge of the 
CS. 

The TOP reconnection case presents some differences with respect to the other 
cases. The three panels of Figure 23 show the same type of graphs as in Figure 22, 
for the TOP reconnection case. The post-shock electron density (almost one order 
of magnitude higher than the ambient density) increases while the plug is in the top 



107 

part of the CS (the tube is being squeezed by the ambient pressure) and decreases 
after crossing the center of the CS. The squeezing increases the electron density after 
the bends while those regions are in the top half-plane of the CS, and a pronounced 
decrease in density occurs after the bends cross the center of the CS. Those regions 
below the center of the CS expand as the ambient pressure decreases. The decreased 
density after the bends can be half that of the background density by the time the 
tube arrives at the arcade. As the tube descends, a very dense and hot central region 
develops and moves downward toward the top of the arcade, while very dark and 
extended regions grow in length to its sides. 

The maximum post-shock temperature is high, over 16 MK. The thermal fronts 
are almost as long as the entire tube, as shown in the middle panel of the figure. 
There, the fictitious force is several times larger than the pressure gradient force and 
has opposite sign. After the bends, plasma is continuously being accelerated toward 
the center of the tube and stopped at the sub-shock. For a given tube element, it is 
possible to roughly estimate the value of the fictitious force, for a short period of time, 
as (|v^|) |A^|, where (|v_l|) represents the average magnitude of the perpendicular 
velocity in the considered period of time, and is the corresponding angle change 
in the parallel unit vector. We have estimated this term for this case, and it is several 
times larger than the corresponding pressure gradient term in Equation (3.13). 

The speed achieved at the bends is almost unchanged over the entire simulation 
as the angle at the bends does not change considerably through the simulation. We 
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Figure 23. DOWNWARD moving tube for the TOP reconnection case, with the same 
format as Figure 22. 

have run simulations with different reconnection angles and this angle seems to always 
remain approximately equal to 180 — C_r degrees. The post-shock parallel velocities are 
reversed with respect to the direction of the post-bend velocities, and the resulting 
speeds are high, almost equal in magnitude with respect to the speeds achieved at 
the bends. 

A zoom near the bends in the right panels of Figures 19-21 would show small 
oscillations in the shape of the tubes. There, the fictitious term is much larger than 
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the second term in the right-hand side of Equation (3.13), as the RDs do not change 
the pressure of the plasma. 

3.3.2. Top of the Arcade 

In the previous sections, the analysis of the evolution of the tubes was done for 
times when the tubes were above the edge of the CS. The DOWNWARD moving 
tubes will arrive at the top of the arcade after a finite time. For example, the middle 
point of the CENTER reconnection tube arrives at the top of the arcade shortly after 
a minute of travel. This arrival is usually referred in the literature as the place of the 
fast-mode shock or termination shock (Forbes & Priest, 1983; Forbes & Malherbe, 
1986; Forbes, 1986). Here, the reconnected flux tubes collide with the underlying 
arcade. 

It is likely that the parallel dynamics will continue evolving along the legs of the 
tube after the tubes had arrived to the top of the arcade. The legs extend downward 
to the solar surface and do not lie in the same plane of the CS, therefore the simple 
two-dimensional analysis we have done so far for the evolution of the tube cannot 
be applied. An exact analysis of the dynamics of the tube after the tube arrives 
at the edge of the CS escapes the scope of this paper. Nevertheless, we can make 
some approximations to study it. The decrease in length of the tubes after they 
arrive at the top of the arcade is probably much smaller than the decrease they 
undergo since reconnection, as the underlying arcade prevents them from shortening. 
By the time the tubes have arrived to the top of the arcade, their total lengths 
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decreased by 7%, 15%, and 28 % for the BOTTOM, CENTER, and TOP reconnection 
cases, respectively. They may also change their length, only by a small fraction, by 
decreasing the angle between the legs of the tube and its flat central part (panel (a) 
of Figure 24 shows a schematic where a reconnected tube (thick line) is lying on top 
of the underlying arcade). 

To simplify the analysis, we can assume that this length is fixed after the arrival at 
the arcade apex. The perpendicular motion of the tube is assumed to have stopped, 
and only the parallel dynamics along the legs would continue, mostly independently 
of how curved the tube is. To study this parallel dynamics, we will approximate the 
tube legs as straight line extensions on the sides of the tubes at the edge of the CS 
and parallel to the CS plane. Therefore, in our simulations an extra straight segment 
is added to the end of the tubes to simulate the legs; then the tubes are extended in 
the horizontal axis near the edge of the CS, as shown in the left panel of Figure 19 
(dashed lines). 

Most of the observed dark voids slow to a stop when they arrive at the top of the 
arcade (Sheeley et al., 2004). Linton et al. (2009) show through three-dimensional 
MHD simulations that reconnected flux tubes decelerate rapidly when they hit the 
Y-lines of a skewed Syrovatskii CS and the sheared arcade beyond them. Therefore, 
we chose in this first analysis to ignore the possibility of full "bounce" back of the 
tube toward the reconnection site. We propose a scenario where tubes are slowed 
down and stopped there by a damping force. We simulate this effect as a restoring 



Ill 



(a) 



A 



\ 



\ / 
\ 



era 



(b) 



(c) 



Before arriving at the top of the arcade 



Tube cross section 



/ 

/ / 
/ / ^ 

' ' ' 



^ \ \ 

- \ \ 

Arcade \ \ \ 



After arriving at the top of the arcade 



i 



Arcade restoring 
force 



\ 



/ 

' ' ' ^ ■ 



Figure 24. Tube arrival at the top of the arcade cartoon, (a) Reconnected tube (dark 
sohd hne) hes on top of the arcade after having descended from the reconnection site. 
Dashed hues represent the underlying arcade. The legs of the tube connect the parts 
of the tube that were close to the CS with the tube's footpoints at the solar surface, 
(b) Cross section of the tube (solid black circle) before its arrival at the top of the 
arcade (dashed parabolic lines), (c) Arcade restoring force. The underlying arcade 
gets compressed and lightly curved due to the arrival of the reconnected tube, and 
these stronger field lines exert an upward force on the tube. 
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force, as shown in panel (c) of Figure 24. The underlying arcade gets compressed 
and lightly curved due to the arrival of the reconnected tube, and these stronger field 
lines exert an upward force on the tube. We add a perpendicular damping and spring 
force densities to Equation (3.5), active only below the edge of the CS 

f = pfi(fi2y-z/fi- v) for < -1. (3.14) 

Here, n indicates a unit vector perpendicular to the tube, fi^ is the spring constant, 
and V is the damping coefficient. This force only acts in the perpendicular direction, 
slowing down the tube and stopping it. 

The above force density will result in overdamped oscillations when the ratio 
between the damping coefficient and the spring force, ^ = u/2Q, is larger than one. 
We opted for this in our simulations in Section 3.3. The cases ^ = 1 and ^ < 1 
correspond to critical and underdamped oscillations, respectively. 

Different arcade damping coefficients lead to different damping times, but do 
not considerably affect the parallel dynamics of the tube. We ran simulations with 
coefficients differing by orders of magnitude and all results were similar. For the 
results presented in this paper, the simulations were carried out with coefficients 
equal to ^ = 2.0 and ~ 0.49 s"^. 

In our simulations, the reconnected tubes retract from the reconnection site to 
the top of the arcade, and slow down due to the damped spring force from the arcade, 
until they lie flat on top of it. Then, the flnal shape of a DOWNWARD moving tube 
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would be a straight line on top of the arcade (see right panels of Figures 19-21 for 
late times in the simulation). 

Since the bends decelerate in the lower part of the CS, they are not the first 
sections of the tube to arrive at the top of the arcade. The first regions to do so 
become flat and grow in size as more tube elements descend atop the arcade (compare, 
for example, t = Q5 and 70 s of the CENTER simulation in Figure 19). We will call 
the outer edges of these fiat regions "second bends." 

Some of the parallel velocity change is due to the perpendicular dynamics captured 
in the first term on the right-hand side of Equation (3.13). After a section of the tube 
arrives at the top of the arcade, its perpendicular velocity is suddenly decreased 
by the damping spring force. Therefore, the first term of Equation (3.13) becomes 
unimportant compared to the other terms in the equation. The top panel of Figure 25 
shows parallel velocity profiles for late times in the CENTER reconnection simulation 
(only the positive x-axis side). The dashed line corresponds to simulation time equal 
to 60 s -the earliest time shown in this graph and the last time shown in Figure 22, 
just before any part of the tube arrives at the arcade. The sharp discontinuity on 
the left shows the deceleration at the sub-shock. The next times shown in the graph 
present a new discontinuity at the locations of the second bend, where the parallel 
velocity is sharply decreased by the sudden reduction in perpendicular speed. The 
"SB" arrows point toward the second bend locations at t = 65 s and t = 70 s. 
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When the second bend and the original bend become close (see arrow labeled "SB- 
B" in Figure 25), the angle at the latter one decreases until it finally becomes 180° 
when the bend arrives at the top of the arcade. Then, the infiow speed is drastically 
reduced because it depends strongly on this angle. The gas becomes rarefied and 
the electron density decreases rapidly behind the bends. Here, previously accelerated 
plasma continues moving toward the center of the tube leaving a density depletion 
behind. By the last time of our simulation, the electron density decrease is more than 
40% of the background density. The bottom panel of Figure 25 shows the electron 
density evolution for the same times shown in the top panel. The almost vertical 
slopes to the left indicate the sub-shock positions that continue moving outwardly, 
even after the bends arrive at the top of the arcade. For the TOP reconnection case, 
the maximum electron density decrease achieved by our simulations is more than 50% 
of the background. 

The sub-shocks will continue moving outwardly and interact with the vacated 
regions. Although our simulations stop before this occurs, we hypothesize that the 
electron density will continue decreasing for \x\ < |xbottom| and the central hot plug 
will mix with these rarefaction waves, generating secondary rarefaction waves that 
would move toward the center of the tube and ultimately will disassemble the hot 
plug, as described by Longcope et al. (2010). 
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Figure 25. DOWNWARD moving tube for the CENTER reconnection case. The 
two panels share the same horizontal axis that corresponds to the right side of the 
reconnected tube. Top panel: the parallel component of the speed for times between 
60 and 150 s, at 5 s intervals. The first time (dashed line) corresponds to t = 60 s, 
which is the last time shown in Figure 22. Time t = 150 s is shown with a thicker line. 
Circles show the x-position of the bends. The "SB" arrows point toward the second 
bend locations at t = 65 s and t = 70 s. The "SB-B" arrow indicates a time when the 
second bend and the original bend are close. Bottom panel: electron density profile 
for the same times as in the top panel. The arrow points to the location where the 
tube becomes tangent to the edge of the current sheet. 
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3.3.3. Emission Measure and Mean Tem- 
perature 

The differences in dynamics between each reconnection case suggest an observa- 
tional signature of where in the CS a flux tube had been reconnected. The tube's 
average temperature, weighted by the differential emission measure per ffux, DEM(T), 
can be calculated as 

J2 DEM{T) T AT 
" E DEM{T) AT ' ^^'^^^ 

T>Te 

where the DEM(T) is computed for small temperature bins of size AT along the tube, 



as 



E EMU) 



\T{l)~T\<-^^ 



DEMiT) = • (3.16) 



Only tube elements whose temperature is inside the range of the bin are included in 
the sum. (T) refers to the heated material of the tube; it excludes segments at or 
below ambient temperature. Figure 26 shows the time evolution of the tube's average 
temperature (dotted line) for the CENTER reconnection case. 

The emission measure per flux of a given tube element, EM(1), is computed in 
the following way 

^ (3.17) 

With the above definition, the tube's total emission measure per flux becomes 

EM 



^ T{l)>T, 



EM{1). (3.18) 
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Figure 26. DOWNWARD moving tube for the CENTER reconnection case. The 
time evolution of the mean temperature of the tube (left vertical axis) is shown with 
a dotted line. The total emission measure per flux of the tube (right vertical axis) 
is also shown for the same times with a dashed line. Arrow "C" indicates the time 
when the center of the tube arrives at the top of the arcade. Arrow "B" indicates the 
time when the bends arrive at the top of the arcade. 

The above quantity considers only the emission measure of the heated parts of the 
tube (T > Te), excluding cool sections like the ones following the bends and preceding 
the heat fronts, shown in the middle panel of Figure 22. Ambient plasma emission 
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measure is not included either. Figure 26 shows the time evolution of the tube's total 
emission measure per flux (dashed line) for the CENTER reconnection case. 

In the figure, the arrow labeled "C" indicates the time at which the center of 
the tube arrives at the top of the arcade, and the one labeled "B" represents the 
arrival of the bends at the top of the arcade. Initially, the mean temperature of the 
tube increases at a faster rate than the emission measure. They both achieve a local 
maximum before time "C". The following negative slope for both curves indicates 
that these quantities would rapidly decrease if the tube were left to retract without 
the arcade interrupting its motion. Remarkably, they increase after this point. The 
length of the tube does not change considerably while lying on top of the arcade, but 
the heat fronts continue moving along the tube, increasing its average temperature. 
Our simulations end before the heat fronts arrive at the footpoints, but we hypothesize 
that the average temperature of the tube will continue increasing until they do so. 

The central hot plug lasts longer than a conductive cooling time of a region 
with similar temperature and size (Longcope et al., 2010) as it is maintained by the 
inflows. These quantities are therefore long lasting. The TOP reconnection case 
presents some differences. The mean temperature and emission temperature grow 
together, as shown in Figure 27. Here, the arcade also prevents a decrease in both 
quantities. 

In order to provide a general method to determine where in the CS reconnection 
happened for a given flux tube, we normalized the mean temperature and emission 
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Figure 27. DOWNWARD moving tube for the TOP reconnection case, with the same 
format as Figure 26. 



measure by their maximum achieved value during the entire simulation. For example, 
the maximum achieved mean temperature during the simulation for the CENTER 
case is (T)max — 6.2 MK (see Figure 26). Then, for each time in the simulation, the 
normalized mean temperature is defined as 



(T) (normalized) 



(T) 



(T) 



(3.19) 



MAX 



The normalized emission measure is defined in an analogous way. 

Figure 28 shows the dependence of normalized emission measure with normalized 



mean temperature for each of the reconnection cases. In this figure, for each time in 
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the simulation, the value of the normalized emission measure is plotted as a function 
of the value of the mean temperature at that same time. BOTTOM and CENTER 
cases (dotted and solid lines, respectively) are very similar to each other. For these 
cases, the mean temperature achieves its maximum before the emission measure does. 
On the other hand, the TOP case presents a different dependence. Both quantities 
initially increase simultaneously and curl at the top end of the curve. 

3.4. Discussion 

We have presented a model of tube dynamics, accompanied by simulations, fol- 
lowing transient and localized magnetic reconnection in a realistic coronal background 
configuration (Syrovatskii CS). We have shown that retracting reconnected flux tubes 
may present elongated regions devoid of plasma, as well as long lasting, dense central 
hot regions. The latter are created by GDSs at the center of the tube, consisting 
in long thermal fronts followed by an isothermal sub-shock. In general, the jump in 
density across the shock exceeds the maximum value predicted by Rankine-Hugoniot 
(Rankine, 1870; Hugoniot, 1887) conditions. These jump conditions are calculated 
assuming the shock is in steady state, which is not the case here. For the TOP 
reconnection case, the jump in density can be almost an order of magnitude. This 
descending plasma plug, although very thin, would be extremely bright. It is long 
lived compared to a free expansion of similar temperature and size region since it is 
maintained by the inflows generated at the bends. 
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Figure 28. Normalized emission measure as function of the normalized mean temper- 
ature of the tube for DOWNWARD moving tubes. The solid line corresponds to the 
CENTER reconnection case, the dotted line to the BOTTOM reconnection case, and 
the dashed line to the TOP reconnection case. 



GDSs are also present in reconnected tubes sliding through a CS with uniform 
skewed fields (Guidoni & Longcope, 2010). However, in our present model, these 
hot plugs respond to the change in background magnetic pressure that compresses 
or expands them laterally, depending on which region of the CS the tube is sliding 
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through. For Syrovatskii CSs, the magnetic field decreases toward the edges of the 
CS (Y-points). There, its magnitude is minimum (it is zero if there is no guide 
field). Therefore, it is possible for the tubes to move through regions of a decreasing 
magnetic field as they descend toward the Sun. This seems counterintuitive because 
usually magnetic fields in the corona are assumed to increase at lower heights. 

We have presented only three reconnection locations with a reconnection half 
angle of 45° for illustrative effects. We used generic dimensional values, but results 
scale with Be, Pe, and Lg. For all the cases presented, the maximum plasma-/3 achieved 
never exceeded unity. 

RDs (bends) move at the local Alfven speed along the legs of the tubes. Elongated 
low density regions are generated by lateral expansions behind the bends, when tubes 
move through the lower half plane of the CS. There, background magnetic pressure 
decreases toward the edge of the CS. The amount of density depletion depends on 
how curved the tube is, which is directly related to the perpendicular gradient in the 
background magnetic field (reconnection angle). On the contrary, in uniform back- 
ground fields, as shown in GuidonilO, tubes present only straight sections, regardless 
of the reconnection angle. No density depletion is expected in this case. When there is 
no significant density depletion, and the GDSs are allowed to extend for a detectable 
length, an observed signature of these tubes would be a descending hot plug with 
high density. A signature in spectra would also be expected from the parallel flows 
that might be detectable as Doppler shifts. 
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The achieved decrease in density could be as much as 30%-50% of background 
(pre-flare) values. This level of depletion agrees with observations (McKenzie & 
Hudson, 1999; Sheeley & Wang, 2002). These percentages are lower bounds since we 
have not considered any increase in the surrounding plasma density by, for example, 
chromospheric evaporation from previous reconnection episodes. 

The appearance of the retracting tubes changes considerably with the direction 
of the line of sight and with the location of the reconnection episode. If reconnected 
tubes (their observable parts) are seen from the same view as in panel (a) of Figure 
17, they would have loop shapes, although their orientation would be perpendicular 
to the arcade since they move in a plane that is parallel to the arcade axis. On the 
other hand, if the reconnected tubes are seen from the view shown in panel (b) of 
Figure 17, they would be seen as "hairpins" because they are tangent to the plane of 
the CS. 

For the TOP reconnection case, the concavity of descending loops when they move 
in the top part of the CS plane is the one that would be expected from a cusp-shaped 
loop. In this region, these loops would be bright at their sides (they are being squeezed 
by ambient plasma) and even brighter at their center. If the line of sight is parallel 
to the x-direction, this loop would appear like a bright small region, descending at 
Alfvenic speeds. After the bends cross the center of the CS, dark regions will develop 
behind them, and this view will change to vertical dark regions preceded by a bright 
region, until the center catches up with the bends and the situation is reversed. The 
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dark and bright regions are hot since the heat fronts move at speeds similar to the 
bend speeds. Their sizes are comparable to the entire tube. 

For BOTTOM and CENTER reconnection cases, their concavity is always U- 
shaped. If seen from the x-direction, a bright region is followed by an elongated 
vertical dark one that might be hot or not, depending on the location of the heat 
front. 

Observed coronal inflows can be bright or dark, or both (Savage & McKenzie, 
2011); they may also have loop- or tadpole-shapes (McKenzie, 2000; McKenzie & 
Hudson, 2001; Sheeley & Wang, 2002). It is possible that many of them are the man- 
ifestation of the same three-dimensional phenomena observed from different lines of 
sight. The location of the tube with respect to the CS also matters. If an observation 
occurs, for example, while a tube is in the top half-plane of the CS, only bright regions 
would be observed. On the other hand, if a tube is moving through the lower portion 
of the CS, it would present dark regions, sometimes preceded by bright regions, or 
sometimes followed by bright regions, as described above. 

A dense, hot region may not always be observed due to its small size, or its high 
temperature. For reconnection angles larger than 45°, it is possible to achieve post- 
shock temperatures that are much higher than the usual instrument passbands. On 
the other hand, after the tube's arrival at the top of the arcade, this region continues 
growing in time, and may become visible. The pile up of several newly reconnected 
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tubes may also increase the size of this region to observable size and would manifest 
stationary bright region on top of the arcade. 

Remarkably, signatures of reconnection persist in a single tube, longer than the 
Alfvenic transit-time required for the tube to relax. As the tubes move downward, 
they encounter the top of the arcade that lies beneath the bottom edge of the CS. 
Here, the downward motion is halted by the arcade. We simulated this effect with 
a perpendicular damped spring force exerted by the compressed arcade. With this 
force, the tube comes to rest on top of the arcade as a straight line. The parallel 
dynamics continues along the legs. After they stop, the tubes cannot decrease their 
length any longer, and the RDs are shut down. Then, the gas gets rarefied even more 
because the already accelerated plasma continues moving toward the center of the 
tube, increasing the density depletion behind it. 

Although our simulations stop before this, we hypothesize that the hot post-shock 
regions will continue moving along the tube and interact with the rarefaction waves, 
subsequently disassembling the hot plug. The legs of the tube would also brighten 
up as the thermal fronts descend toward the footpoints. 

The tube's parallel and perpendicular dynamics are nonlinearly coupled to each 
other, which can be seen, for example, in the oscillations of the parallel velocity pro- 
files. In addition, when the tube arrives at the arcade, the perpendicular dynamics 
is halted, and changes in the parallel speed are evident at the "second bend". The 
parallel velocity profiles are determined by pressure and fictitious forces due to the 
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motion of tube elements along curved paths. This can have important consequences 
for Doppler-shift observations. For example, the post-shock region in the TOP recon- 
nection case has very strong flows in a direction almost parallel to the surface of the 
Sun. This would manifest as coincident blue- and red-shifts along the line of sight, 
if seen from the x-direction, and perpendicular to the line of sight if seen from the 
z-direction. Perpendicular velocities can be higher or lower than the local background 
Alfven speed, which can be related to the generation of secondary shocks as the tubes 
move (the surrounding plasma may be shocked by the passing tubes). 

We described the temporal behavior of the total emission measure and mean 
temperature of the heated parts of the tube and provide an observational method 
that may indicate where in the CS the tube has been reconnected. Tubes that were 
reconnected in the TOP half of the CS present total emission measure that grows 
simultaneously with the mean temperature of the tube. On the other hand, tubes 
that have been reconnected in the BOTTOM part of the CS achieve a maximum 
in their mean temperature much earlier than the emission measure does. For all the 
reconnection cases, the emission measure and mean temperature duration is extended 
by the arrival of the tubes at the top of the arcade. After the arrival, the length of 
the tube remains the same but the heat fronts continue heating the plasma as they 
move along the legs. 
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If observed downflows (SADs and SADLs) are reconnection outflows, tlie question 
of why their speeds are lower than the assumed Alfven speed still remains. We re- 
stricted ourselves to standard reconnection scenarios where the outflows are Alfvenic. 
Nevertheless, if dark voids are related to the dark regions generated behind the bends, 
their speeds near the edge of the CS are decreased considerably from the one they 
had at the center of the CS. For instance, in the CENTER reconnection case, the 
bend speeds are reduced by a fourth of what it was at the reconnection site. This 
decrease in speed of the bends is not enough to explain the observed descending voids 
that move at half of the presumed reconnection Alfven speed. However, if the half 
reconnection angle is greater than 60°, the ratio of the Alfven speed at the center of 
the CS to the edge of the CS is larger than 1/2. 

HXR sources have been observed near top of SXT arcades, at the same time as 
the SADs descend to the top of the arcade (Krucker et al., 2003; Asai et al., 2004). 
This may be the result of the interaction between downflows and the fast shock, or 
could be related to reconnection at the arcade apex between the retracting tubes and 
the arcade. Apex reconnection is a possibility in our three-dimensional model, but 
one we have not yet explored. Field lines that arrive at the top of the arcade have a 
different angle than the underlying arcade (panel (a) of Figure 24) and arrive there 
at Alfvenic or super- Alfvenic speeds. This is not an option in purely two-dimensional 
models where the new reconnected field lines arrive at the arcade with magnetic field 
direction parallel to the arcade. 
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Maintaining a fast shock on top of the arcade may not be feasible at sub-Alfvenic 
speeds. The standing fast shocks exist as long as the reconnection jets exists and is 
super-magnetosonic, therefore it is a feature of a flare gradual phase (Forbes, 1986). 
Most of the SADs have been reported during the decay phase (there are some excep- 
tions to this: see Asai et al. (2004)), but their speeds are not Alfvenic. If SADs are 
the only reconnection outflows, it is not clear if the termination shocks would still be 
located at the top of the SXT loops. 

Another possibility could be that the termination shock is located further up from 
the top of the SXT loops. If loops hotter than the SXT loops (that cannot be seen 
in pass bands at lower temperatures) are located above them, the termination shock 
could be located higher in the corona. In this case, the downflows could be localized 
downstream flows from the fast shock, which would explain their low speeds. The 
models for fast shocks are two dimensional and it is not clear how this would apply 
to localized and time-dependent flux tubes like the ones described here. The down 
side to this explanation is that a fast shock certainly would compress the plasma 
downstream instead of vacating it. 

We have not considered other interactions between the reconnected tubes and 
their ambient plasma, besides pressure balance between them. The addition of a 
drag force to our equations could be an important improvement to our model that 
would contribute to the decrease in speed, as well as to make the model more realistic. 
We expect to include this effect in future work. The assumption of initial uniform 



129 

density could also be improved by assuming a stratified atmosphere. Here, tubes 
that have been reconnected in the corona will descend in regions of higher density, 
and seen darker than the background, as suggested by Savage et al. (2010), although 
this effect has not been corroborated by theoretical analysis. Another effect that 
we have not considered is the possible tangling of field lines (Klimchuk, 1997). For 
example, an UPWARD moving tube reconnected at a given location might run into 
a DOWNWARD moving tube retracting from a higher reconnection site. 

We have neglected any non-fiuid effects and therefore neglected non-thermal par- 
ticles and their acceleration. 

We do not claim that the depleted regions from our model can completely explain 
the observed dark voids, but it is suggestive that several observed phenomena have 
some similarities with our model. Hudson & McKenzie (2001) stated that "Dark 
outflows seem inconsistent with the idea of heating by reconnection" . We believe our 
work demonstrates the contrary, and thus makes the role of reconnection in flares still 
more plausible. 
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4. TESTING THE THIN FLUX TUBE MODEL 

Our model for reconnected thin flux tubes appears to explain some observed 
phenomena in the Solar Corona, as described in the previous chapter. However, it 
is still not clear how realistic or restrictive its assumptions are. The TFT model 
assumes that, if the reconnection region is small and short-lived, a reconnected tube 
will retract from this region and bends will move along its legs at the local Alfven 
speed. The model also predicts that shocks will travel outward from the center of 
the tube, heating and compressing its plasma. The tube is assumed to be isolated 
from its surroundings and to slide freely between the current sheet layers. The only 
interaction with its surroundings is through pressure balance. 

One way to test the validity of the foregoing assumptions is to compare results 
of the model with those from fully three-dimensional solutions of the general MHD 
equations. In this chapter, we present three-dimensional MHD simulations carried out 
in a volume filled with plasma with a Green-Syrovatskii (Green, 1965; Syrovatskii, 
1971) CS, like the one depicted in Figure 17. The reconnection episode is triggered 
by imposing a short-lived and localized magnetic resistivity patch on the top half of 
the current sheet. The results of these simulations are compared to the corresponding 
DEFT simulation. 

For such comparison, the possible options from the parameter space are limitless: 
reconnection location, field line angle at each side of the current sheet, plasma unitless 
numbers, etc. However, given the complexity of the MHD simulations, we only chose 
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one particular case for our comparison. Typical simulations keep track of millions of 
cells that describe a plasma in a given volume. They are computationally demanding, 
and can only be executed in reasonable amounts of time on super-computers. In the 
near future, more of such comparisons will be performed. 



The fully three-dimensional MHD simulations were carried out by Dr. Mark 
Linton at the Naval Research Laboratory. The computer program used was the 
Adaptively Refined Magnetohydrodynamics Solver (ARMS) developed by Dr. C. 
Richard Devore. The ARMS computer program is a time-explicit, conservative, mono- 
tone MHD model. The MHD equations are advanced using flux-corrected transport 
techniques (Devore, 1991), and the adaptive mesh refinement is managed by the 
PARAMESH toolkit (MacNeice et al., 2000). The grid adapts to the spatial varia- 
tions of the magnetic field magnitude and electric current density, therefore the region 
of the current sheet is refined at the highest resolution, and it is effectively one cell 
wide. Newly created grid points are linearly interpolated with a monotone limiter. 

The three-dimensional Syrovatskii current sheet is described by the following 
equation (Priest, 2000) 



4.1. 3D MHD Simulations 




(4.1) 
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where i is the positive value of V— 1, and u = y + iz. In this configuration, the 
current sheet plane lies on the z = plane. To allow for a finite angle between field 
lines on each side of the current sheet, a uniform guide field Bg^ is added. 

Extrapolative, zero gradient, open boundary conditions are imposed in the y- 
and z-directions, while periodic boundary conditions are imposed in the guide field 
direction. The maximum resolution in the x, y, and z directions is 256 x [3, 1, 1] cells, 
respectively. 

The program solves ideal MHD Equations (2.1), (2.2), and (2.3), with the addition 
of an explicit magnetic resistivity (A) term in the right hand side of induction equation 
(2.4) 

-V X (AV X B) , (4.2) 

and the corresponding ohmic heating term added to the energy equation (2.3) 

j^lVxBr (4.3) 

The simulations are in non-dimensional form. The density and gas pressure are 
initially set to be uniform at = 1/2 and Pi = 20/3. The chosen magnitude of the 
magnetic field at {x,y,z) = (0,0, ~ 0) is Bi = 20. With these initial conditions, the 
Alfven speed and the plasma-/? at the center of the current sheet are Vai = 20/\/2ti, 
and /3j = 0.42. The resistivity patch (reconnection region) is centered at = 2/3, 
close to the upper Y-line, so the upfiowing reconnected lines quickly hit that line. 
The chosen angle at the center of the CS is (q ~ 45°. This implies a reconnection 
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angle equal to (fi ~ 37°, and a plasma-/? there equal to 0.54. The unitless half-length 
of the current sheet is Lj = 1. 

Simulations are run with uniform background resistivity, Ae, locally enhanced at 
the reconnection location within a sphere or radius 6 = 0.175. This enhancement is 
short lived, for the first t* = 0.1 of the simulation, and has the form 

A = Ae(l + Ce-'-'/^'), (4.4) 

for r = \J x"^ + {y — yuf + z'^ < 26, and I + C = 10^. The background Lundquist 
number, Sx = 5vai/Xe, is approximately equal to 1.4 x 10^, and the Lundquist number 
for the peak resistivity at the center of the reconnection region is approximately 14. 

Including realistic transport coefficients, such as viscosity and thermal conduc- 
tion, in a fully three-dimensional simulation is extremely difficult. The viscous stress 
tensor components in the direction parallel to the magnetic field are many orders of 
magnitude larger than the ones in the perpendicular direction. The same is true for 
the heat conduction. The ARMS computer program does not include these coeffi- 
cients explicitly, but there is numerical resistivity due to the finite grid-scale of the 
simulation. This is analogous to having transport coefficients that are temperature- 
independent with a Prandtl number equal to 1. These limitations are quite drastic 
since this transport coefficient temperature dependence is non trivial in establishing 
the thickness of the shocks, as shown in chapter 2. 
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A set of Lagrangian test particles can be followed dynamically during the simula- 
tion, and field lines are traced from these particles (Linton et al., 2009). In particular, 
field lines that originally belonged to one side of the current sheet, are anchored at the 
bottom of the arcade, and intersect the sphere of high resistivity at the start of the 
simulation, can be traced in time. The particles that initially lie in the high resistiv- 
ity region are not frozen onto field lines, but they follow the field line evolution after 
they leave the reconnection region, or after the resistive sphere is turned off. Figure 
29 shows the evolution of one such field line, where each solid color corresponds to 
a different time. The 2;-coordinate has been enhanced for clarity, so the figure is not 
to scale. For concreteness, we will present our results assuming a current sheet half- 
length of Le = 100 Mm, initial uniform electron density ng = 10^ cm"'^, temperature 
Te = 1 MK, and magnetic field magnitude at the reconnection site, B^. ~ 1.1 G. In 
this units, the reconnection sphere radius is 5 = 17.5 Mm. 

The dashed lines on the coordinate box walls are the projections of the field lines 
on those planes for the corresponding colors. The thick black field line is the field 
line at t = 0. At this time, the field line is completely on the front side {z > 0) of 
the current sheet plane. This is easy to see on its projection on the bottom box wall. 
This field line intersects the reconnection sphere with center at {xr, yn, zji) = (0, 66, 0) 
Mm, and its bottom footpoint is anchored at {x,y, z) = (—251, —125,45) Mm. From 
this point, the field line is traced for different times. We will call this field line, LineO. 
Its top endpoint is initially located at [x, y, z) = (123, 125, 45) Mm. 
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Figure 29. LineO evolution for early times in the three-dimensional simulation. The 
colored solid lines are the field line curves for times at 30 s intervals: s (black), 
30 s (cyan), 60 s (magenta), 90 s (yellow), 120 s (blue), 150 s (green). The plane 
of the current sheet is located at 2; = 0. LineO bottom footpoint at t = s is 
point {x,y,z) = (—251, —125, —45) Mm. Its top endpoint at t = is located at 
{x, y, z) = (123, 125, —45) Mm. Colored dashed lines are the field line projections on 
the box walls of the same color three-dimensional curve. 
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While the resistive sphere is active, the top footpoint moves as the field line 
"diffuses" . This artificial motion of the footpoint continues, and later on the field line 
is connected to the other side of the current sheet (the blue curve has one footpoint 
in the z > region). Figures 30 and Figures 31 show the same process for field lines 
with bottom footpoints at {x,y,z) = (—251,-125,44) Mm (we will call this line, 
Linel) and {x,y,z) = (-251,-125,-43) Mm (Line2), respectively. These field lines 
also cross the reconnection sphere. While Figures 29, 30, and 31 have an enhanced 
2;-axis, Figure 32 shows field fines LineO, Linel, and Line2 to scale at t = s (before 
reconnection) and at t = 120 s (after reconnection); they change their connectivity 
as they reconnect across the current sheet. The red circles are the projections on the 
box walls of a sphere of radius 6 = 17.5 Mm, representing approximately the region 
of high resistivity. After reconnection, the field lines are symmetric around the x = 
plane. The sections of these field lines whose ^/-coordinates lie between y = —100 
Mm and y = +100 Mm (the edges of the current sheet) are almost fiat and parallel 
to the plane of the current sheet. The sections outside these edges are the tube legs 
and are not parallel to the current sheet. 

Once a field line is completely reconnected, its newest footpoint remains almost 
stationary, and the field line retracts from the reconnection region, decreasing its 
length. The evolution of LineO after the resistive sphere has been turn off is shown 
in Figure 33. Most of the retraction occurs in the plane of the current sheet, as can 
be seen from the field line projection on the bottom box wall. It was a principal 
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Figure 30. Linel evolution for the same times as in Figure 29. Linel bottom footpoint 
at t = is point (x, y, z) = (—251, —125, 44) Mm. Its top endpoint at t = is located 
at {x,y,z) = (132,125,44) Mm. Colors represent the same times as in Figure 29. 
Dash lines are also the field line projections on the walls. 
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assumption of the TFT model underlying previous chapters that post-reconnection 
motion was confined to the current sheet. 




Figure 31. Line2 evolution for the same times as in Figure 29. Line2 bottom footpoint 
at t = is point (x, y, z) = (—251, —125, 43) Mm. Its top endpoint at t = is located 
at {x,y,z) = (144,125,43) Mm. Colors represent the same times as in Figure 29. 
Dash lines are also the field line projections on the walls. 



If LineO is observed from the x-direction as it moves, it would look like a retracting 



hairpin (this can be seen in its projections on the background walls in Figures 29 and 
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32). On the other hand, if the field fine is observed from the z-direction, it would look 
like a retracting loop (as can be seen from the projections on the other background 
wall). This is very similar to what we described in our conclusions in the previous 
chapter. 




Figure 32. LineO, Linel, Line2 to scale. Solid lines represent field lines at t = 0, as 
well as their projections on the box walls, and dashed lines are the same color field 
line at t = 120 s, as well as their corresponding projections. Magenta color is used 
for LineO, blue for Linel, and black for Line2. Red circles are the projection on the 
box walls of a sphere with radius 6 = 17.5 Mm. 
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The limit of Equation (4.1) when z — )■ corresponds to the magnetic field com- 
ponents in the plane of the current sheet, and it is equivalent to the one described in 
the previous chapter by Equation (3.1) with Bg = Bey, and uniform guide field B^x- 
Therefore, the axes of the initial three-dimensional tubes that cross the reconnection 
region on each side of the current sheet are well described by the two-dimensional 
field lines (3.3), at least while the tubes are close to the CS plane. The tubes have 
legs that do not lie on this plane and have footpoints at the base of the arcade, but 
their retractions occur mostly while the field lines lie within the plane of the current, 
until they arrive to the top of the arcade, as described in the previous chapter. The 
left panel of Figure 34 shows the evolution of LineO projected on the plane of the 
current sheet. The solid (dotted) black line is the theoretical initial field line on the 
front (back) side of the current sheet (Equation (3.3) for a reconnection location 
and Co = 45°). Each color represents a different time as described in the caption. 
Circles indicate the positions of the bends on the front side of the current sheet. 

To be compared with the above three-dimensional simulations, we run DEFT sim- 
ulations with the following parameters: reconnection location y'^ = 2/3, reconnection 
angle Ci? = 37°, and plasma-/? there equal to 0.54. We implement temperature- 
independent viscosity and thermal conduction. The Reynolds number of this simula- 
tion is 300, and the Prandtl number is 1. We divided the tube in 1000 segments. We 
will present our results in the same units used for the three-dimensional simulation. 
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Figure 33. LineO retraction from the reconnection site. The colored sohd hnes are the 
field line curves for times at 90 s intervals: 180 s (red), 270 s (green), 360 s (blue), 540 
s (yellow), 630 s (magenta), 720 s (cyan), and 810 s (black). Colored dashed lines are 
the field line projections on the box walls of the same color three-dimensional curve. 
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The right panel of Figure 34 shows the reconnected flux tube retraction for the DEFT 
simulation. The initial shape of the tube is similar to the reconnected shape of LineO. 

As times goes by, the concavity of the DEFT tube is less pronounced than in the 
run for LineO. This appears to be due to a smaller initial velocity of the center of 
the tube in the three-dimensional case and could be related to the finite time it takes 
to create the reconnected tube by the diffusion process. In the DEFT simulations, 
reconnection is assumed to happen instantaneously. The bends in both cases move 
at the local Alfven speed. The circles in the right panel of Figure 34 were drawn 
assuming they do so, and they coincide with the actual bends. The top panel of 
Figure 35 shows the x- and y-positions of the bend on the front side of the current 
sheet, as well as the theoretical motion they should follow if the move at the local 
Alfven speed along the legs of the field line (black dashed lines). 

The delay in the motion of the central part of the tube can be seen in the middle 
panel of Figure 35. Magenta color represents the downward motion of the center of 
the tube for LineO, and the cyan color shows the motion of the DEFT tube. The 
dashed lines represent quadratic fits. Their accelerations are similar, but the initial 
speed for LineO is significantly smaller than the one for the DEFT tube. In our thin 
flux tube model, the shape of the tube at the reconnection position has a sharp angle 
and it immediately retracts after being reconnected at Alfvenic speed. The center 
part of LineO stays in the top part of the current sheet longer than the DEFT tube 
does. This has some consequences on the parallel dynamics since the plasma responds 
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Figure 34. Current sheet plane. Both panels share the same y-axis. Left panel: MHD 
LineO evolution projection on the current sheet plane. Solid (dotted) black line is the 
theoretical initial field line on the front (back) side of the current sheet. Orange 
dashed field line corresponds to t = s and violet dashed field line corresponds to 
t = 90 s. The rest of the colored dashed lines correspond to the same times and 
colors as in Figure 33. Circles indicate the positions of the bends for each field line 
on the front side of the current sheet. No circles have been plotted for the initial field 
line. Crosses indicate the locations of the gas-dynamic shocks. Right panel: DEFT 
fiux tube simulation. Each line represents the shape of the tube for the same times 
as in the left panel. Circles are the positions of the bends and crosses indicate the 
locations of the gas-dynamic shocks. 
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to the changes in background magnetic pressure, as described in the previous chapter. 
Despite the mentioned differences, the retraction is surprisingly similar for both cases. 

4.1.1. Parallel Dynamics 

The comparison between the parallel dynamics for both cases can be seen in 
Figures 36 and 37. The first figure shows the evolution of the electron density for both 
cases. The left panel corresponds to LineO. The central density depression is caused 
by the initial reconnection episode. Ohmic heating from the enhanced resistivity 
creates a hot spot at the center of the tube that cannot be efficiently smoothed 
out due to the lack of thermal conduction in the ARMS computer program. Initial 
pressure equilibrium is quickly established through sound waves, leaving a density 
void in the center of the tube that persists for a long time. This is an artificial effect 
of the three-dimensional simulations. We hypothesize that this hot spot will quickly 
disappear in a real plasma due to thermal conduction. The thin flux tube model does 
not present this behavior since the effects of the ohmic heating are assumed to have 
already ceased by the time the simulation starts. 

Sharp discontinuities in the electron density can be seen moving outward from the 
center of the tube. These are the gas-dynamic shocks. Crosses show their locations 
for each time, also shown in Figure 34. The density jumps by approximately 40%. 
The post-shock density increases while the center part of the tube is in the top part 
of the current sheet and decreases while this region moves through the lower half of 
the plane (cyan and black colors). This is as expected from our conclusions from the 
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Figure 35. LineO speed comparison. All panels share the same horizontal axis (time 
in seconds). Top panel: the green symbols (left coordinate axis) show the x-position 
of the bend on the front side of the current sheet as time goes by for the MHD 
simulation of LineO. The dashed black line close to the green symbols represents the 
theoretical x-positions of the bend if it were moving along the legs at the local Alfven 
speed. The magenta symbols (right coordinate axis) show the ^/-position of the bend 
on the front side of the current sheet as time goes by for the MHD simulation of 
LineO. The dashed black line close to the magenta symbols represents the theoretical 
^/-positions of the bend if it were moving along the legs at the local Alfven speed. 
Times before approximately t = 150 s are not shown since the position of the bends 
are not well-defined while the tube reconnects. Middle panel: the magenta symbols 
show the y-position of the center of the tube as time goes by for the MHD simulation 
of LineO. The dashed black line close to the magenta symbols represents a quadratic 
fitting of the magenta curve, y{t) = 65 Mm —0.08 Mm s"^ t — 4.8 x 10~^ Mm s~^ t^. 
The magenta symbols show the ^/-position of the center of the tube as time goes by 
for the DEFT simulation. The dashed black line close to the cyan symbols represents 
a quadratic fitting of the cyan curve, y{t) = 66 Mm —0.13 Mm s"^ t — 5.9 x 10~^ Mm 
s~^ t"^. Bottom panel: the magenta symbols show the length of the post-shock region, 
/, as time goes by for the MHD simulation of LineO. The dashed black line close to 
the magenta symbols represents a linear fitting of the magenta curve, y{t) = 3.14 
Mm +0.23 Mm s"^ t. The cyan symbols show the length of the post-shock region, /, 
as time goes by for the DEFT simulation. The dashed black line close to the cyan 
symbols represents a linear fitting of the cyan curve, y{t) = 10.0 Mm +0.27 Mm s"^ 
t. 
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previous chapter. A decrease in density can also be seen following the bends (circles 
in the figure) for late times as the bends move through the lower part of the current 
sheet. Due to the increased concavity of the tube, this decrease in density is not as 
pronounced as in its DEFT counterpart. 

The electron density evolution for the DEFT simulation is shown in the right 
panels of the figure. The jump in density is very similar to the one for the MHD 
simulation. Its value increases as the central part of the tube moves through the top 
half of the current sheet and decreases during its motion through the lower half of the 
plane. The density voids after the bends are more pronounced in this case than in 
its three dimensional counterpart since the sections of the tube following the bends 
move through regions of the current sheet with lower magnetic pressure. Although 
the MHD simulations were stopped before the field lines arrived to the top of the 
arcade, we expect them to come to a rest at the top of the arcade. Therefore, the 
sections after the bends will eventually decrease their density as they move through 
the bottom half of the current sheet. The slowing of the tubes on top of the arcade 
for three-dimensional simulations has been shown by Linton et al. (2009), although 
no analysis of the parallel dynamics was attempted in that paper. 

The hot spot at the center of the tube for the three-dimensional simulations 
can be seen in the left panel of Figure 37. This panel shows temperature profiles 
for different times. The right panel shows the temperature profiles for the DEFT 
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Figure 36. LineO electron density comparison. Both panels share the same y-axis. 
Left panel: each solid curve represents the electron density at a different time for 
LineO. Colors represent the same times as in Figure 34. Circles and crosses are also 
the same as for the left panel of Figure 34. Right panel: each solid curve represents 
the electron density at a different time for the DEFT simulation. Circles and crosses 
correspond to the same ones as in the right panel of Figure 34. 
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simulation. The jump in temperature across the gas-dynamic shocks is comparable 
for both simulations. 

MHD simulation DEFT Simulation 



1 .6 
1 .4 



1.2 

1 .0 

-150 -100 -50 50 100 150-150 -100 -50 50 100 150 
X (Mm) X (Mm) 

Figure 37. LineO temperature comparison. Both panels share the same y-axis. Left 
panel: each solid curve represents the temperature at a different time for LineO. Colors 
represent the same times as in Figure 34. Circles and crosses are also the same as for 
the left panel of Figure 34. Right panel: each solid curve represents the temperature 
at a different time for the DEFT simulation. Circles and crosses correspond to the 
same ones as in the right panel of Figure 34. 

The speed of the shocks along the tubes is shown for both cases in the bottom 
panel of Figure 35. The magenta symbols show the LineO post-shock region length, 
/, as function of time. For each time shown in Figure 37, this corresponds to the 
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distance between crosses. The cyan color corresponds to the post-shock region length 
for the DEFT simulation. The speeds are similar. The slight difference could be 
attributed to the fact that the shock speed is related to the pre-shock temperature. 
These temperatures are different for each case since the regions after the bends have 
different temperatures, depending on where in the current sheet these sections are 
moving through. 

LineO behaves like a typical field line from a flux tube that reconnected at the high 
resistivity sphere. If all the field lines that cross the high resistivity sphere followed 
this behavior, the DEFT simulation would approximately reproduce their evolution. 
The evolutions of Linel and Line2 are shown in Figures 38 and 39, respectively. 
Apparently, for these field lines the disturbance caused by the initial reconnection 
event spontaneously excites patchy reconnection events later in the simulation. This 
can be seen in the artificial motion of the footpoint on the z < region. The source 
of this secondary reconnection may be numerical resistivity (Linton et al., 2009). We 
do not know if this would occur in a real plasma. 

4.2. Conclusions 

We have presented a comparison between three-dimensional MHD simulations 
and DEFT simulations. It seems that the assumptions and approximations of our 
thin flux tube model are not too restrictive and describe in reasonable terms the 
behavior of a flux tube that has been reconnected in a patch, as long as there are 
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Figure 38. Linel retraction from the reconnection site. The colored sohd hnes are the 
field line curves for times at 90 s intervals: 180 s (red), 270 s (green), 360 s (blue), 540 
s (yellow), 630 s (magenta), 720 s (cyan), and 810 s (black). Colored dashed lines are 
the field line projections on the box walls of the same color three-dimensional curve. 
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Figure 39. Line2 retraction from the reconnection site. The colored sohd hnes are the 
field line curves for times at 90 s intervals: 180 s (red), 270 s (green), 360 s (blue), 540 
s (yellow), 630 s (magenta), 720 s (cyan), and 810 s (black). Colored dashed lines are 
the field line projections on the box walls of the same color three-dimensional curve. 
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not secondary reconnections along the tube. The non-trivial assumption that a small 
enhanced diffusion region will disconnect field lines from one side of the current sheet 
and connect them to field lines on the other side is satisfied for the MHD simulations. 

As predicted, reconnected field lines retract from the reconnection sphere in an 
almost two-dimensional motion parallel to the current sheet, and gas-dynamic shocks 
move along the field lines, heating and compressing the plasma. The jump conditions 
across the shock predicted by our model seem to be accurate. 

There are some expected differences in the field line shapes and parallel dynamics 
between the simulations, but they are surprisingly small, given the simplifications of 
our model. The main differences are a consequence of the initial velocity of the center 
of the tube. In three-dimensional simulations, the formation of the tube is delayed 
during the finite amount of time that reconnection takes to occur, and therefore the 
center region of the field lines spends a longer time on the top half of the current sheet. 
Since the bends move at the local Alfven speed along the field lines, the concavity of 
the field lines is much larger than for the DEFT simulations. 

We have only presented one case for the comparison between our model and MHD 
simulations. The chosen conditions are quite generic, and we expect agreement also 
for other cases. We hope that in the near future we can extend our analysis to more 
points in the parameter space. 

The thin flux tube equations are simpler than the fully three-dimensional MHD 
equations. This has allowed us to predict approximately the strength of the shocks 
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generated by reconnection. The thickness of the shocks, determined by the trans- 
port coefficients, is non-trivial for solar coronal conditions. This thickness can be 
calculated using our model more accurately than through three-dimensional simu- 
lations since anisotropic, temperature-dependent viscosity and thermal conduction 
cannot be easily implemented in three-dimensional MHD simulations. In addition, 
due to the simplicity of our DEFT code, the exploration of the parameter space is 
less computationally-expensive and time-consuming. 
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5. CONCLUSIONS 

The foregoing chapters have given an account of several investigations and devel- 
opments into the dynamics of magnetic loops following their formation by magnetic 
reconnection. 

We have not only refined and improved a novel, time-dependent patchy reconnec- 
tion model, but also presented general equations that govern the evolution of thin flux 
tubes. The ideal part of these equations could be applied to model thin flux tubes 
in several non-relativistic astrophysical situations. For example, there is evidence for 
localized and transient nature of magnetic reconnection in the magnetosphere. The 
first such evidence was found in so-called flux transfer events identified in situ mag- 
netospheric observations (Russell & Elphic, 1978). There, reconnection starts and 
stops in a few minutes resulting on the ripping off of flux tubes. It is also believed 
that magnetic fields generated by the solar dynamo are convectively collapsed into 
coherent flux tubes (through processes which are still not well understood), become 
buoyant in the solar convection zone, and consequently erupt as active regions (Schus- 
sler et al., 1994; Nandy, 2006). If non-ideal processes are unimportant, our equations 
could be used to study the mentioned phenomena. 

Through our investigations, we have learned that even for plasmas that are 
strongly magnetized (like the ones in the lower solar corona), where thermal pres- 
sure is much smaller than magnetic pressure (low plasma-/3), the pressure-driven 
parallel term in the thin flux tube momentum equation cannot be neglected since 
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magnetic forces are only perpendicular to the magnetic field. This term modulates 
tube's acoustics and prevents unlimited mass pile up inside the tubes. 

Standard models of reconnection assume perfectly anti-parallel magnetic fields. 
We believe this to be too restrictive for the solar corona, where photospheric motions 
produce coronal fields with varied directions. In addition, the observation of transient 
and localized features, like SADs and SADLs, and possibly individual loops (like the 
ones seen with TRACE) seem to indicate that, at least for these cases, reconnection 
occurs in short-lived patches. The traditional steady-state, two- or 2. 5- dimensional 
reconnection scenarios are not adequate to describe them. 

We believe our work has shown that the simple addition of a guide field extends 
considerably the range of possible retracting tube geometries. If only anti-parallel 
fields are assumed, the shape for post-reconnection outflows corresponds to coUi- 
mated jets. On the other hand, patchy reconnection with skewed fields not only 
accommodate jet geometries (if the current sheet is observed edge-on), but also loop- 
shape geometries (if the plane of the sky is almost parallel to the current sheet plane). 
SADs and SADLs seem to be examples of each case, respectively. 

The parallel dynamics inside the tubes are as rich as the perpendicular dynamics. 
Super-sonic parallel flows are generated at the RDs that move along the legs of the 
tubes, "transmitting" at the local Alfven speed the information of what has happened 
at the reconnection location (change in topology). The super-sonic flows collide at the 
center and result in gas-dynamic shocks that also move along the legs of the tubes. 
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heating and compressing the plasma. This generates two long-hved (confined by the 
infiows) and hot plasma plugs that move outwardly in opposite directions from the 
reconnection region. 

The reconnection patch creates two V-shaped reconnected tubes that shorten 
as they retract in opposite directions, due to magnetic tension. One of them moves 
upward toward the top edge of the CS, and the other one moves downward toward the 
top of the underlying arcade. The plasma inside each tube is far from equilibrium and 
the densities achieved in the plasma plugs can be much higher than the steady-state 
Rankine-Hugoniot post-shock values. This is only possible for very strong shocks 
and time-dependent models like the one we have presented in this dissertation. The 
strength of the shocks is determined by the Mach number of the parallel infiows 
inside the tubes (Equation (2.43)). For strongly magnetized plasmas (small plasma- 
P), the jump in density across the shocks may be more than one order of magnitude, 
depending on the half reconnection angle value. 

We have also included realistic non-ideal terms in our equations, namely thermal 
conduction and viscous momentum exchange. We restricted ourselves to strongly 
magnetized plasmas, for which the magnetic field confines the transport processes to 
be mostly aligned with its direction. For fully ionized plasmas, like the ones found 
in the solar corona, the transport coefficients have a strong temperature dependence 
(^ T^/^). The addition of these diffusive processes has proven to be essential for a 
self- consisting generation of the gas-dynamic shocks. 
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It is generally believed that shocks require a kinetic treatment. We have shown 
that a simple fluid treatment is adequate whenever the diffusive coefficient's growth 
with temperature is faster than the particle's mean free path growth with tempera- 
ture. This condition is satisfied for fully ionized plasmas where Coulomb potentials 
determine the particle mean free paths (~ T^) and Spitzer's transport coefficients 
(~ T^/^) are used. In this case, the thickness of the shocks exceeds the particle mean 
free paths by at least an order of magnitude. When temperature independent trans- 
port coefficients are used, the results are several orders of magnitude smaller than the 
corresponding mean free paths, which in the past created some concern about the use 
of ffuid equations for shock studies. 

The validity of assuming ffuid equations, as opposed to kinetic theory, may also 
be a concern if the ratio between the diffusion heat ffux and the value for streaming 
particles is larger than a certain value (Campbell, 1984). We tested the collisionality of 
the plasma in our simulations, and obtained ratios that support the use of continuum 
ffuid equations. 

We theoretically predicted the thickness of the shocks, as well as their internal 
structure. Our analysis is novel in its focus on solar ffare conditions, including very 
low Prandtl numbers (ratio of viscosity to thermal conductivity). In this case, the 
shocks present a relatively thin isothermal sub-shock where most of the plasma com- 
pression occurs, and thermal fronts that can be as long as the entire tubes and where 
temperature increases considerably. These fronts move along the legs of the tube and 
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may be relevant to chromospheric evaporation, although we have not yet analyzed 
this effect. Our model can be used to couple the reconnection regime high in the solar 
corona with evaporation dynamics at the footpoints. 

Traditional models have not been able to explain why some of the observed SADs 
appear to be hot and relatively devoid of plasma. We have shown that plasma deple- 
tion naturally occurs in flux tubes that retract along regions of decreasing external 
magnetic pressure. We have demonstrated this effect on non-uniform Syrovatskii-type 
CSs with skewed magnetic fields, that seem to be ubiquitous in the solar corona. In 
this kind of CSs, the background magnetic pressure has its maximum at the center of 
the CS plane and decreases toward its edges. The hot temperature of the SADs can 
be explained by the long thermal fronts that heat the plasma in regions of density 
depletion. 

We have modeled the interaction of the DOWNWARD moving tubes with the 
underlying arcade by means of a simple damping force. A DOWNWARD moving tube 
will arrive at the top of the arcade that will slow it to a stop. Here, the perpendicular 
dynamics is halted, but the parallel dynamics continues along its legs; the RDs are 
shut down, and the gas is rarefied to even lower densities. After the tubes lie on top of 
the arcade, rarefaction waves continue decreasing plasma density. The hot post-shock 
regions continue evolving, determining a long lasting hot region on top of the arcade. 
We have also provided an observational method based on total emission measure and 
mean temperature, that indicates where in the CS the tube has been reconnected. 



159 

The DEFT computer program we developed solves the thin flux tube equations 
and simulates the retraction of post-reconnection tubes for a wide range of back- 
ground conditions. The program seems to be very robust and resolves the gas- 
dynamic shocks efficiently. Its output can be used as input for simulations that 
deal with chromospheric evaporation (for DOWNWARD moving tubes) and with so- 
lar ejecta (UPWARD moving tubes). A comparison of DEFT simulations with fully 
three-dimensional magnetohydrodynamic simulations has been presented to assess 
the validity of the thin flux tube model. The seemingly restrictive thin flux tube 
assumptions reproduce the general behavior of the three-dimensional plasma. 

The correct answer to the disparity between observations and models regarding 
the speed of localized reconnection outflows is still not known, but most likely the 
slowing down is due to the interaction between tubes and their surroundings, which we 
have not considered in the present work. Here, the reconnected tubes are assumed to 
be completely isolated from their background, therefore we have restricted ourselves 
to standard reconnection scenarios where the outflows are Alfvenic. We do not claim 
that we have explained all SADs' observables, but our model seems to support the 
idea that they are a by-product of reconnection. 

There are many improvements to the model that can be incorporated and we 
hope that this work has proven its usefulness and opens the door to such additions. 



160 



REFERENCES CITED 



161 

Aly, J. J., & Amari, T. 1989, A & A, 221, 287 

Asai, A., Yokoyama, T., Shimojo, M., & Shibata, K. 2004, ApJ, 605, L77 
Aschwanden, M. J. 1999, Sol. Phys., 190, 233 

— . 2005, Physics of the Solar Corona. An Introduction with Problems and Solutions 
(2nd edition) (Praxis Publishing Ltd., Chichester, UK; Springer, New York, Berlin) 

Biernat, H. K., Heyn, M. F., & Semenov, V. S. 1987, J. Geophys. Res., 92, 3392 

Biernat, H. K., Semenov, V. S., & Rijnbeek, R. P. 1998, J. Geophys. Res., 103, 4693 

Birn, J., & Priest, E. R. 2007, Reconnection of magnetic fields : magnetohydrody- 
namics and coUisionless theory and observations (Cambridge University Press) 

Birn, J., et al. 2001, J. Geophys. Res., 106, 3715 

Biskamp, D. 1986, Physics of Fluids, 29, 1520 

Biskamp, D., & Schwarz, E. 2001, Physics of Plasmas, 8, 4729 

Braginskii, S. I. 1965, Reviews of Plasma Physics, 1, 205 

Campbell, P. M. 1984, Phys. Rev. A, 30, 365 

Cargill, P. J., Mariska, J. T., & Antiochos, S. K. 1995, ApJ, 439, 1034 
Carmichael, H. 1964, NASA Special Publication, 50, 451 

Caspi, A. 2010, PhD thesis. Department of Physics, University of California, Berkeley, 
CA 94720-7450, USA 

Caspi, A., & Lin, R. P. 2010, ApJ, 725, L161 

Cassak, P. A., Shay, M. A., & Drake, J. F. 2005, Physical Review Letters, 95, 235002 

Christensen-Dalsgaard, J. 2002, Reviews of Modern Physics, 74, 1073 

Colburn, D. S., & Sonett, C. P. 1966, Space Sci. Rev., 5, 439 

Cook, J. W., Cheng, C, Jacobs, V. L., & Antiochos, S. K. 1989, ApJ, 338, 1176 

Coroniti, F. V. 1970, Journal of Plasma Physics, 4, 265 

Courant, R., Friedrichs, K., & Lewy, H. 1967, IBM Journal of Research and Devel- 
opment, 11, 215 

Craig, I. J. D. 2008, A & A, 487, 1155 

Craig, I. J. D., & Litvinenko, Y. E. 2009, A & A, 501, 755 



162 

Craig, I. J. D., Litvinenko, Y. E., & Senanayake, T. 2005, A & A, 433, 1139 

Crosby, N. B., Aschwanden, M. J., & Dennis, B. R. 1993, Sol. Phys., 143, 275 

Devore, C. R. 1991, Journal of Computational Physics, 92, 142 

Domingo, V., Fleck, B., & Poland, A. I. 1995, Sol. Phys., 162, 1 

Drake, J. F., Swisdak, M., Schoeffler, K. M., Rogers, B. N., & Kobayashi, S. 2006, 
Geophys. Res. Lett., 33, 13105 

Edmondson, J. K., Antiochos, S. K., DeVore, C. R., & Zurbuchen, T. H. 2010, ApJ, 
718, 72 

Emslie, A. C, Dennis, B. R., Holman, G. D., & Hudson, H. S. 2005, Journal of 
Geophysical Research (Space Physics), 110, A11103 

Emslie, A. G., et al. 2004, Journal of Geophysical Research (Space Physics), 109, 
A10104 

Erkaev, N. V., Semenov, V. S., & Jamitzky, F. 2000, Physical Review Letters, 84, 
1455 

Feldman, U., Schiihle, U., Widing, K. G., & Laming, J. M. 1998, ApJ, 505, 999 

Forbes, T. G. 1986, ApJ, 305, 553 

— . 2000, J. Geophys. Res., 105, 23153 

Forbes, T. G., & Acton, L. W. 1996, ApJ, 459, 330 

Forbes, T. G., & Isenberg, P. A. 1991, ApJ, 373, 294 

Forbes, T. G., & Malherbe, J. M. 1986, ApJ, 302, L67 

Forbes, T. G., Malherbe, J. M., & Priest, E. R. 1989, Sol. Phys., 120, 285 

Forbes, T. G., & Priest, E. R. 1983, Sol. Phys., 84, 169 

Forbes, T. G., et al. 2006, Space Sci. Rev., 123, 251 

Fowles, G. R. 1975, Physics of Fluids, 18, 776 

Furth, H. P., Killeen, J., & Rosenbluth, M. N. 1963, Physics of Fluids, 6, 459 

Gallagher, P. T., Dennis, B. R., Krucker, S., Schwartz, R. A., & Tolbert, A. K. 2002, 
Sol. Phys., 210, 341 

Gilbarg, D., & Paolucci, D. 1953, Indiana Univ. Math. J., 2, 617 



163 

Goldstein, H. 1959, Classical mechanics (Addison- Wesley Publishing Company, Inc., 
Reading, Mass., U.S. A, London, England) 

Golub, L., et al. 1999, Physics of Plasmas, 6, 2205 

Gosling, J. T., Asbridge, J. R., Bame, S. J., & Feldman, W. C. 1976, J. Geophys. 
Res., 81, 5061 

Grad, H. 1952, Communications on pure and applied mathematics, v, 257 

Green, R. M. 1965, in lAU Symposium, Vol. 22, Stellar and Solar Magnetic Fields, 
ed. R. Lust, 398-+ 

Guidoni, S. E., & Longcope, D. W. 2010, ApJ, 718, 1476 

— . 2011, ApJ, 730, 90 

Handy, B. N., et al. 1999, Sol. Phys., 187, 229 

Harris, E. 1962, II Nuovo Cimento (1955-1965), 23, 115, 10.1007/BF02733547 

Heyn, M. F., Biernat, H. K., Rijnbeek, R. P., & Semenov, V. S. 1988, Journal of 
Plasma Physics, 40, 235 

Hirayama, T. 1974, Sol. Phys., 34, 323 

Holman, G. D., Sui, L., Schwartz, R. A., & Emslie, A. G. 2003, ApJ, 595, L97 

Hoyt, D. v., & Schatten, K. H. 1998a, Solar Physics, 179, 189, 
10.1023/A:1005007527816 

— . 1998b, Solar Physics, 181, 491, 10.1023/A:1005056326158 

Hu, Y. Q., & Low, B. C. 1982, Sol. Phys., 81, 107 

Hudson, H. S., & McKenzie, D. E. 2001, Earth, Planets, and Space, 53, 581 
Hugoniot, H. 1887, Journal de I'Ecole Polytechnique, 57, 3 
Innes, D. E., McKenzie, D. E., & Wang, T. 2003a, Sol. Phys., 217, 267 
— . 2003b, Sol. Phys., 217, 247 

Jiang, Y. W., Liu, S., Liu, W., & Petrosian, V. 2006, ApJ, 638, 1140 

Katz, N., Egedal, J., Fox, W., Le, A., Bonde, J., & Vrublevskis, A. 2010, Physical 
Review Letters, 104, 255004 

Kennel, C. F. 1988, J. Geophys. Res., 93, 8545 



164 



Klimchuk, J. A. 1997, in Astronomical Society of the Pacific Conference Series, 
Vol. Ill, Magnetic Reconnection in the Solar Atmosphere, ed. R. D. Bentley & 
J. T. Mariska, 319-+ 

Klimchuk, J. A. 2001, Space Weather (Geophysical Monograph 125), ed. P. Song, 
H. Singer, G. Siscoe (Washington: Am. Geophys. Un.), 143 (2001), 125, 143 

Kopp, R. A., & Pneuman, G. W. 1976, Sol. Phys., 50, 85 

Krucker, S., Hurford, G. J., & Lin, R. P. 2003, ApJ, 595, L103 

Krucker, S., & Lin, R. P. 2008, ApJ, 673, 1181 

Landau, L. D., & Lifshitz, E. M. 1959, Fluid mechanics (Oxford: Pergamon Press, 
1959) 

Lin, J., Cranmer, S. R., & Farrugia, C. J. 2008, Journal of Geophysical Research 
(Space Physics), 113, 11107 

Lin, J., & Forbes, T. G. 2000, J. Geophys. Res., 105, 2375 

Lin, J., Ko, Y., Sui, L., Raymond, J. C., Stenborg, G. A., Jiang, Y., Zhao, S., & 
Mancuso, S. 2005, ApJ, 622, 1251 

Lin, R. P., & Hudson, H. S. 1971, Sol. Phys., 17, 412 

— . 1976, Sol. Phys., 50, 153 

Lin, R. P., Schwartz, R. A., Felling, R. M., & Hurley, K. C. 1981, ApJ, 251, L109 

Linton, M. G., Devore, C. R., & Longcope, D. W. 2009, Earth, Planets, and Space, 
61, 573 

Linton, M. G., & Longcope, D. W. 2006, ApJ, 642, 1177 
Litvinenko, Y. E. 2005, Sol. Phys., 229, 203 
Longcope, D. W. 2001, Physics of Plasmas, 8, 5277 
Longcope, D. W., & Bradshaw, S. J. 2010, ApJ, 718, 1491 

Longcope, D. W., Des Jardins, A. C., Carranza-Fulmer, T., & Qiu, J. 2010, Sol. 
Phys., 267, 107 

Longcope, D. W., Guidoni, S. E., & Linton, M. G. 2009, ApJ, 690, L18 

Longcope, D. W., McKenzie, D. E., Cirtain, J., & Scott, J. 2005, ApJ, 630, 596 

Low, B. C. 2003, in Astronomical Society of the Pacific Conference Series, Vol. 
286, Current Theoretical Models and Future High Resolution Solar Observations: 
Preparing for ATST, ed. A. A. Pevtsov & H. Uitenbroek, 335 — h 



165 

Ma, Z. W., & Bhattacharjee, A. 2001, J. Geophys. Res., 106, 3773 

MacNeice, P., Olson, K. M., Mobarry, C, de Fainchtein, R., & Packer, C. 2000, 
Computer Physics Communications, 126, 330 

Mandrini, C. H. 2010, in lAU Symposium, Vol. 264, lAU Symposium, ed. A. G. Koso- 
vichev, A. H. Andrei, & J.-P. Roelot, 257-266 

Martens, P. C. H., Kankelborg, C. C, & Berger, T. E. 2000, ApJ, 537, 471 

Masuda, S., Kosugi, T., Hara, H., Tsuneta, S., & Ogawara, Y. 1994, Nature, 371, 495 

McKenzie, D. E. 2000, Sol. Phys., 195, 381 

McKenzie, D. E., & Hudson, H. S. 1999, ApJ, 519, L93 

— . 2001, Earth, Planets, and Space, 53, 577 

McKenzie, D. E., & Savage, S. L. 2009, ApJ, 697, 1569 

Nandy, D. 2006, Journal of Geophysical Research (Space Physics), 111, A12S01 
Narukage, N., & Shibata, K. 2006, ApJ, 637, 1122 

National Research Council. 1996, Radiation Hazards to Crews on Interplanetary Mis- 
sions: Biological Issues and Research Strategies, editor: National Academy Press, 
Washington, DC 

Nishizuka, N., Takasaki, H., Asai, A., & Shibata, K. 2010, ApJ, 711, 1062 

Nitta, S., Tanuma, S., & Maezawa, K. 2002, ApJ, 580, 538 

Norquist, D. 2011, Solar Physics, 269, 111, 10.1007/sll207-010-9689-y 

Ogawara, Y., Takano, T., Kato, T., Kosugi, T., Tsuneta, S., Watanabe, T., Kondo, 
I., & Uchida, Y. 1991, Sol. Phys., 136, 1 

Parker, E. N. 1957, J. Geophys. Res., 62, 509 

— . 1958, ApJ, 128, 664 

— . 1963, ApJS, 8, 177 

Petschek, H. E. 1964, NASA Special Publication, 50, 425 
Petschek, H. E., & Thorne, R. M. 1967, ApJ, 147, 1157 
Priest, E. R. 1982, Geophysics and Astrophysics Monographs, 21 
— . 1985, Reports on Progress in Physics, 48, 955 



166 

Priest, E. R., & Raadu, M. A. 1975, Sol. Phys., 43, 177 

Priest, E. &Forbes, T., ed. 2000, Magnetic Reconnection (Cambridge University 
Press) 

Pudovkin, M. I., & Semenov, V. S. 1985, Space Science Reviews, 41, 1 
Rankine, W. J. M. 1870, Phil. Trans. R. Soc. Lond., 160, 277 
Rosner, R., Tucker, W. H., & Vaiana, G. S. 1978, ApJ, 220, 643 
Russell, C. T., & Elphic, R. C. 1978, Space Science Reviews, 22, 681 
Savage, S. L., & McKenzie, D. E. 2011, ApJ, 730, 98 

Savage, S. L., McKenzie, D. E., Reeves, K. K., Forbes, T. G., & Longcope, D. W. 
2010, ApJ, 722, 329 

Schettino, G., Poletto, G., & Romoli, M. 2010, ApJ, 708, 1135 

Scholer, M., & Roth, D. 1987, J. Geophys. Res., 92, 3223 

Schrijver, C. J. 2009, Advances in Space Research, 43, 739 

Schrijver, C. J., et al. 1999, Sol. Phys., 187, 261 

Schrijver, C. J. & Siscoe, G. L., ed. 2010, Heliophysics: Space Storms and Radiation: 
Causes and Effects (Cambridge University Press) 

Schussler, M., Caligari, P., Ferriz-Mas, A., & Moreno-Insertis, F. 1994, A & A, 281, 
L69 

Seaton, D. B., & Forbes, T. G. 2009, ApJ, 701, 348 

Semenov, V. S., Heyn, M. F., & Kubyshkin, I. V. 1983, Soviet Astronomy, 27, 660 

Semenov, V. S., Kubyshkin, I., Levedeva, V. V., Sidneva, M. V., Biernat, H. K., 
Heyn, M. F., Besser, B. P., & Rijnbeek, R. P. 1992, J. Geophys. Res., 97, 4251 

Semenov, V. S., Volkonskaya, N. N., & Biernat, H. K. 1998, Physics of Plasmas, 5, 
3242 

Sheeley, N. R., Walters, J. H., Wang, Y., & Howard, R. A. 1999, J. Geophys. Res., 
104, 24739 

Sheeley, Jr., N. R., & Wang, Y. 2002, ApJ, 579, 874 
— . 2007, ApJ, 655, 1142 

Sheeley, Jr., N. R., Warren, H. P., & Wang, Y. 2004, ApJ, 616, 1224 



167 

Shibata, K., Masuda, S., Shimojo, M., Hara, H., Yokoyama, T., Tsuneta, S., Kosugi, 
T., & Ogawara, Y. 1995, ApJ, 451, L83+ 

Shimizu, T., Kondo, K., Ugai, M., & Shibata, K. 2009, ApJ, 707, 420 

Shu, F. 1991, Physics of Astrophysics, Vol. II: Gas Dynamics (University Science 
Books, 648 Broadway, Suite 902, New York, NY 10012) 

Skender, M., Vrsnak, B., & Martinis, M. 2003, Phys. Rev. E, 68, 046405 

Soward, A. M. 1982, Journal of Plasma Physics, 28, 415 

Soward, A. M., & Priest, E. R. 1982, Journal of Plasma Physics, 28, 335 

Spitzer, L. 1962, Physics of Fully Ionized Gases (New York: Interscience (2nd edition)) 

Spruit, H. C. 1981, A & A, 98, 155 

Sturrock, P. A. 1966, Nature, 211, 695 

Sturrock, P. A. 1968, in lAU Symposium, Vol. 35, Structure and Development of 
Solar Active Regions, ed. K. O. Kiepenheuer, 471 — h 

Sweet, P. A. 1958, II Nuovo Cimento SuppL, 8, 188196 

Syrovatskii, S. I. 1971, Soviet Journal of Experimental and Theoretical Physics, 33, 
933 

Thomas, L. H. 1944, J. Chem. Phys., 12, 449 
Tomczak, M. 2009, A & A, 502, 665 

Tsuneta, S., Hara, H., Shimizu, T., Acton, L. W., Strong, K. T., Hudson, H. S., & 
Ogawara, Y. 1992, PASJ, 44, L63 

Tsuneta, S., Masuda, S., Kosugi, T., & Sato, J. 1997, ApJ, 478, 787 

Tsuneta, S., et al. 1991, Sol. Phys., 136, 37 

Ugai, M., & Tsuda, T. 1977, Journal of Plasma Physics, 17, 337 

Svestka, Z., Farm'k, F., Hudson, H. S., & Hick, P. 1998, Sol. Phys., 182, 179 

van Allen, J. A. 1981, in Space Science Comes of Age: Perspectives in the History of 
the Space Sciences, ed. P. A. Hanle & V. D. Chamberlain, 58-73 

Vourlidas, A., Subramanian, P., Dere, K. P., & Howard, R. A. 2000, ApJ, 534, 456 

Wang, T., Sui, L., & Qiu, J. 2007, ApJ, 661, L207 



168 



Wang, T. J., Solanki, S. K., Innes, D. E., & Curdt, W. 2002, in ESA Special Publica- 
tion, Vol. 505, SOLMAG 2002. Proceedings of the Magnetic Coupling of the Solar 
Atmosphere Euroconference, ed. H. Sawaya-Lacoste, 607-610 

Warren, H. P., Bookbinder, J. A., Forbes, T. G., Golub, L., Hudson, H. S., Reeves, 
K., & Warshall, A. 1999, ApJ, 527, L121 

Withbroe, G. L., & Noyes, R. W. 1977, ARA&A, 15, 363 

Yamada, M., Kulsrud, R., & Ji, H. 2010, Reviews of Modern Physics, 82, 603 

Yokoyama, T., Akita, K., Morimoto, T., Inoue, K., & Newmark, J. 2001, ApJ, 546, 
L69 



169 



APPENDIX A 
LIST OF ABBREVIATIONS 
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Al 


— Aluminum 


AlMg 

o 


— Aluminum Magnesium 


ARMS 


— Adaptively Refined Magnetohydrodynamics Solver (numeri- 




cal MHD solver) 


CME 


- Coronal Mass Ejection 


CS 


- Current Sheet 


CSHKP 


— Model of flares/CMEs proposed by Carmichael 1964; Stur- 




rock 1968; Hirayama 1974; Kopp & Pneuman 1976 


DEFT 


- Dynamical Evolution of Flux Tubes 


DEM 


— Differential Emission Measure (per flux) 


DOE 


— Department of Energy 


EM 


- Emission Measure (per flux) 


ESA 


- European Space Agency 


EUV 


— Extreme Ultraviolet 


Fe 


— Iron 


CDS 


- Gas-dynamic Shock 


GEM 


— Geospace Environmental Modeling 


HXR 


— Hard X-ray 


JES 


— Joined Equilibrium Solution 


JPL 


- Jet Propulsion Laboratory 


LASCO 


— Large Area Solar Coronagraph 


LHS 


— Left Hand Side 


MHD 


— Magnetohvdrodvnamics 

J ^^-^^j ^^^^^^^^^ 


MK 


— Mega- Kelvin 


NASA 


- National Aeronautics and Space Administration 


NSF 


— National Science Foundation 


RD 


- Rotational Discontinuity 


RHS 


— Right Hand Side 


SAD 


- Supra-arcade Downfiow 


SADL 


- Supra-arcade Downflowing Loop 


SMS 


- Slow-mode Shock 


SOHO 


- Solar & Heliospheric Observatory 


SXR 


- Soft X-ray 


SXT 


- Soft X-ray Telescope on Yohkoh 


TFT 


- Thin Flux Tube 


TRACE 


- Transition Region and Coronal Explorer 



